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Summary. — Neutrinos play a crucial role in the collapse and explosion of massive 
stars, governing the infall dynamics of the stellar core, triggering and fueling the 
explosion and driving the cooling and deleptonization of the newly formed neutron 
star. Due to their role neutrinos carry information from the heart of the explosion 
and, due to their weakly interacting nature, offer the only direct probe of the dynam¬ 
ics and thermodynamics at the center of a supernova. In this paper, we review the 
present status of modelling the neutrino physics and signal formation in collapsing 
and exploding stars. We assess the capability of current and planned large under¬ 
ground neutrino detectors to yield faithful information of the time and flavor depen¬ 
dent neutrino signal from a future Galactic supernova. We show how the observable 
neutrino burst would provide a benchmark for fundamental supernova physics with 
unprecedented richness of detail. Exploiting the treasure of the measured neutrino 
events requires a careful discrimination of source-generated properties from signal 
features that originate on the way to the detector. As for the latter, we discuss 
self-induced flavor conversions associated with neutrino-neutrino interactions that 
occur in the deepest stellar regions; matter effects that modify the pattern of flavor 
conversions in the dynamical stellar envelope; neutrino-oscillation signatures that 
result from structural features associated with the shock-wave propagation as well 
as turbulent mass motions in post-shock layers. Finally, we highlight our current 
understanding of the formation of the diffuse supernova neutrino background and 
we analyse the perspectives for a detection of this relic signal that integrates the 
contributions from all past core-collapse supernovae in the Universe. 

PACS 14.60. Pq, 14.60.Lm, 97.60.Bw- . 


1. — Introduction 

Core-collapse supernovae are among the most powerful sources of neutrinos in our 
Universe. During a supernova explosion, 99% of the emitted energy 10 53 erg) is re¬ 
leased by neutrinos and antineutrinos of all flavors, with energy of several MeV, which 
play the role of “astrophysical messengers,” escaping almost unimpeded from the super¬ 
nova core. The supernova neutrino flux has been extensively studied as a probe of both 
fundamental neutrino properties and core-collapse physics. Therefore, supernova neu¬ 
trinos represent a truly interdisciplinary research field at the interface between particle 
physics, nuclear physics and astrophysics. 

While Galactic supernovae are rare, existing or proposed large neutrino detectors will 
allow collection of a high-statistics neutrino signal from the next Galactic explosion. The 
supernova neutrino detection will be crucial to test the explosion mechanism and thus to 
compare current supernova models with direct empirical information from the supernova 
core. Originating from deep inside the core, neutrinos are affected by flavor conversions 
in the dense supernova matter on their way through the stellar mantle and envelope. 
Therefore, the neutrino fluxes reaching the detectors will carry intriguing signatures of 
oscillation effects in the deepest supernova regions, depending on the unknown neutrino 
mass hierarchy. In this sense, the dense supernova interior represents a unique laboratory 
to probe neutrino flavor mixing under high-density conditions. 

Matter effects in a supernova can be truly dramatic as neutrinos propagate through a 
dense turbulent environment. Furthermore in the deepest supernova regions, the neutrino 
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density is so high that it dominates the flavor evolution, producing a fascinating collective 
behavior associated with v-v interactions [1-3]. In the recent past, our description of self- 
induced neutrino oscillations has seen substantial progress. The seminal studies started 
almost a decade ago [4-9] stimulated a still-ongoing torrent of analytical and numerical 
works to clarify several aspects of this unusual flavor dynamics (see [10] for a review). 

It is thus evident that the physics potential of a supernova neutrino detection is 
enormous. Stellar collapse neutrinos were observed for the first (and so far only) time 
in 1987. Within a decade since the advent of large underground neutrino detectors, 
Nature was kind enough to provide a neutrino burst associated with the collapse of 
Sanduleak -69° 202 in the Large Magellanic Cloud, 51 kpc away. Two kiloton-scale water 
Cherenkov detectors, one in Japan, Kamiokande-II [11], and one in the United States, the 
Irvine-Michigan-Brookhaven (IMB) experiment [12], observed about 20 events within 13 
seconds, with timing consistent with the optical observation of the SN 1987A explosion. 
The positron spectra produced by inverse beta decays of supernova v e in the two detectors 
are represented in Fig. 1. Two smaller scintillator detectors, Baksan [14] and the Mont 
Blanc Liquid Scintillation Detector (LSD) [15] also reported several events ( 1 ). 

The detection of SN 1987A neutrinos has long been taken as a confirmation of the 
salient features of our physical comprehension of the core-collapse supernova phenomenon 
and of the associated neutrino emission [16]. This observation allowed us to put strong 
constraints on exotic neutrino properties (e.g., decays, neutrino charge) that would have 
altered the supernova neutrino emission. Most importantly, the total energy of v e and 
the inferred cooling time scale of a few seconds of the proto-neutron star put severe 
limits on non-standard cooling mechanisms associated with new particles emitted from 
the supernova core, notably right-handed neutrinos and axions (see [17] for a review). 
At the same time, with the poor statistics of the SN 1987A events, determining both 
supernova and neutrino parameters is impossible, although some hints could be obtained, 
e.g., the signal of SN 1987A was analysed in the light of neutrino oscillations [16,18,19]. 

The Galactic supernova rate can be estimated with different techniques (see, e.g., Ta¬ 
ble X in [20]). Typically one expects 1-3 core-collapse supernovae per century. However, 
except for SN 1987A in the Large Magellanic Cloud, no stellar collapse has been observed 
over more than 30 years of neutrino experiments. This absence of a signal allowed differ¬ 
ent neutrino experiments to place non-trivial upper bounds to the rate of collapses and 
failed supernovae [21-23], which confirm that a Galactic supernova explosion is a rare 
event. 

However, there are ~ 10 supernova explosions per second in the visible Universe. The 
cumulative emission of neutrinos from all the past core-collapse supernovae produces a 
cosmic background of (anti)neutrinos, the so-called Diffuse Supernova Neutrino Back¬ 
ground (DSNB), whose existence was predicted already before the observation of SN 
1987A [24]. Although weak, the DSNB is a “guaranteed” signal that can also probe 
physics different from a Galactic explosion, including processes which occur on cosmo¬ 
logical scales in time or space. The DSNB covers a wide range of physics, including 
the cosmic star formation rate, the stellar dynamics and fundamental neutrino proper¬ 
ties. Forecasts of the DSNB can be obtained using the neutrino spectra predicted by 
supernova simulations or the ones reconstructed from SN 1987A data [25], as shown in 
Fig. 2. The stringent observational upper limit on the DSNB flux, obtained by Super- 


( 1 ) The reported LSD burst is more controversially associated with the supernova, because the 
events were recorded some hours earlier than the others. 
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Fig. 1. - Positron spectra detected at Kamiokande-II (upper panel) and IMB (lower panel) 
in connection with SN 1987A. (Reprinted figure from [13]; copyright (2005) by the American 
Physical Society.) 


Kamiokande [26,27], is only a factor of ^ 2 higher than typical theoretical estimates [28]. 
This suggests an imminent detection of the DSNB in current and planned detectors. 

Nearly three decades after SN 1987A, we are eager for the next supernova neutrino 
signal. Meanwhile our understanding of supernova neutrinos has grown significantly. 
Remarkably, our modeling of the supernova neutrino emission and of flavor conversions 
in the stellar matter has experienced several breakthroughs. On the other side there is 
vivid experimental activity in low-energy neutrino astronomy, whose main goal is the 
detection of Galactic and diffuse supernova neutrino signals. Given the recent advances 
both on the theoretical and experimental fronts, in this review we aim to present state- 
of-the-art supernova neutrino physics and astrophysics. Our plan is as follows. 

In Section 2 we review the neutrino-emission phases of core-collapse supernovae, the 
characteristic signal properties expected during these phases, and our current understand¬ 
ing of the delayed neutrino-driven mechanism based on multi-dimensional modeling of 
stellar collapse and explosion. Besides reporting new and unexpected neutrino-emission 
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Fig. 2. - The DSNB detection spectra based on the neutrino spectra inferred from either the 
Kamiokande-II or IMB data sets alone or their combination, compared to a model (shaded shape) 
with canonical neutrino emission parameters. (Reprinted figure with permission from [25]; 
copyright (2007) by the American Physical Society.) 


features associated with the nonradial hydrodynamic flows during the accretion phase, 
we also present long-timescale proto-neutron star cooling calculations with state-of-the- 
art treatment of neutrino transport, in which —for the first time— convective transport 
effects were taken into account by a mixing-length description for spherically symmet¬ 
ric simulations. In Section 3 we discuss the detectability of (extra)galactic supernova 
neutrinos in current and planned detectors. We describe the physics potential of the dif¬ 
ferent detection techniques proposed to measure the neutrino burst and present the most 
relevant applications of a supernova detection to neutrino astronomy and astrophysics. 
Section 4 is devoted to the characterization of supernova neutrino flavor conversions. 
Self-induced effects associated with the neutrino-neutrino interactions in the deepest su¬ 
pernova regions and matter effects in the dynamical supernova environment are discussed. 
Finally, observable signatures of flavor conversions imprinted on the neutrino burst are 
presented. In Section 5 we focus on the DSNB. We discuss our current knowledge of 
the cosmic supernova rate and present estimates of the expected DSNB signal in large 
underground neutrino detectors by taking into account the effects of flavor conversions 
as well as delicate background issues, and including core-collapse and invisible supernova 
progenitors. Finally, conclusions and future perspectives are discussed in Section 6. 

We remark that in the current review we assume a standard three-neutrino oscillation 
framework. We neglect exotic neutrino properties such as a neutrino magnetic moment, 
neutrino decays or extra sterile neutrino states that would have a potential impact on the 
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supernova neutrino signal [17]. Also, aspects of neutrino-induced and neutrino-affected 
nucleosynthesis will only be mentioned in passing. This field is linked to a large diversity 
of questions that reach far beyond the scope of the topic of our review and which are 
addressed in Ref. [29]. 
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2. — Neutrino Signals from Stellar Collapse and Supernova Explosions 

Authors: H.-Th. Janka, R. Bollig, L. Hildepohl 

In a sequence of hydrostatic nuclear burning stages massive stars build up degenerate 
iron or oxygen-neon-magnesium cores, whose gravitational collapse terminates the stellar 
life. These events are the strongest cosmic sources of MeV neutrinos, comparable to 
neutron-star mergers, but roughly 1000 times more frequent. An intense burst of order 
10 58 neutrinos is released on a time scale of several seconds when a neutron star (NS) 
or black hole (BH) is born, possibly accompanied by a supernova (SN) explosion that 
expels several solar masses of stellar debris with velocities up to some ten percent of the 
speed of light. 

Neutrinos are the main agents for the transport of energy and lepton number during 
the infall of the stellar core and the formation of its compact relic. Therefore they 
play a decisive role during all stages of such an event. Electron neutrinos {y e ) are 
produced by electron captures on nuclei and on free protons and thus accelerate the 
initial implosion. Their continuous release drives the evolution from the lepton-rich post¬ 
collapse configuration to the final deleptonized NS, while neutrinos and antineutrinos 
of all flavors carry away the gravitational binding energy of the assembling remnant. 
Neutrino-energy deposition behind the stalled bounce shock can revive the shock and can 
thus initiate the SN explosion. The interaction of electron neutrinos and antineutrinos 
with the neutrino-heated outflow sets the neutron-to-proton ratio and hence determines 
the nucleosynthesis in the innermost SN ejecta. The intense neutrino radiation that 
escapes from the central regions can trigger nuclear spallations, whose free nucleons can 
seed interesting chemical element formation even in the outer stellar shells. Matter and 
neutrino-induced flavor oscillations inside the exploding star can affect these processes, 
and vacuum oscillations as well as the matter effects in the Earth contribute in shaping 
the neutrino signal that is detectable in experimental facilities on earth. 

In this Section we summarize the present status of the numerical modeling of stel¬ 
lar core collapse and explosion and of the associated neutrino emission by spherically 
symmetric and multi-dimensional simulations. The effects of neutrino flavor oscillations, 
which are usually ignored in the source modeling and only computed in a post-processing 
treatment, will be discussed in Section 4. 

2T. Neutrinos from supernovae and SN 1987A in retrospect. - The Nobel prize 
awarded detection of neutrinos from SN 1987A on February 23, 1987, sets a landmark 
for the beginning of extragalactic neutrino astronomy. It is the first direct empirical 
proof of the formation of a hot NS in the gravitational collapse of the core of a massive, 
evolved star. The overall properties of the measured neutrinos [11,12,14], i.e., their in¬ 
dividual particle energies, signal duration, and total signal energy, were compatible with 
theoretical predictions of the neutrino emission from the birth of NSs [30]. Although the 
limited statistics of only two dozen registered events did not allow for high significance, 
numerous relevant constraints on particle properties and fundamental physics could be 
derived on the basis of this neutrino observation in the months and years after SN 1987A 

[ 17 ]- 

The lack of signal statistics also appears as the most likely explanation of some puz¬ 
zling features of the neutrino detection, in particular the gap of roughly 7 s between ^2 s 
and ^9 s in the Kamiokande-II data, and the clear excess of events pointing away from the 
source. Moreover, the Kamiokande and IMB measurements are only marginally consis- 
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tent with each other. A joint analysis of both experiments also yields a best-fit value for 
the spectral temperature of the detected electron antineutrinos (p e ) that is considerably 
lower than expected for the time-integrated neutrino emission on grounds of the most 
detailed models existing at the time of SN 1987A. This tension is even enhanced given 
the neutrino mixing parameters indicated by the current phenomenology (cf. Sec. 4) and 
the high mean energies that had been predicted for the radiated heavy-lepton neutrinos 
(y x = z/ r , v T ) by those early models [19]. 

Meanwhile, however, this conflict has disappeared because state-of-the-art simula¬ 
tions of the Kelvin-Helmholtz cooling of nascent NSs yield significantly softer heavy- 
lepton neutrinos and a smaller difference of their spectrum compared to that of electron 
antineutrinos [31-34]. These trends are along the lines of model adjustments that had 
been identified in Ref. [19] as a necessary implication of the SN 1987A neutrino data if 
solar neutrino experiments would bear out a large mixing angle between v e and a heavy- 
lepton neutrino. The improvements in modern treatments of neutrino transport in SNe 
are connected to the introduction of nucleon-nucleon bremsstrahlung as the most impor¬ 
tant pair-production process [35,36], whose relevance for neutrino-spectra formation in 
NS cooling models had first been pointed out in Ref. [37]. The similarity of v e and v x 
spectra is further enhanced by the inclusion of energy transfers in neutrino-nucleon scat¬ 
terings [31,38]. A detailed reanalysis of the SN 1987A neutrino detections on the basis 
of signal predictions by modern hydro dynamical models of SN 1987A, also applying our 
current knowledge of SN neutrino oscillations, still needs to be performed. This would 
remove some of the uncertainties in the two-phase parametrizations used in Refs. [39-41] 
and would allow for a closer assessment of the remaining discrepancies between model 
predictions and SN 1987A neutrino measurements. 

One of the biggest unsolved riddles in connection to SN 1987A is a cluster of 5 neu¬ 
trinos of 7-11 MeV registered within 7 s by the LSD scintillator experiment in the Mont 
Blanc Laboratory roughly 4.5 hours before the Kamiokande II, IMB and Baksan events 
occurred [15]. The fact that none of these other experiments reported any significant de¬ 
tection at the time of the Mont Blanc measurement might be understood by the relatively 
low energies of these neutrinos and the higher detection thresholds of the experiments, 
whereas the small active mass of only 90 tons of the LSD detector has been used as an 
argument why this facility did not see any neutrinos when the other three experiments 
captured their events. However, the detection of five low-energy neutrinos by the small 
LSD mass requires a total energy in the neutrino burst that is several times higher than 
the typical energy release from NS formation [42]. No truly convincing scenario has been 
proposed as an explanation so far. In Refs. [43, 44] a two-stage collapse scenario was 
suggested in which the fragmentation of the collapsing stellar core leads to the forma¬ 
tion of a binary NS whose inspiral is driven by gravitational radiation and whose final 
merger causes the second neutrino burst. While core fragmentation seems possible in 
the presence of extreme rotation [45], the long delay to the second burst can hardly be 
attributed to gravitational-wave emission, because the collapse fragments are embedded 
in the dense, infalling gas of the SN progenitor and must be expected to dissipate the 
energy of their orbital motion by hydrodynamic effects (pressure waves and shocks) much 
faster. 

Another long-standing question concerns the possible formation of a BH in SN 1987A. 
Considering the fact that recent observations set the lower limit of the maximum mass of 
nonrotating, cold NSs to more than 2 M 0 [46,47], BH formation from the collapse of an 
18-20 Mq progenitor of SN 1987A [48] appears very unlikely from the theory perspec¬ 
tive. Neither the stellar iron core of such a star nor the fallback mass expected for an 
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explosion energy of more than 10 51 erg, which was deduced from SN 1987A observations, 
are sufficiently large to push the NS beyond the BH threshold after the emission of the 
detected neutrino burst. From the observational perspective there is also no real problem 
with the lack of clear evidence for a pulsar so far. The current luminosity of the ejecta 
cloud of this SN 10 36 ergs -1 ; [49,50]) is not compatible with the spin-down power of 
a bright pulsar like Crab (~4.5xlO 38 ergs -1 ). However, the ejecta emission is strongly 
affected by the shock interaction with the circumstellar medium, and the corresponding 
brightening is well able to cover the radiation of a thermally cooling NS like the compact 
object in Cassiopeia A with a present bolometric luminosity of about 7 x 10 33 ergs -1 
[51] and an initial luminosity of at most ^10 35 ergs -1 (e.g., [52,53]). Interestingly, very 
recent spectral and morphological analysis of the remnant of SN 1987A with ALMA and 
ATCA data seems to indicate the possible presence of a diffuse synchrotron source in 
the unshocked ejecta at a westward offset from the SN position. If caused by the energy 
release of a pulsar wind nebula it would set an upper limit to the pulsar spin-down power 
of ^10 35 ergs -1 [54]. The pulsed emission of such a spinning NS would so far have re¬ 
mained undetected because of its relative faintness and the obscuration of the scattering 
and absorbing stellar debris around the SN site. 

While the neutrinos from SN 1987A were able to yield basic confirmation of the 
theory of NS formation, their event statistics was too poor to allow for useful insights 
into the dynamics of the beginning explosion and thus into the still debated explosion 
mechanism of core-collapse SNe. A next Galactic SN, however, is likely to provide a 
high-statistics neutrino signal, depending on the availability of running experimental 
facilities. It will therefore carry valuable information about the properties of the neutrinos 
themselves (e.g., [55]), and it will, in particular, help us unravelling one of the most 
nagging problems of stellar astrophysics, namely why and how massive stars achieve to 
reverse their catastrophic collapse to the gigantic blast of the SN. But even if a massive 
star in our Milky Way ends its life without a brilliant explosion, which is a possibility that 
could happen in a fair fraction of up to several 10% of all stellar core collapses [56-59], the 
formation of a BH will be accompanied by a luminous outburst of neutrinos. In this case 
neutrinos (and possibly gravitational waves) might be the only messengers of the “silent” 
stellar death. The magnitude of the neutrino luminosity will be a measure of the rate 
of mass infall from the collapsing stellar core to the accreting NS, and large-amplitude 
luminosity variations will reflect the dynamics of the (unsuccessful) SN shock, caused 
by shock expansion or contraction episodes, which modulate the NS accretion rate. The 
final termination of the neutrino emission will mark the collapse of the transiently stable 
NS to a BH and thus might set constraints to the uncertain equation of state (EoS) at 
super nuclear densities. 

A high-statistics measurement of neutrinos from a future Galactic SN will therefore 
be of paramount importance not only for astrophysics and neutrino physics but for 
nuclear physics, too. Also the measurement of the diffuse supernova neutrino background 
(DSNB) as an integral signal of stellar core collapses in the past will bear a high potential 
for interesting astrophysics. It will be able to set limits on the total rate of stellar core¬ 
collapse events in the local Universe. If the rate of BH formation cases accounts for a 
considerable fraction of these events, they might leave characteristic fingerprints in the 
DSNB spectrum because of pronounced differences between the neutrino luminosities and 
spectra radiated by NS and BH formation. Improved predictions of the properties of the 
DSNB are an important task for future studies to be performed when more detailed stellar 
collapse and SN models with sophisticated neutrino transport have become available for 
large sets of progenitors. 


9 



Supernova Neutrinos: Production, Oscillations and Detection 


2*2. Numerical methods for supernova modeling in this work. - All simulations by 
the Garching group reported in this article, spherically symmetric (ID) as well as multi¬ 
dimensional, were performed with the Prometheus-Vertex SN code [60]. It utilizes a 
two-moment scheme for the transport of neutrinos and antineutrinos of all three flavors, 
which employs a Boltzmann closure for the variable Eddington factor and accounts for 
the full energy and velocity (to order (v/c)) dependence of the transport in the comoving 
frame of the fluid. For multi-dimensional problems a “ray-by-ray-plus” approximation is 
applied [60,61]. Gravity is generalized beyond the Newtonian description by an effective 
relativistic potential [60], adopting “Case A” of Ref. [62], and relativistic reshifting is 
taken into account in the neutrino transport. The code includes the whole set of neu¬ 
trino opacities as detailed in Ref. [60] with the improvements of Ref. [61]; see also the 
summaries in Refs. [63,64]. Recent upgrades of the code include( 2 ): 

(i) The possibility to treat all individual types of muon and tau neutrinos and antineu¬ 
trinos as separate species instead of clustering them together into one representative 
heavy-lepton neutrino [65]; 

(ii) the implementation of self-energy shifts of unbound neutrons and protons in their 
charged-current /^-reactions with neutrinos; 

(iii) an optional treatment of quasi-stationary convection in ID simulations by a mixing- 
length description (see Sect. 2*4.1). 

These upgrades will be mentioned when they are applied in the models described in this 
section. 

Although Prometheus-Vertex makes use of an approximate description of the 
effects of general relativistic gravity (the hydrodynamics solver is Newtonian), it permits 
simulations up to the onset of BH formation. Tests did not only show good compatibility 
with fully relativistic calculations during the whole cooling evolution of nascent NSs 
[66]. Also for proto-neutron stars that approach gravitational instability by continuous 
accretion, test calculations in comparison to cases in the literature [67] revealed an almost 
perfect match of the time scale to reach the collapse of the NS and very good agreement 
with respect to the neutrino-emission properties (with minor differences for v e and v e 
and maximal differences of 10-20% for v x \ [68]). 

2*3. Multi-dimensional phenomena in supernova cores and explosion mechanism. - 

2*3.1. The neutrino-driven mechanism. When nuclear-matter densities are reached at 
the center of a gravitationally imploding stellar core, a NS begins to assemble from the 
infalling matter. The stiffening of the EoS due to repulsive forces between the nucleons 
leads to an abrupt bounce that sends a strong hydrodynamical shock wave into the 
supersonically collapsing outer material of the stellar core. Only milliseconds after core 
bounce, however, the newly formed SN shock is weakened by iron photo-disintegration 
and a prompt burst of electron neutrinos, which drain energy from the postshock material 
and reduce the pressure behind the shock. Therefore the shock stagnates and becomes 
an accretion shock with mass infall instead of expansion characterizing the flow in the 
downstream region. Increasingly more sophisticated numerical models have consolidated 


( 2 ) The implementations of points (ii) and (iii) in the simulation code in an early, preliminary 
version were assisted by Bernhard Muller. 
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Fig. 3. - Top: Three post-bounce snapshots (125, 165, 210 ms p.b.) of a 3D simulation of 
an 11.2 Mq star, whose postshock accretion flow is characterized by convection. (Image from 
Ref. [102]; copyright (2014) by the American Astronomical Society.) Bottom: Three images 
(240, 249, 278 ms p.b.) from a 27 Mq 3D simulation, in which episodes of strong activity by 
the standing accretion shock instability (SASI) alternate with convection-dominated periods. 
(Image from Ref. [70]; copyright (2013) by the American Astronomical Society.) Surfaces of 
constant entropy are displayed in yellow and red; the SN shock is visible as a bluish, semi¬ 
transparent envelope. SASI sloshing or spiral motions show up by large-amplitude unipolar 
or dipolar deformations, whose orientation flips between the hemispheres on time scales of 
milliseconds. 


this nowadays generally accepted failure of the prompt bounce-shock mechanism (see 
[63] for a review). On a much longer time scale of tens to hundreds of milliseconds 
after core bounce, however, neutrinos leaving the nascent NS deposit energy in the gain 
layer between the so-called gain radius and the shock, mainly through charged-current 
reactions with free nucleons: 


(1) V e + n —> p + e , 

(2) v e + p —> n + e + . 

As time goes on, the conditions for this neutrino heating become continually more favor¬ 
able because the neutrinos are radiated with increasingly harder spectra as the neutri- 
nospheric temperature rises in the contracting and compressionally heated NS. 

If the energy transfer by neutrinos is strong enough, it can raise the postshock pres¬ 
sure to trigger the reacceleration of the SN blast wave and to thus initiate a successful 

explosion [69]. At later times after bounce, not only the rising spectral temperatures 

of the radiated neutrinos improve the conditions for this shock revival. Also the rate 
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at which matter of the collapsing stellar core falls into the shock gradually declines and 
with it the ram pressure of the infalling material. This also favors the rejuvenation of the 
stalled shock when the mass accretion rate drops below the critical value at which the 
thermal pressure created by neutrino heating can overcome the ram pressure of matter 
ahead of the shock. The explosion therefore sets in with a considerable time delay after 
the initial shock formation at the moment of core bounce. 

For a period of 30 years, Wilson’s “delayed neutrino-driven mechanism” has sur¬ 
vived as the standard paradigm for explaining how massive stars achieve to explode, 
despite the fact that numerical modeling repeatedly experienced setbacks in its efforts to 
demonstrate the viability of the mechanism and despite the nagging lack of an ultimately 
convincing theoretical or observational confirmation, both of which still nourish scepti¬ 
cism and criticism. However, in the absence of rapid rotation and very strong magnetic 
fields, whose relevance for ordinary SNe is disfavored by current stellar evolution models 
[71], neutrinos are the most efficient way to extract energy from the rich reservoir of the 
hot, newly formed NS and to transfer parts of this energy to the overlying stellar shells 
in order to reverse their infall to an explosion. In fact, neutrinos carry away the huge 
gravitational binding energy of the compact remnant, several 10 53 erg, and less than one 
percent of this energy is well sufficient to account for the observed power of a typical SN. 

The basic functioning of the neutrino-driven explosion mechanism can be understood 
as a global runaway instability of the postshock accretion layer caused by the neutrino- 
energy deposition. According to the “critical luminosity condition” [72] the explosion is 
launched when the neutrino luminosity exceeds a critical threshold value that depends 
on the mass-accretion rate of the stalled shock. This concept is now widely accepted and 
basically consistent with many parametric numerical models (e.g., [73-77]) and analytic 
studies (e.g., [78,79]), although many details are still disputed, e.g. which parameters 
define the critical condition most accurately and can be used as most reliable indicator 
of the threshold for runaway (e.g., [76,79,80]). 

2*3.2. Importance of non-radial hydrodynamic instabilities. Hydrodynamic simulations 
with an increasingly more sophisticated description of the crucial neutrino physics and 
transport, however, have shown that a state-of-the-art treatment of the microphysics 
(in particular of the neutrino interactions and EoS) does not allow for explosions in 
spherical symmetry (i.e., in one dimension; ID). An exception to this are the lowest-mass 
progenitors of SNe, i.e., ~9-10 M© stars with oxygen-neon-magnesium cores [34,81,82] or 
small iron cores [83], which are surrounded by extremely dilute overlying shells with very 
low gravitational binding energy. But even in these cases the neutrino-heated postshock 
layer is convectively unstable and high-entropy, neutrino-heated plasma becomes buoyant 
and rises in Rayleigh-Taylor plumes (see Fig. 3, upper row). These multi-dimensional 
flows assist neutrinos in triggering the explosion earlier, accelerate the expansion of the 
SN shock, enhance the explosion energy, create asymmetries in the explosion ejecta, and 
modify the conditions for element formation in the neutrino-processed gas as well as in 
the layers of explosive nucleosynthesis [83,84]. 

For progenitor stars that are more massive than about 10 M 0 , the support by multi¬ 
dimensional hydrodynamic instabilities is indispensable to reach the critical condition for 
a neutrino-driven explosion. There is a variety of effects by which non-radial flows in the 
postshock layer can aid the onset of the SN explosion [73,85,86]. Non-radial convective 
flows stretch the residence time of matter in the heating layer and thus enhance the 
energy transfer by neutrinos (e.g., [32,74,87]). The buoyant rise of neutrino-heated 
matter also reduces the energy loss by the re-emission of neutrinos (mostly through the 
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Fig. 4. - Evolution of the convective region in the newly formed NS of a collapsing and exploding 
27 Mq star [68] using the EoS of Lattimer &; Swesty [90] with a nuclear incompressibility modulus 
of K — 220 MeV. The NS has a baryonic (final gravitational) mass of 1.776 (1.592) Mq. Time 
is normalized to bounce. The upper panels show a ID simulation without convection, the lower 
panels with a mixing-length treatment of convection. The left column displays the evolution of 
the Brunt-Vaisala frequency, positive values of which indicate instability for Ledoux convection. 
In the model that includes convective energy and lepton-number transport, the convectively 
unstable regions (red in the upper left plot) become convectively neutral (whitish in the lower 
left plot). The right panels display the temperature evolution. The cooling (and deleptonization) 
of the nascent NS is considerably faster when convective fluxes accelerate the neutrino loss. 



inverse reactions of Eqs. (1,2)). Moreover, the expanding plumes push the shock to 
larger radii and thus increase the mass and the volume of the neutrino-heated layer. 
With the larger optical depth of the gain layer the efficiency of neutrino heating grows 
even more, creating ideal conditions for a runaway. Some authors have attributed the 
combination of these explosion-assisting effects, which are associated with the non-radial 
mass motions in the postshock layer, to the generation of turbulent pressure behind 
the shock [88,89]. The corresponding turbulent energy scales with the neutrino-energy 
deposition rate [36,80], and the turbulent pressure in addition to the gas pressure has 
been shown to effectively reduce the critical neutrino luminosity for shock revival [80]. 
Although this conceptual picture seems to be able to capture some basic features found in 
multi-dimensional explosion models, the actual, highly time-dependent gas dynamics in 
collapsing stellar cores is probably too complex to be compatible in all phases and in all 
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Fig. 5. - Left: Post-bounce luminosities of v e , L e , and v x as measured by a distant observer for 
a 3D collapse simulation of a 27 Mq star. The observer position is close to the plane of a spiral 
SASI shock-oscillation mode between ~120ms and ~250ms. For the displayed luminosities the 
neutrino emission was integrated over the hemisphere facing the observer with limb-darkening 
effects taken into account as described in the appendix of Ref. [91], from where the figure was 
taken. (Image copyright (2014) by the American Physical Society.) The SASI modulations of 
the mass-accretion of the NS imprint large-amplitude, near-sinusoidal variations on the neutrino 
emission. There is a second phase of SASI activity starting at ~400ms, whose orientation is 
inclined to the direction of the observer. Therefore the plot does not reflect the full amplitude 
of these later modulations. Right: Relative amplitude of the v e signal variations on a “sky-plot” 
of all observer directions during the first SASI phase (120-250 ms p.b.) in the 27 M© model. 
(Image from Ref. [92]; copyright (2013) by the American Physical Society.) 


aspects with the properties of fully developed turbulent flows at steady-state conditions. 

Theoretical models therefore led to the conclusion that the neutrino-driven SN mech¬ 
anism is a generically multi-dimensional phenomenon. This important insight is in line 
with observations of SN 1987A and of other relatively nearby and well studied SNe, and 
it is further supported by the morphological properties seen in young SN remnants. All of 
these events show large-scale asymmetries, non-spherical deformation, and extended ra¬ 
dial mixing of the chemical elements ejected during the explosion, none of which could be 
understood if the onion-shell stratification of the pre-collapse star was preserved during 
the SN blast. 

The existence of convectively unstable regions in the SN core was recognized soon after 
the first detailed hydrodynamic simulations had revealed the evolution and structure of 
the newly formed NS and its surrounding layers. Besides the negative entropy gradient 
that is built up by neutrino heating behind the stalled shock, and which is unstable 
to convective overturn [93] setting in typically 80-100 ms after bounce, the decelerating 
and weakening bounce shock also leaves a negative entropy profile behind, which can 
decay in an early post-bounce phase of prompt convection lasting some ten milliseconds 
[94-97]. Moreover, in the hot and still lepton-rich proto-neutron star the combination 
of a negative entropy gradient and an unstable lepton-number gradient drives Ledoux 
convection in a thick shell interior to the neutrinosphere [94,98]. Over a period of seconds 
this convective activity penetrates deeper and deeper towards the center and accelerates 
the cooling and deleptonization of the nascent NS (Fig. 4; [30,99-101]). 

In addition to being stirred by convective mass motions, the accretion shock and the 
whole postshock accretion flow are unstable to global, non-radial deformation modes, 
which do not allow the accretion shock to remain spherical and at a stationary radius. 
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Fig. 6. - Time evolution of the lepton-number ( v e minus V e ) flux density normalized by the 
average value over all directions in a 3D core-collapse simulation of an 11.2 Mq progenitor. The 
panels show “all-sky” images for the indicated post-bounce times. One can see the emergence 
of a clear dipole pattern with strong excess of the v e emission in one hemisphere and a reduced 
v e emission or even excess flux of v e in the opposite hemisphere. This phenomenon has been 
found in all 3D simulations of the Garching group and was termed LESA for lepton-emission 
self-sustained asymmetry. The black dot marks the maximum of the dipole, the cross the anti¬ 
direction. The gray line indicates the path described by a slow drift of the dipole direction. 
(Image from Ref. [102]; copyright (2014) by the American Astronomical Society.) 


The possibility of such a non-convective instability in the SN core was first discussed 
in Ref. [103] and was termed “standing accretion shock instability” (SASI). It leads to 
the oscillatory growth of an initially small seed perturbation with the largest growth 
rates being found for the spherical harmonics modes of lowest order, i.e., the dipolar and 
quadrupolar modes [104-106]. This instability leads to large-amplitude shock sloshing 
and spiral motions [107] (see Fig. 3, bottom row). The underlying growth mechanism 
for this instability is a so-called advective-acoustic cycle in the accretion flow between 
stagnant shock and proto-neutron star [108-110] and can also be studied in an inexpensive 
water experiment that shares basic features with the accretion flow in a SN core ([111]; 
for a nice review with focus on this instability and a more complete collection of relevant 
references, see Ref. [112]). 

Both instabilities, postshock convection and the SASI, can provide crucial support 
to the onset of neutrino-driven explosions. A growing number of groups with energy- 
dependent neutrino-transport schemes of different levels of sophistication (and quite a 
variety of differences in other important numerical and physical aspects of the hydro- 
dynamical modeling) have meanwhile obtained successful neutrino-driven explosions of 
progenitors above ^10 M 0 in self-consistent, first-principle simulations in two dimen¬ 
sions (2D), i.e., in simulations that were artifically constrained to axisymmetry (e.g., 
[32,64,87,113-119]). This suggests the basic viability of the neutrino-driven mechanism, 
although there is still no unanimous agreement [120]. Simulations in 2D, however, are 
problematic because of the artifically imposed axial symmetry, which attributes a toroidal 
geometry to all structures, and because of the inverse direction of the turbulent energy 
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Fig. 7. - Post-bounce evolution of the lepton-number emission in 3D core-collapse simulations 
of 11.2 M© (left), 20 M© (middle) and 27M© (right) progenitors. The upper panels show the 
monopole, i.e., the total lepton-number flux (red curve), and the dipole component (blue curve). 
The bottom panels display the polar angles 0 and f> defining the dipole direction (i.e., the direc¬ 
tion of the maximum excess of the v e emission relative to the v e emission) in the polar coordinate 
grid of the star (north and south pole directions are indicated by “N” and “S”, respectively). 
The monopole and dipole amplitudes A mon and Adi p are normalized such that the lepton-number 
flux is given by A mon + Adi p cos$ in coordinates aligned with the dipole direction, if the flux 
distribution contains only monopole and dipole terms. (Image from Ref. [102]; copyright (2014) 
by the American Astronomical Society.) 


cascade (from small to large scales) compared to the realistic three-dimensional (3D) 
situation. It is therefore indispensable to confirm the functioning of the mechanism by 
3D simulations, which have recently become possible because of the increasing power of 
modern supercomputers. 

2*3.3. Status of explosion modeling in 3D. The first self-consistent 3D simulations of 
stellar core collapse and explosion were performed with a relatively schematic neutrino 
treatment by grey diffusion [121,122]. More recently, first-principle 3D simulations with 
energy-dependent neutrino transport (however applying transport schemes with largely 
differing degrees of sophistication) have also obtained successful explosions, although the 
revival of the stalled shock happens somewhat later than in the corresponding 2D models 
[7,123,124,126]. The longer delay of the onset of the explosion may be problematic 
because later explosions tend to be weaker and might be incompatible with observed 
SN energies. However, the detailed effects in 3D and the differences between 2D and 
3D models are still a matter of intense research, and the current models do not yet 
provide final answers, in particular concerning the role of the SASI in 3D (cf. [?]) and 
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concerning the energetics of the explosions (cf. [?]). Moreover, more studies are needed to 
investigate the question of numerical convergence in describing the potentially turbulent 
postshock flow (e.g., [129-131]) and to explore the influence of so far poorly understood 
perturbations and non-sphericities in the convective silicon and oxygen burning layers of 
the pre-collapse star (e.g., [80,132,133]). 

Also the details of the microphysics in the newly formed NS remain a matter of 
concern. Interestingly, in Ref. [?] it was found that a modest reduction of the neutrino- 
nucleon neutral-current scattering opacity by just 10-20%, e.g. due to possible effects of 
strange-quark spins in the axial-vector structure factor of the scattering cross section, is 
sufficient to convert a 3D simulation of a 27 M 0 by the Garching group from failure to 
a successful explosion. This result clearly demonstrates the proximity of current state- 
of-the-art 3D simulations of stellar core collapse to explosions, and it emphasizes the 
sensitivity of the outcome of these simulations to the detailed input in the neutrino- 
opacity sector. A similar sensitivity exists with respect to the NS EoS, because it has 
been known for a while already that “softer” EoSs are favorable for the possibility of 
neutrino-driven explosions [87,114,117]. Both effects, reduced neutral-current neutrino- 
nucleon scattering rates as well as a soft EoS, lead to a faster contraction of the newly 
formed, hot NS in the first ^0.5 s after bounce( 3 ). In the case of a reduced neutrino¬ 
scattering opacity the faster contraction is a consequence of the enhanced emission of 
heavy-lepton neutrinos. A faster NS contraction causes more rapid compressional heating 
of the neutrinospheric layer and thus allows for higher accretion luminosities and a steeper 
rise of the average energies of the radiated neutrinos with time. Both higher luminosities 
and higher mean energies enhance the neutrino heating behind the shock and are therefore 
supportive for neutrino-driven explosions. 

Non-radial hydrodynamic instabilities in the SN core are not only of crucial impor¬ 
tance for the neutrino-driven mechanism; they also impose asymmetries on the beginning 
explosion and therefore set the conditions for NS acceleration through the “gravitational 
tug-boat mechanism”, which can yield NS kick velocities in agreement with those ob¬ 
served for young pulsars (e.g., [134,135]). Moreover, they seed the observable ejecta 
asymmetries that develop by the interaction of the primary non-sphericities with sec¬ 
ondary mixing instabilities that grow during the first day of the SN blast (e.g., [136]). 
The violent hydrodynamic flows associated with convective overturn and the SASI in the 
postshock layer also produce characteristic imprints on the neutrino signal emitted dur¬ 
ing the shock-stagnation phase after bounce. While convection leads to small-amplitude, 
high-frequency (> 100 Hz) fluctuations of the neutrino luminosities and mean energies, 
SASI sloshing and spiral motions create quasi-periodic variations of much bigger ampli¬ 
tudes (up to 10-20 percent in the neutrino luminosities, cf. the left panel of Fig. 5, and 
1-2 MeV in the mean neutrino energies; [91,92]) and with typical frequences of < 100 Hz, 
because the large-scale shock expansion and contraction phases associated with the SASI 
modulate the mass-accretion flow towards the nascent NS massively. During episodes 
of enhanced accretion the emission of neutrinos is boosted by additional accretion lumi- 


( 3 ) Note that here a “softer” EoS is defined by the faster contraction of the newly formed NS 
during the early post-bounce accretion phase and not by a lower maximum mass of a cold, non¬ 
rotating NS, which is often considered as the criterion for discriminating soft from stiff EoSs. It 
is important to realize that neither the radius nor the maximum mass of a cold NS are of direct 
relevance for the properties of the hot, accreting remnant immediately after core bounce and 
therefore for the conditions that influence the onset of the SN explosion during this phase. 
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nosity, and the neutrinos escape with higher mean energies because of the compressional 
heating of the accretion layer. 

Measurements of neutrinos from a future Galactic SN with Cherenkov telescopes will 
well be able to detect these signal modulations [91,92,137,138] (see Sect. 3’3.3 for more 
information). The observation of SASI modulations for SN neutrinos would provide 
a very important confirmation of our present, purely theoretical picture of the shock 
dynamics in the SN core. The presence of one or more SASI episodes is more probable 
in cases of more massive progenitor stars (maybe with masses beyond ~15 M 0 ), where 
the shock expansion comes to a longer halt and even shock contraction can occur before 
the onset of runaway. However, the possibility to see the corresponding neutrino signal 
modulations will depend on the viewing angle relative to the main direction or the plane of 
the SASI sloshing or spiral motions. Observers with positions close to the SASI sloshing 
axis or near the main plane of the SASI spiral mode will receive bigger modulation 
amplitudes (Fig. 5, right panel; see also Sect. 3‘3.3) and will have the better chance to 
diagnose them in the detection. 

2*3.4. LESA: A dipolar neutrino-emission asymmetry as new phenomenon. Another 
interesting and novel phenomenon was recently discovered in the neutrino emission of 
the first 3D stellar core-collapse models with state-of-the-art multi-group three-flavor 
neutrino transport. Instead of persistently picturing the high-order multipole pattern 
of buoyant plumes and sinking downdrafts in the convection zones in the interior of the 
proto-neutron star and in the neutrino-heated postshock layer, the neutrino emission 
develops a strong dipolar asymmetry on a time scale of about 100-150 ms after bounce 
(Fig. 6; [102]). The observable hemispheric luminosity difference of electron neutrinos, 
z/ e , and electron anti-neutrinos, z/ e , individually can reach up to ^20% and even higher 
for the lepton-number (y e minus P e ) luminosity (Fig. 6), whereas it is only of order 1- 
2% for the sum of v e and v e and for heavy-lepton neutrinos, v x . This astonishing and 
completely unexpected phenomenon was found in all 3D SN simulations conducted by the 
Garching group so far, i.e., in a convection-dominated 11.2 Mq model as well as in SASI- 
dominated 20 and 27 M 0 cases (see Fig. 7) and also in a more recently published model 
of a 9.6 M q progenitor, in which the explosion starts on a short post-bounce time scale of 
only ^130 ms [83]. The development of a strong lepton-emission dipole therefore seems 
to be a generic instability of the neutrino transport in the convectively stirred proto¬ 
neutron star. The exact growth conditions of this instability are not yet understood, 
but the main contribution to the emission dipole builds up in the convective layer well 
inside of the neutrinosphere, while the outer accretion develops a dipolar asymmetry, 
too, and enhances the hemispheric difference of the lepton-number emission. In fact, 
the accretion asymmetry might trigger the growth of the dipole and seems to stabilize 
its existence over periods of at least hundreds of milliseconds. The phenomenon was 
therefore named “lepton-emission self-sustained asymmetry” or LESA (see Ref. [102] for 
a detailed discussion). 

After its steep growth phase until roughly 200 ms post bounce, the amplitude of the 
LESA lepton-number emission dipole can become even larger than the corresponding 
monopole (Fig. 7). The dipole direction is amazingly stable if one compares its slow 
drifting with the much shorter life time of convective cells in the inner and outer con¬ 
vection regions, which collapse and are regenerated on a time scale of typically ~10ms 
only. SASI and LESA are two different phenomena, which can be well distinguished also 
with respect to their effects on the neutrino emission [91]. SASI modulates the v e and 
v e emission coherently, whereas LESA amplitudes of the v e and v e fluxes have opposite 
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signs. The LESA dipole direction and SASI mass motions can have arbitrary orientations 
relative to each other, for example the LESA dipole and normal vector of the plane of 
SASI spiral mode are roughly parallel during the first SASI episode in the 20 M 0 model 
and nearly perpendicular in the 27 M 0 case. The LESA amplitude and dipole orientation 
can reflect the SASI-imposed modulations of the neutrino emission (see Fig. 7), but the 
presence of the LESA is not overruled by the SASI, and the LESA dipole direction is 
also remarkably stable during the phases of violent SASI sloshing and spiralling activity. 

The LESA phenomenon is a new type of instability, which occurs as a combined 
effect of multi-dimensional hydrodynamics and neutrino transport. Presently, however, 
it cannot be rigorously excluded that this stunning phenomenon is a numerical artifact, 
e.g. connected to the “ray-by-ray-plus” approximation used for the neutrino transport in 
the multi-dimensional simulations [120]. It is also possible that the growth of the LESA 
dipole is fostered by some special (but not necessarily realistic) aspects of the employed 
microphysics, in particular of the neutrino opacities. Nevertheless, even in these cases a 
lepton-number emission dipole would still be an interesting result of the 3D simulations. 

If happening in nature, LESA would have far-reaching observable consequences for 
SNe. It would not only imply that the detectable signal of electron antineutrinos is a 
function of the viewing direction during the LESA-active phase of the SN evolution. De¬ 
pending on the duration of the asymmetry, a neutrino-emission dipole (as sum over all 
neutrino species) with an amplitude of about one percent could also lead to a NS kick 
of tens of kilometers per second. Such velocities, however, are far too low to account for 
the typical space velocities observed for young pulsars, and they are subdominant com¬ 
pared to the velocities that can be achieved by the “gravitational tug-boat mechanism” 
associated with asymmetries of the mass ejection in SN explosions ([134,135] and refer¬ 
ences therein). Also the nucleosynthesis conditions in the innermost, neutrino-processed 
ejecta, whose neutron-to-proton ratio is set by the reactions of Eqs. (1) and (2), will 
depend on the direction, with proton excess being present on the side of the stronger 
v e emission and neutron excess in the opposite hemisphere. 3D models that track the 
explosion considerably beyond the onset of shock revival are needed to study the exact 
implications of these possibilities for a varity of progenitors. 

2‘4. Neutrino signals from proto-neutron star cooling. - Full 3D modeling of SN ex¬ 
plosions has only just begun and many aspects of the neutrino-driven mechanism are still 
heavily debated and not generally agreed on (see Sect. 2*3). In particular, it is presently 
still uncertain which progenitor stars explode, when exactly the explosions set in, which 
energies they develop, and which stars collapse to give birth to BHs. Parametric ex¬ 
plosion simulations are therefore useful to explore the landscape of possibilities and to 
predict progenitor dependent variations of the neutrino signal as input to the investiga¬ 
tion of neutrino oscillations in SNe (Sect. 4), to calculations of the DSNB (Sect. 5), and 
to studies of neutrino detection in connection with future Galactic SNe (Sect. 3). 

Spherically symmetric, i.e. one-dimensional (ID), simulations are a very efficient way 
to achieve this goal, in particular if the long-time cooling evolution of the newly formed 
NS over seconds to tens of seconds shall be followed, too. However, ID simulations do not 
explode in a self-acting way (except for the lightest progenitor stars of core-collapse SNe, 
cf. Sect. 2‘3) and therefore the explosions have to be triggered artificially. Moreover, ID 
models lack in the generically multi-dimensional effects in the neutrino emission associ¬ 
ated with convection, the SASI and the LESA during the accretion phase that precedes 
the runaway acceleration of the shock. Also the transition from the accretion phase to 
the Kelvin-Helmholtz cooling phase of the remnant cannot be reliably represented by ID 


19 



Supernova Neutrinos: Production, Oscillations and Detection 



Fig. 8. - Post-bounce evolution of the NS radius (defined by a density of 10 11 gem -3 ; top panels ), 
luminosities of z/ e , z'e, and (a single representative of) v x (given in units of Bs -1 = 10 51 ergs -1 ; 
middle panels ) and corresponding mean energies (defined as the ratio of energy flux to number 
flux; bottom panels) in the laboratory frame for the proto-neutron star formed during the SN 
explosion of a 9.6 Mq progenitor star. The baryonic (final gravitational) mass of the compact 
remnant is ~1.363 (1.252)M©. The left panels show a comparison of ID models with mixing- 
length treatment of NS convection (SFHo-z9.6co-z/-lD; thin black lines) and without (SFHo- 
z9.6o-lD; thick blue lines). The right panels demonstrate the excellent agreement between the 
ID model with convection and a corresponding 2D simulation, in which the convection is treated 
hydrodynamically. Both models explode self-consistently at nearly the same time in the ID and 
2D cases, which offers ideal conditions for a clean and accurate comparison. 


models, because such models are unable to describe the simultaneous presence of shock 
expansion and accretion that continues even after the revival of the stalled SN shock. 
This post-runaway accretion phase might last for hundreds of milliseconds [?, 87] and 
will not only exhibit an enhanced level of v e and v e emission due to the persistent accre¬ 
tion contribution to the neutrino luminosity; it may also show short, spike-like eruptions 
caused by bursts of v e and v e production when particularly massive accretion downflows 
hit the proto-neutron star surface and dissipate their kinetic energy in shock-heated 
plasma [116]. 

2*4.1. Mixing-length treatment of proto-neutron star convection. After accretion has 
ended, the proto-neutron star is left behind as a quasi-spherical object (if rotation does 
not play an important role). However, as discussed already in Sect. 2*3, a convective 
layer eats deeper and deeper into the NS and accelerates the transport of energy and 
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Fig. 9. - Neutrino-signal evolution for the ~1.776 (~1.592)M© baryonic (final gravitational) 
mass NS formed in the explosion of a 27 M© progenitor. The red lines correspond to a ID 
simulation without convection, the black lines to a ID model with a mixing-length treatment of 
proto-neutron star convection (indicated by the letter “c” in the model name). The left column 
shows the phase of the z/ e -burst at shock-breakout, the middle column displays the accretion 
phase, the right column the Kelvin-Helmholtz cooling phase of the proto-neutron star. The 
top panels display the neutrino luminosities for z/ e , v e and a single species of v x (given in units 
of Bs -1 = 10 51 ergs -1 ), the middle panels show the radiated mean energies (ratios of energy 
fluxes to number fluxes), and the bottom panels the shape parameters a of the neutrino spectra, 
all given as functions of time after bounce and as measured by an observer in the laboratory 
frame. Note that the explosions of the displayed ID models were artificially triggered at ~0.5s 
p.b. (in contrast to the explosions of the low-mass progenitor of Fig. 8, which developed self- 
consistently). The initiation of the explosion was done slightly differently in both models, which 
explains the differences that are visible between ~0.5s and ~0.8s during the transition from 
the accretion to the cooling phase. 


lepton number compared to neutrino diffusion on its own. In contrast to the highly time- 
dependent and aspherical mass flows associated with convective overturn and the SASI 
in the neutrino-heating layer, the quasi-stationary convection inside of the proto-neutron 
star offers the possibility to describe the effects of convective energy and lepton-number 
transport in ID simulations by a mixing-length approximation (e.g., [139]), because the 
time scale of hydrostatic equilibration is much smaller than the time scale of convective 
transport, which again is much smaller than the time scale of evolutionary changes of 
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the NS. For the convective lepton-number and energy fluxes one can thus write [68]: 


( 3 ) 

( 4 ) 


F le P —cm A ^^ lep 
L conv H u c / 'c ^ 5 

= a. 


where r is the radius, p the baryonic-mass density, P the pressure and e the specific inter¬ 
nal energy. In the neutrino-trapping regime at densities above the neutrinospheres, these 
thermodynamic quantities include the contributions from neutrinos. Correspondingly, 
liep is the total electron-lepton number including electron neutrinos and antineutrinos. 
The mixing length is coupled to the pressure scale height, A c = £P|dP/dr| -1 , where 
£ is a dimensionless free parameter of order unity. The convective velocity is deter¬ 
mined as v c = \J2gAp p~ 1 X c with g > 0 being the local gravitational acceleration, and 
A p = A c CLedoux is the density contrast between surrounding medium and convective fluid 
elements after travelling a distance A c when Cbedoux is the Ledoux-convection criterion 
given by 
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Here, s is the entropy per baryon, c s = (9P/3p) s? y lep is the adiabatic sound speed, and 
the thermodynamic T-indices are given by the derivatives 
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Instability to Ledoux convection requires CLedoux > 0, in which case the Brunt-Vaisala 
frequency co>bv is real and co>bv = sign (CLedoux 1 CLedoux | > 0 denotes the growth 

rate of convective perturbations in the linear regime. In the numerical simulations a 
value of £ = 1 is used, but tests showed no noticeable changes for values between 1 and 

10 . 

Figure 8 testifies the good agreement that can be achieved between ID simulations 
utilizing this mixing-length treatment of convection and corresponding 2D calculations. 
Results of a proto-neutron star formed in the collapse of a 9.6 M 0 progenitor are dis¬ 
played. The compact remnant assembles a baryonic mass of about 1.363 M 0 , corre¬ 
sponding to a final gravitational mass of ^1.252 M 0 for the nuclear EoS used in the 
simulation. Proto-neutron star convection sets in typically 30-50 ms after bounce, which 
is the time when differences between convective and non-convective models become vis¬ 
ible first. The NS radius shows the usual increase associated with interior convection, 
and the radiated neutrino luminosities and mean energies reveal the well known time and 
neutrino-species dependent relative changes compared to the ID, non-convective case (for 
a detailed discussion of the pre-explosion effects, see Ref. [32]; for post-explosion differ¬ 
ences, see Refs. [30,100,101]). The neutrino luminosities during the displayed phase of 
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the first 1.3 s after bounce are enhanced by convective energy and lepton-number trans¬ 
port by up to about 20%, and the mean neutrino energies are initially reduced and at 
later times increased by up to nearly 2 MeV with the largest effects for heavy-lepton neu¬ 
trinos. The 2D model and the ID simulation with convection exhibit basically identical 
time evolutions of the NS radius and neutrino parameters (right panels of Fig. 8). 

2*4.2. Neutrino-emission phases and characteristics. The effects of post-bounce ac¬ 
cretion on the nascent NS are not very prominent in the neutrino signal emitted from 
SNe of low-mass progenitors like the 9.6 Mq star of Fig. 8. In these events the shock 
moves outward essentially in a continuous expansion and the explosion sets in early after 
bounce. In the 9.6 Mq star this happens already at ^130 ms in the convective models 
and at roughly 300 ms p.b. in the non-convective case [83]. The accretion phase is there¬ 
fore very short in such a low-mass progenitor, and the mass delivered on the accreting 
NS during this time interval is relatively small. 

The influence of longer and stronger accretion and of convection in the proto-neutron 
star interior on the neutrino emission of a more massive compact remnant can be seen 
in Fig. 9. The NS there is born in the collapse of a 27 Mq progenitor and has a baryonic 
mass of 1.776 Mq and a final gravitational mass of ^1.592 Mq. Note that in contrast to 
the 9.6 Mq model the explosion of the 27 M 0 star cannot be obtained self-consistently 
in ID but requires artificial initiation. This was done at 0.5 s after bounce in both of the 
models shown in Fig. 9 by applying slightly different methods. This is the reason for the 
differences in the neutrino emission between 0.5 s and ^0.8 s after bounce. The features 
seen during this time and the differences between both models should not be interpreted 
by physics. 

Figure 9 displays the neutrino signal properties for the three main phases of neutrino 
production after core bounce, which are 

(a) the high-luminosity v e burst at shock breakout from the neutrino-trapping regime, 

(b) the subsequent accretion phase of the stalled shock until the beginning of the ex¬ 
plosion, 

(c) the long-time Kelvin-Helmholtz cooling of the newly formed NS on its way to a 
final, deleptonized and cold compact remnant. 

The neutronization burst radiates away the electron neutrinos that are abundantly pro¬ 
duced by rapid electron captures on protons when the hot postshock matter becomes 
transparent to neutrinos. It exhibits universal properties concerning its width (several 
milliseconds) and height (around 4 x 10 53 ergs -1 = 400bethes -1 = 400Bs -1 ), which 
depend only weakly on the progenitor star and uncertain nuclear physics in the SN core 
([140]; cf. Fig. 37 in Sec. 4). Also the faster rise of the v x luminosities compared to the 
slower increase of the v e emission in the first 10-20 ms after bounce is a generic feature of 
state-of-the-art neutrino transport simulations. It can be understood by the high electron 
and v e degeneracy in the postshock material right after shock breakout, which suppresses 
the rapid production of v e . These neutrino-emission features, which are predicted to be 
independent of the progenitor star, offer interesting possibilities for neutrino diagnostics 
and neutrino astronomy (see Sec. 4). Differences of the neutrino emission from different 
collapsing stars grow as time after core bounce progresses. This is a consequence of the 
gradual infall of the stellar shells that surround the initially collapsing inner core and 
whose density profiles vary greatly between progenitors of different masses. 
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In contrast to the 9.6 M 0 model of Fig. 8, the 27 M e star shows evidence for prompt 
post-bounce convection, which causes small differences in the neutrino luminosities and 
energies even during the first ^50 ms after bounce. Towards the end of this phase the 
convective activity interior to the neutrinosphere begins to develop. Since the explosion 
in ID simulations of such a massive star has to be triggered artificially, the duration 
of the accretion phase is uncertain (in fact, it is prescribed by the modeler), and the 
transition from the accretion to the cooling phase is not reliable. Moreover, neither the 
accretion phase nor the transition to the cooling phase can exhibit special features of 
the neutrino emission associated with multi-dimensional flows in the accretion region 
between the stalled shock and the NS (e.g., convective overturn, SASI, or sporadic, post¬ 
explosion accretion episodes), which were discussed in Sect. 2‘3 and at the beginning of 
this section. Nevertheless, the neutrino signal of ID models still reflects some generic 
properties of the accretion phase. These include the following aspects: 

• The luminosities of heavy-lepton neutrinos can rather well be expressed in terms 
of a Stefan-Boltzmann relation, 

(7) L Vx = 4w<jx7 v lfRl, 

where o v = 4.751 x 10 35 ergMeV -4 cm -2 s -1 . The luminosity, L Ux , effective spectral 
temperature, T v (measured in MeV), and neutrinospheric radius, R u , are measured 
at infinity, and the radiated neutrino spectrum is assumed to be represented by a 
Maxwell-Boltzmann distribution, in which case (E u ) = 3 T u . The “greyness factor” 
0 is determined by numerical simulations to range between ~0.4 and ^0.8 [33,116]. 

• The luminosities of v e and v e are enhanced by an accretion component that strongly 
depends on the progenitor star through the mass-accretion rate of the infalling 
stellar core and the mass and radius of the assembling proto-neutron star (see 
Fig. 37). The sum of v e and v e luminosities can be expressed as 

( 8 ) L Ve + L= 2 + p 2 GM ™ M . 

-tins 

The first term on the r.h.s. represents the “core component” of the luminosity 
carried by neutrinos diffusing out from the high-density inner region of the proto¬ 
neutron star. This component can be assumed to be similar to the heavy-lepton 
neutrino luminosity, because the core radiates neutrinos in roughly equal numbers 
from a reservoir in thermal equilibrium, which is confirmed by the close similarity 
of the luminosities of all neutrino species found after the end of accretion. The sec¬ 
ond term on the r.h.s. of Eq. (8) stands for the accretion component of neutrinos 
lost from the hot, inflated accretion mantle around the neutrinosphere. It can be 
approximately represented by the product of mass accretion rate, M, and Newto¬ 
nian surface gravitational potential of the NS, 4> = GM ns /R ns , with neutron-star 
mass M ns and radius R ns . By a least-squares fit to a large set of ID results for the 
post-bounce accretion phase of different progenitor stars, the values /?i « 1.25 and 
/?2 ~ 0.5 were deduced in Ref. [68]. These values apply later than about 150 ms 
p.b., when the postshock accretion layer has settled into a quasi-steady state, and 
they depend only weakly on the nuclear EoS. In Ref. [116] the authors used a form 
similar to Eq. (8) with f3\ = 1; they found ~ 0.5-1 for the pre-explosion phase 
of the 2D models they investigated. 
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• During the accretion phase the mean energies of all neutrino species show an over¬ 
all trend of increase, which is typically steeper for v e and v e than for v x . Interest¬ 
ingly, for progenitors more massive than ~10 M©, which possess extended accretion 
phases, this leads to a violation of the canonical hierarchy of the mean neutrino 
energies (defined as ratios of energy to number fluxes), i.e., one observes a moment 
when the usual order (E„ e ) < (E Pe ) < (E Ux ) changes to (E Ue ) < (E Ux ) < (E Pe ) 
(see Ref. [141]). The reason for this behavior is a local temperature maximum 
somewhat inside of the neutrinospheres of v e and z/ e > which forms because of the 
compressional heating of the growing accretion layer, whose cooling becomes inef¬ 
fective when the radius of the proto-neutron star shrinks and the density profile 
becomes steeper [68,116]. The hierarchy inversion is enhanced and shifted to a con¬ 
siderably earlier time when energy transfer in neutrino-nucleon scattering is taken 
into account. Non-isoenergetic neutrino-nucleon scattering reduces the mean ener¬ 
gies of v x in a “high-energy filter” layer between the v x number sphere and the v x 
scattering sphere [31,38]. The corresponding energy transfer to the stellar medium 
also raises the luminosities and mean energies of v e and v e . Different from the mean 
energies, (E„), the rms energies, \J (E^), always obey the canonical hierarchy. 


• The secular rise of the mean energies of the radiated v e and v e during the accretion 
phase is rather well captured by the proportionality (E u ) oc M ns (t ) between the 
mean energies and the growing mass of the proto-neutron star. The proportionality 
constant depends on the neutrino type but is only slightly progenitor-dependent 
[116]- 


During the accretion phase, convection inside of the massive NS of Fig. 9 causes 
differences compared to the non-convective case similar to those already mentioned for 
the low-mass remnant of a 9.6 M© star above. It is stressed again that the transition 
from accretion to proto-neutron star cooling exhibits features that differ between the two 
models shown in Fig. 9, because in both cases the artificial initiation of the explosion 
was handled slightly differently. During the Kelvin-Helmholtz cooling the influence of 
convection leads to growing effects and prominent differences in the long-time evolution of 
the neutrino emission properties. Convection enhances the luminosities over the first few 
seconds and thus reduces the late-time luminosities (Fig. 9, right upper panel), because 
the NS cooling and deleptonization are accelerated by convective transport effects; see 
also Fig. 4 for the faster evolution of the temperature profile in the convective model. 
Correspondingly, the mean neutrino energies are increased by up to ~ 1.5 MeV in the first 
seconds, followed by a faster decline because of the cooler NS at late times (Fig. 9, right 
middle panel). 


2*4.3. Spectral properties of the neutrino emission. The spectrum of radiated neutrinos 
is usually somewhat different from a thermal spectrum. Since neutrino-matter interac¬ 
tions are strongly energy dependent, neutrinos of different energies decouple from the 
background medium at different radii with different temperatures of the stellar plasma. 
Nevertheless, the emitted neutrino spectrum can still be fitted by a Fermi-Dirac distribu¬ 
tion, f(E) oc E 2 [l + exp (E/T — 77 )] —1 , with fit temperature T and effective degeneracy 
parameter p [142]. A mathematically more convenient representation was introduced in 
Ref. [31], who proposed the following dimensionless form for the energy distribution at 
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a large distance from the radiating source: 

(9) /.(£) *(A), 

where 


/ 0 °°d EEUE) 

J 0 °°dEf a (E) 


is the average energy( 4 ). The parameter a represents the amount of spectral pinching 
and can be computed from the two lowest energy moments of the spectrum, (. E ) and 
<£ 2 >,by 


( 11 ) 


(E 2 ) _ 2 + a 

Je )2 “ 1 + a ' 


Higher energy moments (E £ ) for i > 1 are defined analogue to Eq. (10) with E £ under 
the integral in the numerator instead of E. Besides its analytic simplicity, this functional 
form has the advantage to also allow for the representation of a wider range of values 
for the spectral (anti)pinching than a Fermi-Dirac fit. A Fermi-Dirac spectrum with 
vanishing degeneracy parameter (77 = 0) corresponds to a ~ 2.3, a Maxwell-Boltzmann 
spectrum to a = 2, and a > 2.3 yields pinched spectra whereas a < 2.3 gives antipinched 
ones. Comparing to high-resolution transport results, the authors of Ref. [143] showed 
that these “a-fits” also reproduce the high-energy tails of the radiated neutrino spectra 
very well. 

In Fig. 9 the shape parameters a as functions of time are given in the bottom panels. 
The evolution is similar in the convective and non-convective models, with moderate 
quantitative differences in some phases. The modest effects of convection are not as¬ 
tonishing, because the convective layer is located well inside of the neutrinospheres and 
the neutrinospheric layer itself is convectively stable. The region where the radiated 
spectrum is shaped is therefore dominated by neutrino radiative transport and not by 
convective transport. In particular neutrinos of higher energies, around and beyond the 
spectral peak, which have the strongest influence on the higher energy moments (E £ ) 
with t > 2, are still diffusively coupled to the stellar medium outside of the convective 
layer. Therefore a major change of the shape and time evolution of the emitted spectra 
due to convection cannot be expected. 

Figure 9 shows that the v e spectrum during the shock breakout burst has a consid¬ 
erable pinching (a up to nearly 7), but the pinching decreases with time. The opposite 
trend occurs for v e during the early post-bounce evolution (f p b>3ms), while heavy- 
lepton neutrinos have spectra closest to the Maxwell-Boltzmann form. As the growing 
and contracting accretion mantle of the proto-neutron star heats up during the accretion 
phase, the pinching of the spectra of all types of neutrinos decreases. During the subse¬ 
quent cooling phase, v e and v x exhibit fairly stable shape parameters (a ~ 2.7 for v e and 
a « 2.2 for v x \ whereas the spectra of v e become more similar to those of v x as time goes 


( 4 ) With a being determined, a normalization factor A relates the spectrum to the neutrino 
luminosity: L = 4tt r 2 A f o °° d EEf a (E) = 4tt r 2 A(E) d Ef a (E). 
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Fig. 10. - Similar to Fig. 9 but for ID collapse and explosion simulations of the 27 M© star that 
include nucleon self-energy shifts in the charged-current neutrino-nucleon interactions in one 
case (model LS220-s27.0co) while ignoring these self-energy potentials in the other case (model 
LS220-s27.0c). Both models were performed with a mixing-length treatment of convection. 
Note that similar to the models shown in Fig. 9, also here the explosions of the two cases were 
triggered slightly differently, for which reason the simulations exhibit artificial differences in 
the time interval from ~0.5s to ~0.8s p.b. without any deeper physical meaning and without 
any significant perturbing consequences for the subsequent evolution. The effects of nucleon 
self-energy shifts during the z/ e -burst and accretion phases are small as pointed out before [144]. 
During the proto-neutron star cooling evolution, there are visible differences of up to about 10% 
in the neutrino luminosities. The main, persistent effect of the nucleon mean-field potentials 
is a reduction of the average energy of the emitted v e by up to ~0.7MeV. At very late times 
(tpb 6 s) the relative enhancement of the luminosity of v x for the case with nucleon self-energies 
increases and the mean energies of v e and v x become ~lMeV higher than in the other model. 


on. In the convective model, different from the non-convective one, the pinching of the 
v e and v x spectra increases again at late times. This reflects the much faster progression 
of cooling in the convective model (cf. Fig. 4), where the late evolution towards more 
pinching of the v e and v x spectra occurs earlier than in the non-convective case. 

2*4.4. Nucleon self-energy shifts in /3-reactions . Besides convection in newly formed 
NSs, the EoS of hot NS matter and the mean-field potentials of the nucleons play an 
important role for the neutrino emission and cooling evolution of the compact remnant. 
The importance of the EoS-dependent nucleon self-energies has been pointed out in 
the recent past [144-148] after the corresponding physics had been ignored in earlier 
studies of proto-neutron star cooling [33,34,149] and of proto-neutron star accretion to 
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BH formation [67]. Considering neutrons and protons as gases of nonrelativistic quasi¬ 
particles that move in a single-particle mean-field potential, E7, one can write their energy- 
momentum relation E^(pi) in generalization of the non-interacting case as 

2 

(12) Ei(pi) = + WjC 2 + Ui for i = n,p, 


with rrii and m* being the particle rest masses and the (Landau) effective masses, respec¬ 
tively. Both the masses m* and the potentials Ui depend on density, temperature, and the 
neutron-to-proton ratio (or the electron-to-nucleon ratio Y e ). The electron and positron 
energies of the reactions of Eqs. (1) and (2) are therefore related with the energies of 
their corresponding v e and u e by 


(13) 

(14) 


E v 


E P - + 


2m* 


E Pp — E e + + 


2m* 


- (m n - m p )c 2 - (U n - Up), 

2^r + ( m n - m p)c 2 + ( U n - Up) . 


The potential difference A U = U n — U p is directly related to the nuclear symmetry energy 

[148]: 

(15) A U - 4(1 - 2 Y e )E^(p ), 

where E™^ is the interaction part of the symmetry energy at zero temperature, and the 
r.h.s. ignores terms of higher than linear order in Y e . Since A U is positive at the very 
neutron-rich conditions in forming NSs, the contribution from the mean-field potentials 
reduces the energy of v e and increases the energy of v e created in the /^-reactions. More¬ 
over, the mean-free path for v e absorption on neutrons is decreased because the extra 
energy associated with the potential difference A U diminishes the importance of Pauli 
phase-space blocking in the final state of the electron, whereas the mean-free path for v e 
is increased at low energies and decreased at high energies. 

As a consequence of these changes, the authors of Ref. [144] reported a strong impact 
of the nucleon self-energy shifts on the luminosities and spectra of all neutrinos, but 
with an important difference between pre-explosion and post-explosion phases. While 
they found the nucleon mean-field potentials in the charged-current reactions to have a 
negligible influence during the accretion phase prior to explosion, they obtained reduced 
luminosities for all neutrino species and a considerable decrease of the mean energy of 
the emitted v e during the Kelvin-Helmholtz cooling of the proto-neutron star. They also 
observed values for the electron fraction in the neutrino-driven wind (see Sect. 2*4.5) that 
were significantly lower (down to F e ,wind ~ 0.45) than in the simulation without nucleon 
mean-field potentials. Using a more accurate determination of Y^wmd, Ref. [29] obtained 
less neutron-rich ejecta with T e ,wind ~ 0.48 despite having replaced the TM1 EoS [150] 
of the previous work by a DD2 EoS table provided by M. Hempel [147], which predicts 
higher values of AU between the neutrinospheric layer and the nuclear saturation density, 
from where most of the neutrino emission originates. 

Using the STOS TM1 EoS [150] and the LS220 EoS [90], Ref. [68] confirmed the 
basic trends of the results of Ref. [144]. Similar absolute shifts were obtained for both 
investigated EoSs, although the symmetry energy of the STOS EoS is larger. Ref. [68] 
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Fig. 11. — Evolution of the proto-neutron star radius (upper left) and of characteristic parameters 
of the neutrino-driven wind: expansion time scale from 5 x 10 9 K to 2 x 10 9 K (upper right), 
asymptotic electron fraction, P e ,wmd (bottom left), and asymptotic entropy (bottom right)', the 
last two quantities were evaluated at a radius of 1000 km. The two ID simulations for the 
9.6 Mq progenitor (LS220-z9.6co and SFHo-z9.6co) were performed with different hadronic EoSs 
and produced explosions self-consistently, whereas the explosion of the 27.0 Mq star in ID 
was triggered artificially at 0.5 s after bounce. All simulations were performed with a mixing- 
length treatment of proto-neutron star convection and including nucleon self-energy shifts in the 
charged-current neutrino-nucleon reactions. All wind ejecta are proton-rich. 


also found that nucleon self-energy potentials lead to reduced luminosities for all neutrino 
species with the biggest effects for v e . While the average energy of the emitted v e is 
lowered by roughly 0.75 MeV compared to the case without self-energy shifts, the effect 
is about one third of that for v x and even smaller for u e , in contrast to the slight rise 
seen in Ref. [144] ( 5 ). More important, however, is the fact that mixing-length convection 
shrinks the differences between the results with and without nucleon self-energies to 
about half of the values obtained without convection. Because convection speeds up 
the transport of electron-lepton number out of the NS, it enhances the v e emission and 
partly compensates the effects of the mean-field potentials in the neutrino opacities, at 


( 5 ) Our simulations differ from those of Ref. [144] in many aspects of the neutrino opacities. 
In particular, our models take into account weak magnetism and recoil effects in neutrino- 
nucleon interactions, while these relevant corrections were ignored in the calculations discussed 
in Ref. [144] (Meng-Ru Wu; private communication). 
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Fig. 12. - Similar to Fig. 9 but for two ID collapse and explosion simulations of a 9.6 Mq star 
with two different nuclear EoSs. Model LS220-z9.6co uses the EoS of Ref. [90] with a nuclear 
incompressibility modulus of K — 220 MeV, whereas model SFHo-z9.6co employs the SFHo 
hadronic SN EoS of Ref. [151]. The cooling evolution of the proto-neutron star is clearly different 
for both EoSs. Both models were performed with a mixing-length treatment of convection and 
in both cases the explosions developed self-consistently within 170-200 ms after bounce. Note 
that different from Figs. 8-10, the neutrino emission properties are shown here for z/ e , Ve, one 
species of or z/ r , and one species of or z> T , because the neutral-current scattering cross 
sections of neutrinos and antineutrinos differ in the sign of the weak-magnetism corrections [152]. 
The slightly lower nucleon-scattering opacity of leads to slightly higher luminosities and 
up to ~1 MeV higher mean energies. In all other simulations discussed in this Section, the four 
heavy-lepton neutrinos were treated equally by one representative species v x , whose emission 
properties turned out to agree well with the average of those of and z/^r [65]. 


least as long as the deleptonization of the proto-neutron star is proceeding. Figure 10 
displays two simulations for our 27 M 0 progenitor for comparison. As pointed out before 
in Ref. [144], the nucleon self-energy shifts are not relevant during the shock-breakout and 
accretion phases. However, even during the proto-neutron star cooling the differences are 
small until very late times, t p b ^ 6 s. The main consequences are a persistent reduction of 
the mean energy of the radiated v e by at most ^0.7 MeV and no late-time convergence of 
all mean energies as observed in the case without nucleon self-energies. Also the spectral 
pinching of the v e is moderately lowered when the mean-field potentials are taken into 
account. 

2‘4.5. Neutrino-driven outflows from proto-neutron stars. Neutrinos leaving the hot 
proto-neutron star deposit energy just outside of the neutrinosphere and thus launch a 
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Fig. 13. - Mass-shell plots for the evolution of a collapsing progenitor star with a birth mass 
of 40 Mq (model s40s7b2 of [153]) towards BH formation without mixing-length treatment 
of proto-neutron star convection (upper panel ) and including such a treatment (lower panel). 
The simulations made use of the LS220 model of the nuclear EOS ([90] with incompressibility 
modulus of K — 220 MeV). The color coding shows the Brunt-Vaisala frequency as defined in 
Eq. (5). The different lines are explained in the legend of the plot. The thick blue, solid line is 
the stalled and retreating shock, the ingoing solid lines indicate continuous accretion of infalling 
matter, the thin dashed lines the contraction of the proto-neutron star. Convection slightly 
stretches the evolution time until the proto-neutron star becomes unstable to BH formation, 
because the accelerated deleptonization leads to a faster heating of the proto-neutron star interior 
and therefore enhanced thermal support. 


baryonic mass outflow, the so-called “neutrino-driven wind” (e.g., [155,156] and refer¬ 
ences therein). This wind is discussed as a potentially important site of heavy-element 
formation (for recent reviews, see Refs. [157,158]). The neutron-to-proton ratio in the 
ejected matter is set by v e and u e absorption via the processes of Eq. (2). If neutrons 
dominate over protons, the wind environment may provide favorable conditions for the 
formation of neutron-rich nuclei in the rapid-neutron capture process (r-process), de¬ 
pending on sufficiently high entropies and/or sufficiently short expansion time scales of 
the ejecta (cf., e.g., [159-161]). If the outflow material develops a proton excess, nuclei 
on the proton-rich side of the valley of stability might form through the neutrino-proton 
process [162,163]. 

The neutrino-driven winds computed in our models are proton-rich during the whole 
proto-neutron star cooling evolution. Figure 11 displays two simulations for the 9.6 Mq 
progenitor and one for the 27 M 0 star, which bracket the conditions that can be expected 
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Fig. 14. - Same as Fig. 13, but for the 40 M© progenitor s40.0 of [154], which has a lower 
mass-accretion rate and therefore needs a much longer evolution until BH formation. 


for proto-neutron stars in the range of baryonic masses of 1.36-1.8 M© (final gravitational 
masses of 1.25-1.6 M©), which are representative for the far majority of core-collapse SNe. 
All three displayed simulations include the mixing-length treatment of convection and the 
mean-field potentials of nucleons in the /3-processes. The lower-mass case was computed 
with two different nuclear EoSs, both yielding qualitatively similar time evolutions of 
the wind parameters. Instead of an early wind phase of about 3 s with T ejW ind < 0-6, 
which was reported as a consequence of the nucleon self-energy shifts in Ref. [144], the 
simulations shown in Fig. 11 exhibit a prominent maximum of Y^ jW i n d >0.5 during the 
early cooling phase. This maximum is a consequence of proto-neutron star convection, 
which was not included in the models of Ref. [144]. Correspondingly, the maximum 
is also absent in simulations that we performed without the mixing-length treatment of 
convection or for conditions (e.g., nuclear EoSs) that yield only weak convection, in which 
cases T e?w i n d increases monotonically during the proto-neutron cooling [68]. Convective 
enhancement of the lepton-number transport inside of the nascent NS increases the v e 
number flux relative to that of v e and, in particular, leads to increased pinching of the v e 
spectrum (see Fig. 9), thus reducing v e captures on protons compared to v e absorptions 
on neutrons in the outflowing wind material. This favors higher values of T e?w i n d until the 
main period of deleptonization is over and the proto-neutron star convection becomes 
weaker at around 3-4 s after bounce (cf. Fig. 4), at which time T e?w i nc j declines from 
a local maximum. Only subsequently, after having gone through a local minimum on 
the proton-rich side, T ejW i n d resumes the trend found in Ref. [144] with increasing values 
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towards later times. The entropies in the neutrino-driven outflows rise continuously until 
they reach maximum values of 120-140 &b per nucleon with the higher values for the more 
massive NS. Such conditions of entropy and proton excess clearly disfavor the formation 
of r-process elements in this environment, confirming previous results of less complete 
models [33,34]. 

Active-active as well as active-sterile neutrino-flavor oscillations modify the neutrino 
emission from the proto-neutron star and can thus have an impact on the neutron-to- 
proton ratio in the neutrino-driven wind. This requires that the flavor-changing effects 
occur interior to the neutrinosphere or in the close vicinity of the nascent NS, i.e., in 
the region where neutrino interactions set the electron fraction of the outflowing matter. 
While quite a number of studies have explored the implications of matter-enhanced neu¬ 
trino conversions to active (e.g., [164]) and sterile channels (e.g., [165-171]) in relevant 
parts of the space of possibilities for mixing parameters and scenarios, the implications 
of additional neutrino-background induced flavor conversions are not yet satisfactorily 
understood. First investigations have not been able to reveal promising conditions for 
r-processing in the wind ejecta [170,171], but more studies that account for the growing 
insight into the complex physics of collective flavor transformations are urgently needed 
(see also Sec. 4). 

2*4.6. Equation-of-state dependence of the neutrino emission. The neutrino signal from 
the proto-neutron star cooling phase of a future Galactic SN explosion has the potential 
to provide us with extremely valuable information about the nuclear EoS that describes 
the physics in the NS interior. Figure 12 shows interesting EoS dependent effects of 
the neutrino signal radiated by the NS forming in the collapsing and exploding 9.6 Mq 
progenitor. Both displayed simulations were carried out with different models of the 
high-density EoS above 10 11 gcm -3 and both include a mixing-length treatment of NS 
convection as well as the nucleon self-energy effects in the /^-processes. During the shock- 
breakout and the post-bounce accretion phase, which in this low-mass model lasts only 
until ^200 ms after bounce, the differences of the neutrino luminosities and mean energies 
are still modest, with slightly higher values for the LS 220 EoS. At later times, however, 
the differences become considerable and the simulation with the LS 220 EoS yields 10 - 
20% higher neutrino luminosities and mean energies until ^3 s after bounce. This can be 
understood as a consequence of the more rapidly contracting proto-neutron star and the 
smaller NS radius over a period of ^6 s in the case of the LS 220 EoS (see left upper panel 
of Fig. 11 ). Also differences in the efficiency of convective transport in the NS interior 
play a role, because the region of convective instability is sensitive to the properties of the 
nuclear EoS, in particular to the symmetry energy of the matter composed of neutrons 
and protons [172]. At intermediate times during the cooling evolution, 3 s < £ p b ^ 10 s, the 
order of the luminosities and mean energies for the two investigated EoSs is reversed. Now 
the NS computed with the SFHo EoS exhibits considerably higher values of the neutrino- 
emission parameters, because the remnant with the LS 220 EoS has cooled faster and its 
temperatures have become correspondingly lower. As a consequence of the accelerated 
cooling of the latter model, the time £^90 at which 90% of the gravitational binding 
energy of the cold, compact remnant have been radiated away is considerably shorter for 
LS220. The effect is so extreme that for our ^1.36 M 0 (baryonic mass) NS computed 
with the SFHo EoS as well as for our ~ 1 . 8 M 0 (baryonic mass) NS described by the 
LS220 EoS, £^90 ~ 6 s is roughly the same, although the total release of gravitational 
binding energy is 198 B = 1.98 x 10 53 erg in the former case and 335 B = 3.35 x 10 53 erg 
in the latter case. 
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Fig. 15. - Neutrino signals (luminosities and mean energies for an observer in the laboratory 
frame) for v e (top), v e (middle), and one species of v x (bottom) for the 40 M® BH formation 
simulations shown in Fig. 13 (left panel) and Fig. 14 (right panel). The results of the models 
including mixing-length convection in the proto-neutron star (black solid lines) exhibit only 
relatively small differences compared to the models without convection (red dashed lines). In 
the case of the s40s7b2 progenitor (left) convection mainly reduces the luminosities and mean 
energies of heavy-lepton neutrinos, whereas in the s40.0 case (right) it reduces the mean energies 
of v e and v e by up to ~1 MeV. 


Such differences indicate the paramount importance of the high-density EoS model 
for reliable calculations of the neutrino-cooling signal of nascent NSs. They also suggest 
that highly valuable constraints of the still incompletely understood properties of matter 
in the interior of NSs could be deduced from the measurement of the neutrino signal 
from a future Galactic SN. If the distance to the SN is well determined, the energy-loss 
time scale and calorimetric arguments can yield useful information. But even if the SN 
distance is not known, the characteristic time structure of the signal [172] and the spectral 
evolution can provide crucial evidence. This will require most accurate predictions of 
the SN neutrino emission (employing consistent microphysics of the EoS and neutrino 
opacities) and the combined information from detections in all existing, upcoming, and 
planned experimental facilities. All together this will allow one to interpret the detailed 
features of a high-statistics measurement of the time structure and spectral evolution 
of the neutrino emission during all phases of stellar collapse and explosion, despite the 
complexities of the SN physics and of various neutrino oscillation phenomena on the 
way from the SN core to the detectors on earth. Because of the stochasticity of multi¬ 
dimensional hydrodynamics phenomena that assist shock revival and because of the 
limited sensitivity of the early-time neutrino emission to the high-density EoS, the post¬ 
explosion emission, in particular the late-time neutrino signal, will be the most reliable 
probe of the state of matter at nuclear and supernuclear densities in the deep interior of 
the NS. 

2‘5. Neutrino signals from black hole formation. - A considerable fraction of all stellar 
core collapses, possibly similarly frequent as SNe, might fail to explode. In these cases 
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the newly formed NS continues to accrete matter from the infalling stellar core until it 
reaches the mass limit for BH formation and gravitational instability sets in. Modern 
calculations of the neutrino emission of such events, performed with energy-dependent 
neutrino transport schemes, were carried out, e.g., in Refs. [67,173-177] for different 
EoSs of hot NS matter. However, these models did not include convection in the NS and 
also ignored many aspects of the neutrino interactions, e.g., did not take into account 
neutrino-pair conversion between different flavors [178], energy transfers in neutrino- 
nucleon interactions, and nucleon-correlations in the dense medium. 

Here we present models of the Garching group [68] that include all of these effects. Fig¬ 
ures 13 and 14 show mass-shell plots of two 40 M© progenitors (s40s7b2 of Ref. [153] and 
s40.0 of Ref. [154], respectively), both of which were computed with the LS220 EoS. The 
upper panels of the figures display results of simulations without mixing-length treatment 
of convection, the lower ones with convection inside of the proto-neutron stars. Figure 15 
provides the time evolution of the corresponding neutrino luminosities and mean ener¬ 
gies. The differences between convective and non-convective models are relatively small 
because the convective shell is buried below a hot, convectively stable accretion layer 
that is growing in mass and whose neutrino production dominates the properties of the 
radiated neutrinos. The models with active proto-neutron star convection exhibit slightly 
lower values of some quantities of the neutrino emission. Most remarkably, the convective 
models tend to collapse to BHs a bit later than their counterparts without mixing-length 
convection, despite the higher values of the integrated energy and lepton-number loss 
through neutrinos. This delay is explained by the stabilizing influence of thermal pres¬ 
sure, because the accelerated deleptonization associated with convective lepton-number 
transport leads to a faster rise and higher values of the temperature in the deep interior 
of the NS. 

The two considered 40 M© stars differ significantly with respect to the duration of 
their evolution until BH formation. While it takes ^0.57 s after bounce for model s40s7b2 
to reach this moment, the gravitational instability in model s40.0 occurs as late as ^2.1 s 
after bounce, because the mass accretion rate of the collapsing stellar core in this latter 
case is much lower after a composition-shell interface has passed the stalled shock at 
^0.3 s after bounce. Correspondingly, the neutrino luminosities in the s40.0 run are 
only roughly half of those in the s40s7b2 model; also the mean energies of the radiated 
neutrinos rise more slowly and remain lower in model s40.0 when the same times are 
compared. For heavy-lepton neutrinos these differences in the time-dependent behavior 
are most prominent. While their luminosity and mean energy grow steeply in the s40s7b2 
simulation, L Ux in s40.0 is basically constant for t p b ^ 0.3s and (E Ux ) increases only with 
a flatter gradient (Fig. 15). These differences reflect the compression and corresponding 
heating of the rapidly mass-accumulating accretion layer in model s40s7b2, whereas in 
s40.0 neutrino cooling can well keep pace with the inflow of energy by the accretion of 
fresh gas, for which reason the mantle layers of the newly formed NS in the latter model 
remain cooler. 

Failed SNe, i.e., “direct” BH formation events in which the transiently formed NS 
is pushed across the BH formation limit by permanent post-bounce accretion, therefore 
produce neutrino signals with distinct features like a continuous hardening of the emitted 
spectra before the emission terminates abruptly. The two discussed simulations show that 
the detailed properties of the neutrino radiation depend strongly on the survival time 
of the NS. Better predictions of the integrated neutrino emission of such events, whose 
contribution might imprint characteristic structures on the DSNB, will therefore require 
systematic calculations of failed explosions for large model sets of progenitor stars, taking 


35 



Supernova Neutrinos: Production, Oscillations and Detection 


into account the growing understanding of stellar core-collapse cases that are not likely 
to blow up as SNe [57,58,179-181]. 

One should keep in mind, however, that spherically symmetric simulations (even if 
they include a mixing-length treatment of convection interior to the neutrinosphere) are 
unlikely to capture all effects of relevance for the neutrino emission. When the stalled 
shock retreats in reaction to the contraction of the NS, favorable growth conditions of 
the SASI are established and, correspondingly, episodes of violent SASI sloshing and 
spiral-mode activity are found in 3D simulations that follow the post-bounce accretion of 
non-exploding massive stars over hundreds of milliseconds [70,91,92,130,131]. It must be 
expected that large-amplitude neutrino-emission variations caused by SASI modulations 
of the NS accretion (cf. Sect. 2*3) accompany the whole evolution at least until the 
collapse of the NS to a BH sets in. Future 3D simulations will also have to clarify the long¬ 
time development of the LESA dipole that was found for the neutrinos radiated during 
the accretion phase [102]. It is presently unclear at which level these and other multi¬ 
dimensional phenomena (e.g. connected to stellar core rotation, whose effects potentially 
increase on the way to NS collapse) can affect the predictions of the DSNB from failed 
SNe (see Sect. 5). Although most (maybe all) of the direction dependences and short- 
time signal modulations of the observable neutrino emission will average away in the time 
and direction-integrated signals that are relevant for the DSNB, the interpretation of a 
Galactic BH formation event will greatly benefit from a precise knowledge of the signal 
structure until the emission is terminated. 

2*6. Conclusions. - In this Section we discussed the current status of modeling the 
neutrino emission from stellar core collapse events. Hydrodynamical simulations with 
most sophisticated neutrino transport are needed to calculate detailed signal properties, 
which are necessary for studies of neutrino oscillations in super novae, for investigations 
of neutrino-dependent and neutrino-induced nucleosynthesis, and for the exploration of 
detection possibilities of the DSNB and of neutrinos from a future Galactic supernova. 

The neutrino emission from core-collapse supernovae is strongly affected by multi¬ 
dimensional hydrodynamic instabilities inside and around the newly formed NS. Con¬ 
vective overturn and the SASI (standing-accretion shock instability) in the postshock 
accretion layer support the initiation of the explosion and impose large-amplitude time 
variations and large-scale directional asymmetries on the neutrino emission during the 
accretion phase (Sec. 2‘3; [91,92]). Proto-neutron star convection is active interior to the 
neutrinospheres during this phase as well as during the post-explosion Kelvin-Helmholtz 
cooling phase of the compact remnant. 

Recent 3D simulations with detailed neutrino transport have led to the discovery of 
the growth of hemispheric differences in the strength of proto-neutron star convection on 
a time scale of about 150 ms after core bounce (Sec. 2*3). These hemispheric differences 
produce a large dipolar asymmetry of the lepton-number flux with relative enhancement 
of the v e emission compared to v e in one hemisphere and increased v e emission on the 
opposite side. The hemispheric differences of the lepton-number (y e minus D e ) flux nor¬ 
malized to the average can be factors of a few, corresponding to dipole amplitudes of the 
individual emission of v e and v e of 10-20%, whereas the v x emission dipole is only of or¬ 
der percent [102]. While this unexpected and stunning new phenomenon is not yet fully 
understood and still needs further exploration and confirmation by independent results, 
the first successful explosions in 3D simulations with sophisticated, energy-dependent 
neutrino transport have become available [7,83,124,126]. They provide support for the 
viability of the neutrino-driven mechanism for driving and powering supernova explo- 
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sions. 

During the Kelvin-Helmholtz cooling phase convection interior to the neutrinospheres 
accelerates the lepton-number and energy loss of the nascent NS. The basic effects of this 
quasi-stationary convective transport can be well described in spherically symmetric sim¬ 
ulations by a mixing-length treatment, whose results we showed to compare nicely with 
those of direct hydrodynamic modeling of convection in 2D calculations (Sec. 2‘4). We 
also presented here the first time-dependent, hydro dynamical simulations with state-of- 
the-art spectral neutrino transport, in which a mixing-length treatment was applied dur¬ 
ing the whole cooling evolution of proto-neutron stars. The calculations were performed 
with different nuclear equations of state and include the recently intensely discussed 
effects of nucleon self-energy shifts in the charged-current /3-reactions of free nucleons. 

Convection interior to the neutrinospheres facilitates the faster deleptonization of the 
high-density core of the NS and thus enhances the emission of v e relative to v e dur¬ 
ing the first seconds, thus leading to more proton-rich conditions in the neutrino-driven 
wind. The possibility of a z/-p process [162,163] during the early phase of proto-neutron 
star cooling has to be reinvestigated. However, with the neutrino physics used (which 
excludes effects of active or sterile neutrino-flavor oscillations in the hydrodynamic mod¬ 
eling), neutron-rich conditions in the wind are not found at any time during the cool¬ 
ing evolution of compact remnants that bracket the mass range of NSs formed in the 
far majority of all supernova explosions; models with NS baryonic (final gravitational) 
masses of ^1.363 (1.252) M© and 1.776 (1.592) M 0 were discussed. Spherically symmet¬ 
ric simulations of proto-neutron star formation, however, require artificial initiation of 
the supernova explosion. Even in stars near the low-mass end of supernova progenitors, 
where multi-dimensional hydrodynamic effects are not crucial for the onset of the explo¬ 
sion [81,83], convective overturn stirs the postshock accretion layer. One-dimensional 
models therefore cannot correctly describe the transition from the accretion phase to the 
explosion and proto-neutron star cooling phase, because simultaneous inflows and out¬ 
flows coexist in the postshock layer during this transition. Two-dimensional explosion 
models show the ejection of matter with moderate neutron excess in rapidly expanding 
plumes of buoyant plasma during the early expansion phase of the shock front in par¬ 
ticular in low-mass oxygen-neon-magnesium and iron-core progenitors [84,117]. These 
ejecta might provide the conditions for a weak r-process with production of neutron-rich 
nuclei up to mass numbers of A ~ 100-110 [84]. 

Continuing accretion of matter onto the nascent NS in the case of unsuccessful ex¬ 
plosions unavoidably leads to BH formation. Recent insights suggest that these failed 
supernovae might be more common than previously thought (e.g., [56-59,179-181]) and 
might therefore account for a significant contribution to the DSNB. Neutrino signals from 
such events up to the moment when the NS becomes gravitationally unstable were com¬ 
puted with energy-dependent neutrino transport before [67,173-177]. Here we presented 
state-of-the-art ID hydrodynamic simulations that include proto-neutron star convection 
by a mixing-length treatment and the full set of currently discussed neutrino interactions 
of relevance in the core of collapsing stars (Sec. 2*5). While convection inside of the neu¬ 
trinospheres does not have a major impact on the calculated neutrino emission (except 
for slightly stretching the life time of the transiently existing NS), it must be expected 
that violent SASI sloshing and spiral-mode activity induces large-amplitude modulations 
on the neutrino luminosities and spectra radiated by the accreting NS all the way until 
BH formation. Three-dimensional simulations of the evolution from the onset of stellar 
core collapse to the gravitational instability of the NS will therefore have to be per¬ 
formed also for obtaining reliable predictions of the neutrino (and gravitational-wave) 
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signals from failed supernovae in order to get prepared for the —not unlikely— case of 
a future Galactic event. 

Neutrino data and structural profiles from the Garching simulations can be made 
available upon request for download at www.mpa-garching.mpg.de/ccsnarchive. 
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3. — Perspectives for Future Supernova Neutrino Detection 

Authors: K. Scholberg, A. Mirizzi, I. Tamborra 

In this Section, we will consider the physics of neutrino detection in the regime relevant 
for SN neutrinos, and discuss the interaction channels for relevant detector types. We 
will survey main existing and planned detectors sensitive to SN neutrinos, and touch on 
some topics related to real-time SN neutrino astronomy and astrophysics. 

3*1. Supernova neutrino detection. - Few detectors are constructed with the primary 
aim of SN neutrino detection; most neutrino detectors are built for studies of neutrino 
oscillation with solar, atmospheric, reactor and beam neutrinos, or for high-energy as- 
trophysical neutrino observation, or for nucleon decay searches. Happily, however, many 
large neutrino detectors — especially those sited underground — have excellent capabil¬ 
ities for the capture of a core-collapse SN neutrino burst. 

Detectors suitable for SN neutrinos must in general be sensitive to the products of 
interactions of neutrinos in the few to few tens of MeV range. Typically the final-state 
particle energies are also in that energy range, although sometimes can be lower. Neutri¬ 
nos interact with detector materials via both charged-current (CC) and neutral-current 
(NC) channels. For standard-model CC processes, the resulting lepton type depends on 
the incoming neutrino flavor, and the lepton charge is negative for a neutrino and posi¬ 
tive for an antineutrino. Supernova neutrino spectra are, however, almost entirely below 
threshold for CC production of a muon or a tau particle (the thresholds are 110 MeV 
and 3.5 GeV respectively for quasi-elastic production). Therefore CC interactions are 
accessible only for v e and v e from core collapse. The muon and tau components of the 
flux are accessible only via NC interactions. 

The neutrino interaction products depend also on nature of the target; the observabil¬ 
ity of the interaction products depends on the nature of the detection technology. The 
following subsections describe the different channels, and relevant detector technologies, 
for the different components of the SN flux. Detection of SN neutrinos is reviewed 
in [182]. 

3*1.1. Neutrino interactions with matter. The interaction products of SN neutrinos 
are detected using conventional means, typically via the energy loss of the primary par¬ 
ticle out of the interaction, or via energy loss of secondary particles. Detectors have 
different capabilities for energy resolution, angular resolution, time resolution, and parti¬ 
cle identification. They are limited usually by some combination of the intrinsic physics 
of the particle energy loss, the detection mechanism, and the parameters of the detector 
technology (i.e., by technical sophistication and cost). 

Real detectors are plagued by backgrounds. For surface or near-surface detectors, 
cosmogenic backgrounds are a serious issue, typically inducing event rates at or above SN 
neutrino energies at rates far higher than the rate expected from a burst from all but the 
nearest SNe. Cosmogenics include direct cosmic rays (muons being the most penetrating) 
and secondaries, some of which can be long-lived and produce radioactive decays in the 
SN event energy range (e.g., [183-185]). Nevertheless, some surface detectors can expect 
to collect useful SN burst events when triggered by some external observation. With even 
modest overburden, cosmogenic rates can be reduced significantly, to the point where 
the expected background counts within a typical few-tens-of-second burst are negligible. 
Cosmogenic background rates depend on overburden, detector size and geometry, and 
penetrating particle veto capability. Intrinsic detector background from radioactivity 
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from detector components and nearby materials can also be a problem, especially at 
low energy. In the regime below about 5-10 MeV of energy deposition, radioactivity 
can be dominant. Radioactive background rates depend sharply on the composition 
of the nearby environment and the care taken with detector radio-cleanliness. Finally, 
instrumental backgrounds (electronic noise, calibration sources, etc.) can sometimes be 
a problem in practice. In general, underground detectors designed to measure low-energy 
neutrinos from steady-state sources such as the Sun, the atmosphere, or reactors, are able 
to achieve background levels low enough that a Galactic SN burst will have negligible 
background. Backgrounds become more of an issue for detection beyond the Milky Way, 
and for the DSNB (see Sec. 5). 

Different detector materials have different flavor sensitivities in part by virtue of the 
relative interaction rates of CC and NC channels in the material. The ability of a detector 
to determine flavor content by distinguishing different interactions on an event-by-event 
basis (“tagging capability”) also affects the information content that will be possible to 
extract from a burst observation. We summarize here the main detection channels for 
different flavors. 

• Electron antineutrinos: The detectors that were running during the SN 1987A 

neutrino burst, as well as most current detectors, are sensitive primarily to the v e 
component of the neutrino flux. The reason is that the main detector materials for 
large underground detectors, water and liquid scintillator (C n H 2 n ) are rich in free 
protons, which have a rather large, and well known, cross section for interaction 
with v e via inverse beta decay (IBD) [186], v e —>> n + e + . This interaction has 

a threshold of E Uthr = 1.8 MeV. For typical expected SN neutrino spectra, IBD 
dominates in water and scintillator. The IBD interaction products are detected 
via the energy loss of the positron. In some cases, the neutron can be detected 
after moderation and capture (with well-defined delay on the tens or few hundred 
microsecond scale, depending on materials available for capture). For example, a 
neutron capture on a free proton, n + p —>> d + 7 , produces a 2.2-MeV gamma- 
ray, which can be detected via its Compton-scattering energy loss. Other nuclei 
(notably Gd) can also capture neutrons with high cross section, and the subsequent 
deexcitation gammas can be detected. If the time-delayed neutron can be detected 
via capture, it often provides an experiment ally-useful tag of the IBD interaction. It 
is also possible to observe the energy loss of Compton-scattered 0.511-MeV gamma- 
rays resulting from annihilation of the positron. 

Charged-current v e interactions on protons bound in nuclei of common detector 
materials (carbon, oxygen, argon, iron, lead) can also occur (y e + (TV, Z) (N — 
1, Z + l) + e - ), and may result in ejected nucleons and nuclear deexcitation gamma- 
rays in addition to the positron (e.g., [187,188]). However CC v e interactions on 
nuclei are usually subdominant, as the production of a neutron in the final state is 
often Pauli-suppressed. 

• Electron neutrinos: Electron neutrinos interact via CC channels on neutrons in 
nuclei (alas, no free-neutron neutrino target exists), according to v e + (N,Z) —)> 
(.N + 1, Z — 1) + e~k* This interaction will occur for oxygen in water, and on 12 C in 
liquid scintillator, and will represent an observable component of a large-statistics 
sample; however IBD will still be dominant for these detector materials, and it may 
be difficult to do an event-by-event tagging. As for CC v e interactions, prompt and 
delayed final-state particles besides the electron may be produced, depending on 
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the nature of the target nucleus. One of the most promising targets for clean v e 
detection is liquid argon. This has a reasonably high-cross-section CC v e channel 
v e + 40 Ar —>> e~ + 40 K* [189,190]. In a time projection chamber, the resulting elec¬ 
tron track can be reconstructed. Furthermore, the resulting characteristic cascade 
of gamma rays, for which Compton scatters can be detected, can serve to tag the 
channel as v e . This channel should dominate in argon for the expected core-collapse 
spectrum. 

Lead and iron also have a large cross sections (especially lead) for CC v e interaction 
[191-196], for which observable products are (in principle) the lepton, but also single 
or double emitted neutrons (for which there are different thresholds). 

• Muon and tau neutrinos: As noted above, in order to access the and v T (and 
antineutrino) components of the signal, NC sensitivity is required. Because NC 
reactions are flavor-blind, all flavors will be entangled in the observed NC signal; 
nevertheless, since under many scenarios the mu/tau component is hotter (and 
furthermore represents two-thirds of the luminosity), the NC signal will have a 
significant v T component. NC interaction channels occur in every detector type, 
although detectors differ in their ability to tag the channel. The NC interaction 
results in a scattered nuclear target, and is either elastic or inelastic. Inelastic 
NC scattering on nuclei produces excited final-state nuclei which can deexcite via 
gamma emission, or ejection of protons or neutrons. These deexcitation products 
are the experimental observables. While neutrino differential spectral information 
is lost, some integral spectral information may be available by counting observed 
events over threshold. In water, resulting protons and neutrons will always be 
below Cherenkov threshold; gammas may however produce Compton scatters that 
are observable. Neutrons may also be captured on protons or dopants such as Gd 
to produce above-Cherenkov-threshold light. In scintillator, all ejected products 
are in principle observable, and neutrons are generally visible via the gamma-rays 
emitted following their captures. A particularly attractive 15-MeV NC-induced 
deexcitation gamma (the cross section for which has been measured [197]) occurs 
for 12 C. 

In argon, some NC deexcitation channels may be observable [198]; however this is 
currently rather unexplored in the literature. Lead also has prominent neutron- 
ejection NC channels which can be exploited by detectors. (Note there is also a NC 
component to elastic scattering on electrons, which is discussed separately below.) 

Elastic NC scattering on protons or nuclei is also a possibility for detectors with low 
energy thresholds. Elastic scattering produces simple recoil spectra, the shapes of 
which retain source spectral information [199]. Elastic scattering on protons [200], 
which produces hundreds-of-keV recoils, is observable in scintillator, and coherent 
elastic neutrino-nucleus scattering (CEz/NS) [201-203], will also yield events with 
recoils in the tens to hundreds of keV scale. This regime is out of reach for very large 
detectors, but is attainable by detectors designed to search for weakly-interacting 
dark matter. 

• Elastic scattering on electrons: A final category of interaction, the simplest phys¬ 
ically, is sensitive to all flavors of SN burst neutrinos; this is elastic scattering on 
electrons (ES). It proceeds via both CC and NC channels [204]. The detectable 
signature is the energy loss of scattered electron. The ES cross section is low com¬ 
pared to cross sections on nuclei, typically representing a few percent of a SN burst 
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signal. This interaction has directional sensitivity— the electron is scattered with 
respect to the incoming neutrino direction within an (energy-dependent) angle of 
about 30 degrees. Because all materials have electrons, this interaction will occur 
in all detector types, although detectors vary in the efficiency for observing the 
electron and in ability to track the scattered electron direction. 

We note that although the cross sections for the simplest targets — elastic scattering 
on electrons and inverse beta decay of v e on protons— are well understood, and have been 
employed as low-energy detection interactions for decades with e.g., solar and reactor 
neutrinos, the interactions of neutrinos on heavier nuclei are quite poorly understood. 
At present, the only other interactions studied experimentally with better than ^ 10% 
precision is 12 C [197,205]. The targets of interest for SN neutrino detection have never 
been irradiated with a well-understood neutrino source to determine their cross sections 
and the properties of the interaction products (particle types, angular distributions, etc.). 
Furthermore, there are relatively few theoretical studies in the relevant energy regime. 
We note promising opportunities to study these cross sections using a pion decay-at-rest 
source, such as the Spallation Neutron Source at Oak Ridge National Laboratory in 
Tennessee, within the next decade [206]. 

3‘1.2. Main detector types and considerations. There are several kinds of relevant 
detector technologies, some currently running and others planned for the future. The 
most important and the main instances of these are described below. 

• Water Cherenkov (WC): These employ water in liquid or ice form as the detec¬ 
tor material. Optical-frequency Cherenkov radiation from particles moving faster 
than the speed of light in water is collected by photomultiplier tubes. The Super- 
Kamiokande detector, a 50-kton detector in Japan [22] is the main instance of 
this class. Among detectors running at the time of this review, Super-Kamiokande 
will collect the largest number of individually-reconstructed SN neutrino events. 
A possible enhancement to Super-Kamiokande currently under investigation is the 
potential loading with Gd to improve the neutron capture cross section for better 
tagging of v e [207] (and v e [208]). The tagging efficiency for v e via neutron capture 
on free protons is about 18%. A Gd tagging efficiency of at least 67% has been 
demonstrated [209], but it can likely be improved beyond that. We will discuss 
the impact of Gd tagging in the background reduction for DSNB in Sec. 5. The 
next-generation anticipated large WC detector is Hyper-Kamiokande, with half a 
Mton of mass [210]. 

An important special case of WC detector for SN detection is the “long-string” 
type. This type of detector comprises a very large volume of water or ice (part of a 
natural body), instrumented with strings of photomultiplier tubes. These kinds of 
detectors are nominally intended for ultra-high-energy neutrino detection (>GeV), 
and they are too sparsely instrumented to resolve individual SN neutrino inter¬ 
actions. However, they are still able to observe a SN neutrino burst as a diffuse 
glow of Cherenkov photons in the ice or water, via coincident increase in count 
rate of all the photomultiplier tubes within the time window of the burst [211,212]. 
Although event-by-event energy, angle and flavor information is not available, the 
photon statistics are high enough to provide good time resolution for the overall 
time structure of the burst [213]. Instances of this class are IceCube, which has a 
SN trigger; in its complete configuration, IceCube has 5160 optical modules [214] 
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and about 3 Mton effective detection mass, representing the largest current detec¬ 
tor for SN neutrinos. Others, which are noisier [215], may still retain information. 
Increased photomultiplier tube density, such as that proposed for the PINGU in¬ 
fill [216], will improve SN sensitivity by enhancing the probability to for multiple 
photons per interaction to be recorded. 

• Liquid scintillator: These detectors employ organic hydrocarbons in liquid form. 
Like water, the material has a large fraction of free protons. Scintillator detectors 
are therefore also primarily sensitive to IBD and hence to v e . Scintillators emit 
photons in response to energy loss of charged particles. Rather large light yields 
per energy loss (^50 times higher than light yields for Cherenkov radiation) can 
be obtained, resulting in good event-by-event energy resolution and potentially low 
energy thresholds. Typically neutrons can be detected via capture on free pro¬ 
tons with high efficiency; gadolinium doping is often employed to enhance capture 
rate and deexcitation gamma energy release. In order to take advantage of low 
energy thresholds, scintillator detectors must tackle the difficulties of radioactive 
contamination which are typically fierce below a few MeV. 

The currently-running instances of large scintillator detectors are LVD [23,217,218], 
Borexino [219,220], and KamLAND [221]; the SNO+ [222] experiment should turn 
on in the near future. A future planned 20-kton detector is JUNO [223,224]; an¬ 
other well-developed concept is the 50-kton LENA detector [225]. Surface reactor- 
neutrino scintillator detectors (all Gd-loaded), although relatively small (tens of 
ton scale) and unprotected by significant overburden, also represent a SN burst 
detection opportunity. Daya Bay in particular can defeat cosmogenic background 
by high IBD selection efficiency and by the requirement of coincidence between 
multiple, spatially-separated detectors [226]. Because scintillation light is emitted 
isotropically, the direction sensitivity of these detectors is limited (although see 
Sec. 3*3.2). Thanks to the low energy threshold, scintillator also has sensitivity to 
NC elastic scattering of neutrinos on protons. “Quenching” of scintillation light 
from heavy recoil products suppresses observability of this channel; nevertheless 
there should be some capability for existing and future detectors. Furthermore, 
information about the neutrino spectrum is available via the proton recoil spec¬ 
trum [199]. 

• Liquid argon (LAr): An extremely interesting possibility for the future is liquid 
argon time projection chamber (LAr TPC) technology. These detectors drift ion¬ 
ization charge from energy loss of charged particles in an electric field, and collect 
this charge on a planar anode (via different mechanisms- either collection or induc¬ 
tion wires, or multiplication at a liquid-gas interface). Two dimensions of the track 
can be reconstructed from collected-charge projection on the anode plane, and 
the third dimension is determined from drift at known constant velocity; there¬ 
fore three-dimensional tracks can be reconstructed. Products of low-energy SN 
neutrino interactions can be in general reconstructed as short (^10-cm scale) lep¬ 
ton tracks. In principle energy loss from Compton scatters from deexcitation and 
bremsstrahlung gamma rays can also be recorded, and used to tag interaction 
channels. The primary interaction from a SN neutrino burst will be v e absorption 
which gives liquid argon excellent sensitivity to the electron flavor component of 
the flux [227]. The detector also has intrinsic directional capability for anisotropic 
interaction products. Liquid argon time projection chamber technology has been 
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demonstrated for the Icarus [228] detector, as well as smaller prototypes such as 
ArgoNEUT (although with little emphasis on low-energy events). MicroBooNE, 
a 200-ton fiducial surface detector at Fermilab [229], will be online soon. The 
prospects for a large detector, DUNE, underground in the United States are very 
good [230]. 

• Lead: Another possible method of SN neutrino detection is to employ a heavy 
material, such as lead or iron. These materials have a high cross section for neu¬ 
trino interactions in the few tens of MeV range. In principle, detectors can observe 
energy loss of leptons, gammas, or other interaction products. However a poten¬ 
tially inexpensive mode of operation is to exploit the high probability for ejection 
of neutrons subsequent to a CC or NC neutrino interaction; counting single and 
double neutrons gives crude spectral information [231]. The ejected neutrons can 
be detected with dedicated neutron-detection technology. There is currently one 
operating instance of a lead-based SN neutrino detector, the HALO experiment at 
SNOLAB [232], which incorporates 76 ton of lead in conjunction with leftover 3 He 
neutron counters from the SNO experiment. The HALO-2 detector of kton scale 
is an envisioned future upgrade. 

• Dark matter detectors: Another category of detector will have some sensitivity to 
SN neutrinos via CEz/NS, the neutral-current interaction of a neutrino with an 
entire nucleus [203]. The cross section for this process is relatively high, produc¬ 
ing a handful of events per ton of detector material for a 10-kpc core collapse. 
However recoil energies are very low (few tens to hundreds of keV), making this 
channel experimentally challenging. Dark matter direction detection detectors are 
optimized to observe these low energy recoils and hence will have sensitivity to the 
total SN flux [233]. Current-generation detectors of less than ton scale will observe 
small statistics, but with potential scale-ups to multi-ton scales, this category of 
detectors will provide useful information about a SN burst. 

3*2. Detector summary and future prospects. - The possible distribution of core¬ 
collapse SNe in the Galaxy must follow the regions of star formation, notably in the 
spiral arms. The expected distance distribution for a simple model is shown in Fig. 16. 
One realizes that the distribution is very broad and that 10 kpc is probably a reasonable 
benchmark value. However, every distance between 2 and 20 kpc has similar probability. 
Therefore, one would expect a factor of ~100 for variations in the predicted neutrino 
event rate for a future Galactic SN explosion. For definiteness, Table I summarizes neu¬ 
trino event rates for current and future detectors assuming a typical SN at d = 10 kpc. 
Despite the scaling of the event rates with l/d 2 , most of the existing detectors would 
be massive enough to assure a sensitivity to a SN event in the whole Galaxy. Moreover, 
a 0.5 Mton water Cherenkov detector would collect ~ 10 events from an extragalactic 
SN in nearby galaxies at d ~ 1 Mpc, such as M31 and M33. The prospects for detect¬ 
ing mini-neutrino bursts from nearby galaxies is quite exciting, since about one SN per 
year is expected within distances up to 10 Mpc [235,236]. In this regard, the patient 
accumulation of data on an event-by-event basis from SN in nearby galaxies (up to 10 
Mpc) would be an additional possibility for study of the average supernova neutrino 
spectrum [235,236]. 

Figure 17 shows some examples of expected event rates as a function of time for several 
models in proposed future large detectors for a 11 M© and a 27 M 0 SN progenitors with 
LS EoS (see Sec. 2*4), and Fig. 18 shows the time profile for the same models observable 
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Fig. 16. - SN probability vs. distance from the Sun for a simple model of progenitor distribution 
(continuous curve). In comparison also the SN distribution for Type la SNe is shown (dotted 
curve). (Reprinted figure from [234]; copyright (2003) by the Institute Of Physics Publishing.) 


by IceCube. Among the current detectors IceCube is the one having the largest volume 
for SN neutrinos. Moreover, before two decades have passed we will possibly have a 
40-kton argon detector (DUNE), a 560-kton water detector (Hyper-K) and a 20-kton 
scintillator detector (JUNO). For water and scintillator, the event rates are dominated 
by z/ e , and for argon, the event rates are dominated by v e . This figure shows that 
all detectors can clearly distinguish the accretion phase, resulting in a “hump” in the 
neutrino signal, lasting ^ 0.5 s and the following cooling on a timescale of ~ 10 s. In the v e 
signal collected by the argon detector one would also clearly see the neutronization peak, 
lasting ~ 50 ms just after the core-bounce, unless suppressed by oscillations (see Sec. 4). 
We note that for large detectors, provision must be made when designing electronics 
and data acquisition to prevent data loss for the very high event rates expected for core 
collapses at < kpc distances. 

3*3. Neutrino astronomy and astrophysics. - The neutrino burst from a core collapse 
is also useful to astronomers and astrophysicists. In the following we present some of the 
applications discussed in literature. 

3*3.1. Finding the supernova: Early alert. Because neutrinos are released from the 
stellar collapse on a timescale of seconds, starting within milliseconds of infall, detection 
of a SN burst offers an opportunity to provide astronomers with an early alert of a 
SN occurrence. The relative delay between the neutrino wavefront reaching Earth and 
the first photons from shock breakout is expected to be on the timescale of hours ( 6 ). 
The Super Nova Early Warning System (SNEWS) [240, 241] is a world-wide network 
designed to alert astronomers, or any other interested scientists, of a burst. At the 
time of this review, the SNEWS network involves Super-Kamiokande, LYD, Borexino, 


( 6 ) The delay of three hours between the v e burst and the optical signal from SN 1987A implies 
that the velocity of electron antineutrinos is equal to that of light with an accuracy ^ 2 x 
10 -9 [237,238]. This bound represents an accurate test of special relativity. 
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Fig. 17. - Examples of total event rates per time bin for 27 M© and 11 M© SN progenitors (see 
Sec. 2'4), for water, argon and scintillator detectors, respectively from the top to the bottom 
panel, assuming a Galactic SN at d = 10 kpc. Core bounce is at 0.35 seconds. The dotted 
(dashed dotted) lines are rates for a simplified MSW oscillation assumption for the normal 
(inverted) hierarchy case (see Sec. 4'11). Note that expected events in non-equal logarithmic 
time bins are plotted, not rates. 
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Fig. 18. - Expected background-subtracted counts per time bin for IceCube, for IBD events 
(following [239]), for the same models as in Fig. 17. Error bars are statistical. Again note 
non-equal logarithmic time bins. 


KamLAND, IceCube and Daya Bay; sensitivity is to core collapses within the Milky 
Way. A central computer at Brookhaven National Lab (with backup at Bologna) receives 
datagrams sent by individual experiments and sends an automated alert to a mailing list 
if a coincidence is found within ten seconds according to the burst UT time stamps. 
The SNEWS network has been operational for an automated alert since 2005, although 
no coincidences have been detected. It is worth noting that gravitational waves also 
potentially provide a prompt signal of core collapse, and correlation with neutrino burst 
can potentially enhance the early alert [213,242]. Moreover, assuming that the angular 
location of the SN burst can be determined quickly, there is a very high probability for 
detection of the SN shock breakout in the infrared [243]. Such cross-correlations would 
be extremely useful against false triggers. 

It has been estimated that the silicon-burning phase preceding a SN explosion can 
release an energy of about 5 x 10 50 erg for a few days in neutrino-antineutrino pairs with 
average energy ^ 2 MeV [244,245]. IBD events from silicon burning have very low (below 
nominal detection threshold) positron energies, and the best prospects for detection are 
via neutron capture by adding Gd, provided that the events can be statistically distin¬ 
guished from background fluctuations. For a very nearby SN at a distance d = 1 kpc, one 
would expect ~ 500 events per day in a 0.5-Mton WC detector [246]. The kton liquid 
scintillator detector KamLAND will be able to detect pre-SN neutrinos for progenitors 
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Table I. - Current and proposed SN neutrino detectors as of the time of this writing. Event rate 
estimates are approximate for 10 kpc; note there may be significant variation by SN model. The 
“Flavors” column indicates the dominant flavor sensitivity (note that other flavor components 
may be detectable, with varied tagging quality). Not included are smaller detectors (e.g., reactor 
neutrino scintillator experiments) and detectors sensitive primarily to coherent elastic neutrino- 
nucleus scattering (e.g., WIMP dark matter search detectors). An asterisk indicates a surface 
detector, which may not be self-triggering due to background. Numbers in parentheses indicate 
that individual events will not be reconstructed; see text. 


Detector 

Type 

Mass (kt) 

Location 

Events 

Flavors 

Status 

Super-Kamiokande 

h 2 o 

32 

Japan 

7,000 

Ve 

Running 

LVD 

C n H 2n 

1 

Italy 

300 

Ve 

Running 

KamLAND 

C n H 2n 

1 

Japan 

300 

he 

Running 

Borexino 

On H 2n 

0.3 

Italy 

100 

h e 

Running 

IceCube 

Long string 

(600) 

South Pole 

(10 6 ) 

h e 

Running 

Baksan 

C n H 2n 

0.33 

Russia 

50 

h e 

Running 

MiniBooNE* 

On H 2n 

0.7 

USA 

200 

h e 

(Running) 

HALO 

Pb 

0.08 

Canada 

30 

Ve,V X 

Running 

Day a Bay 

On H 2n 

0.33 

China 

100 

h e 

Running 

NOz/A* 

On H 2n 

15 

USA 

4,000 

h e 

Turning on 

SNO+ 

C n H 2n 

0.8 

Canada 

300 

h e 

Near future 

MicroBooNE* 

Ar 

0.17 

USA 

17 

V e 

Near future 

DUNE 

Ar 

34 

USA 

3,000 

Ve 

Proposed 

Hyper-Kamiokande 

H 2 0 

560 

Japan 

110,000 

h e 

Proposed 

JUNO 

O n H 2n 

20 

China 

6000 

h e 

Proposed 

RENO-50 

C n H 2n 

18 

Korea 

5400 

h e 

Proposed 

LENA 

On H 2n 

50 

Europe 

15,000 

h e 

Proposed 

PINGU 

Long string 

(600) 

South Pole 

(10 6 ) 

h e 

Proposed 


at distances within d = 660 pc in the most optimistic scenario [247]. Detection prospects 
for detectors under construction or proposed, e.g., SNO+, JUNO, and LENA, look also 
excellent [248]. Such detection would make it possible to foresee the death of a massive 
star a few days before the stellar core collapse; eventually this could allow discrimination 
of the progenitor type [249]. Possible pre-supernova candidates that could explode at un¬ 
predictable future times include Betelgeuse, 3 Ceti, Antares, Epsilon Pegasi, Pi Puppis, 
NS Puppis, and Sigma Canis Majoris. 

3*3.2. Pointing to the supernova with neutrinos. It is of clear value to know the position 
in the sky of the SN, for several reasons. First, obviously astronomers need to know where 
to look— given that the visible SN may not show up in electromagnetic wavelengths for 
some time, a direction will aid in observing the very early stages. Furthermore, it is not 
clear that the core collapse always results in a bright SN, and direction information will 
help to potentially locate a weak explosion, or even a vanished progenitor. Knowledge 
of distance to the progenitor is valuable for precise determination of neutrino luminosity. 
Knowledge of the direction will aid as well in interpretation of the signal itself, in that 
the specific matter effects undergone by the neutrino flux will depend on the pathlength 
through the Earth [234]. 

The ability of a detector to determine the direction of a signal depends both on the in¬ 
trinsic anisotropy of the interaction, and on the capability of tracking the resulting prod¬ 
ucts. CC and NC interactions with nuclei, including IBD, have some energy-dependent 
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Fig. 19. - Angular distribution of IBD events (green) and ES events (blue) of one simulated 
SN. (Reprinted figure with permission from [253]; copyright (2003) by the American Physical 
Society.) 


anisotropy (e.g., [250]), but this is generally relatively weak for the most important in¬ 
teractions. Elastic scattering of electrons is however fairly strongly peaked forward. The 
angular distribution of IBD and ES events for a simulated SN is shown in Fig. 19. Of the 
interactions relevant for current detectors, argon and oxygen interactions do have some 
potentially observable anisotropy (e.g., [251]). 

Of the detector types currently available, Cherenkov and LAr tracking detectors are 
capable of determining the direction of charged particle tracks. Unfortunately scintilla¬ 
tor produces mostly isotropic light with little chance of determining scattered particle 
direction (some information may be available from reconstructed relative positions of 
IBD-produced positron and neutron energy loss [252]). Of detectors currently running, 
Super-Kamiokande has the only plausible chance of determining the SN direction using 
ES events (which should constitute a few percent of the observed events). The pointing 
accuracy in the worse case would be ^ 8° [253] for a SN at d = 10 kpc. This ability 
should improve with addition of Gd by reducing the near-isotropic background of IBD 
events, reaching ~ 3° [253]. Hyper-Kamiokande’s statistics would improve determination 
even further (with an accuracy as good as ^ 0.6°). LAr should also have decent pointing 
capabilities via lepton tracking. 

Another potential method for source direction determination is via triangulation, us¬ 
ing timing of several signals observed in different locations around the Earth [254,255]. 
This is challenging, given that the spread over time of the burst is longer than the travel 
time through the Earth; therefore, detection statistics of the burst need to be such that 
a sharp feature of the signal, such as the risetime or other feature, can be determined. 
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IceCube may be the best prospect for contributing to this kind of determination. With 
current detectors, triangulation in real time may be difficult, although with large future 
detectors prospects are more promising. Another possible pointing method is via the 
matter oscillation effect pattern, which might be possible to combine with timing infor¬ 
mation [256]. However this requires large statistics and good energy resolution, and is 
unlikely to do better than ES in a single directional detector. 

3*3.3. Probing supernova hydrodynamical instabilities through neutrinos. As discussed 
in Sec. 2, the first hydrodynamical 3D SN simulations with sophisticated neutrino trans¬ 
port have recently become available. The more massive cases of the simulated SN progen¬ 
itors (with 27 and 20 M© progenitor masses) show pronounced SASI phases interrupted 
by episodes of dominant convective overturn activity. The neutrino signal carries imprints 
of these hydrodynamical instabilities (see Fig. 5) and the detection of such modulations 
will therefore offer a unique chance to probe the core-collapse mechanism. It has been 
recently shown [91,92] that measurements of neutrinos from a future Galactic SN with 
neutrino Cherenkov telescopes (i.e., IceCube and Hyper-Kamiokande, see Sec. 3*2) will 
indeed be able to discriminate the SASI neutrino modulation. 

Figure 20 (top panel) shows the expected IceCube rate as a function of time for an 
observer located close to the SASI plane( 7 ) for a 3D SN simulation of a 27 M© star 
(see Sec. 2). Large-amplitude sinusoidal modulations of the IceCube event rate signal 
appear in correspondence to the first SASI episode at 120-260 ms. After a phase where 
convective motions dominate, a second SASI episode begins at about 410 ms. Given our 
scarce knowledge of neutrino flavor oscillations in the presence of SASI, two extreme cases 
have been plotted: One in which flavor conversions have been neglected and therefore the 
signal is caused by the non-oscillated v e (blue curve) and another one where complete 
flavor swap was assumed so that the signal is caused by v x (red curve). Both cases exhibit 
large-amplitude signal modulations with clear periodicity. 

The possibility to detect neutrino signal modulations depends on the viewing angle 
relative to the plane of the SASI sloshing or spiral motions. Observers located along the 
direction orthogonal to the plane of the first SASI episode will only detect very weak 
SASI modulations, as shown in the second panel of Fig. 20 (see also right panel of Fig. 5 
for an estimation of the favorable locations of the observer with respect to the SASI 
plane); observers placed on opposite directions away from the source will detect almost 
the same signal modulations but with opposite phases. 

The third panel of Fig. 20 shows the IceCube signal in 5 ms time bins, including a 
random shot noise realization, and the IceCube background fluctuations are plotted in 
black for comparison. For a SN up to 15 kpc, the SASI modulations in the neutrino signal 
will be clearly visible. The correspondent rate in a 560-kton water Cherenkov detector is 
plotted in the bottom panel of Fig. 20. A 560-kton water Cherenkov telescope will have 
no background in contrast to IceCube and it will be able to collect event-by-event energy 
information, but its expected rate should be about 1/5 of the IceCube one. Similar 
modulations of the neutrino signal due to SASI are also clearly detectable for the 20 M© 
SN progenitor [91,92]. 

The Fourier power spectrum of the IceCube rate exhibits a prominent peak at about 


( 7 ) Note that for all angular locations of the observer, the detected neutrino rate has been 
computed by integrating the neutrino emission emitted from the hemisphere facing the observer 
including projection effects associated with limb darkening. We refer the reader to Ref. [91] for 
details on the projection procedure. 
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80 Hz both for the 20 M 0 and the 27 M 0 SN progenitors. Such a frequency corresponds 
to the sloshing frequency of the shock front and is a function of the neutron star radius 
(similar for both progenitors as they were modeled with the same EoS) and of the shock 
radius (again comparable for both progenitors as their mass-infall history in the collapsing 
stellar core is fairly similar) [92]. 

In synthesis, the observation of SASI modulations by Cherenkov neutrino telescopes 
will provide a confirmation of our current theories of stellar core collapse. However, 
the possibility to discriminate SASI signatures will be dependent on the location of the 
observer relative to the plane where SASI sloshing or spiral motions develop. 

3*4. Outlook. - We are reasonably well prepared world-wide for the next SN burst, 
with a number of detectors online, of which Super-Kamiokande will give the largest event 
signal, and potentially provide pointing information. Additional scintillator detectors, 
some surface detectors, and a small lead detector, will enhance the yield, and IceCube 
will provide excellent burst time structure information. However, the current sensitivity 
is overwhelming for the v e component of the flux. The planned next generation of very 
large detectors will provide even more abundant v e statistics and, in addition, should 
also have broader flavor sensitivity via v e detection in argon, and high-statistics NC 
channels in water and scintillator. We look forward to recording the next burst with 
these enhanced capabilities. 

A high-statistics detection of a Galactic SN v burst would allow accurate studies of 
the signal in energy and time domains and in different interaction channels. In particular, 
we have shown how the observable SN v light curve would carry interesting information 
concerning the explosion mechanism and the hydrodynamical instabilities occurring in 
the deepest SN regions. Moreover, by reconstructing the duration of the different phases 
of the v burst one could have a direct empirical test for the predictions of the SN sim¬ 
ulations. Furthermore, the duration of the neutrino signal is an important diagnostic 
parameter to constrain exotic energy-loss processes in a SN, associated with the emission 
of novel light weakly-interacting particles [17]. These would shorten the duration of the 
observable SN neutrino signal. Finally, the SN v light curves and energy spectra would 
also carry intriguing signatures associated with the flavor conversions occurring in the 
deepest stellar regions, as we will discuss in the next Section. 
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Fig. 20. - Detection rate in IceCube (three upper panels) and in a 560-kton WC detector (bottom 
panel) for a 27 M® SN progenitor at 10 kpc for v e (no flavor conversions) and u x (complete 
flavor swap). The observer direction is chosen close to the SASI plane during the first SASI 
episode (120-260 ms), where a strong modulation of the neutrino signal is expected, except 
for the second panel where a direction of minimal modulation of the neutrino signal is shown. 
The lower panels include a random shot-noise realisation in 5 ms bins; the IceCube background 
fluctuations in the absence of a SN are shown in black. (Figure adapted from [92].) 
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4. — Supernova Neutrino Flavor Conversions 

Authors: A. Mirizzi, N. Saviano, I. Tamborra, S. Chakraborty 

In this Section, we will present an overview of the neutrino flavor transitions in core¬ 
collapse SNe. At first we introduce the equations of motion (EoMs) for SN neutrinos in 
the formalism of the density matrix, describing the vacuum, matter and v-v refractive 
terms entering the neutrino Hamiltonian. Then the phenomenology of self-induced flavor 
conversions, associated with the large v-v interaction potential, is illustrated through 
representative examples of its rich phenomenology. The Mikheyev-Smirnov-Wolfenstein 
(MSW) matter effects are also discussed as well as the impact of the SN shock-wave 
propagation and of matter turbulence on the neutrino conversion probabilities. Finally, 
observable signatures of the SN neutrino oscillations imprinted on the neutrino burst are 
characterized. Particular attention is devoted to the sensitivity of these observables to 
the unknown neutrino mass hierarchy. 

4T. Three-neutrino oscillation framework. - There is compelling experimental evi¬ 
dence [257] that the three known active neutrino states with definite flavor v a (a = e, /i, r) 
are linear combinations of the mass eigenstates Vi (i = 1, 2, 3), and that the Hamiltonian 
of neutrino propagation in vacuum [258-260] and matter [261] does not commute with fla¬ 
vor. The evidence for flavor non-conservation (i.e., “neutrino oscillations”) comes from a 
series of experiments performed in about four decades of research with very different neu¬ 
trino beams and detection techniques. Namely, the solar neutrino [262-264] experiments: 
Homestake [265], SAGE [266], GALLEX-GNO [267,268], Kamiokande [269], Super- 
Kamiokande (SK) [270,271], the Sudbury Neutrino Observatory (SNO) [272-274], and 
Borexino [275]; the long-baseline reactor-neutrino experiment KamLAND [276,277]; the 
atmospheric neutrino experiments like Super-Kamiokande [278-280], MACRO [281], and 
Soudan-2 [282]; the long-baseline accelerator neutrino experiments: KEK-to-Kamioka 
(K2K) [283], MINOS [284], OPERA [285] and Tokai-to-Kamioka (T2K) [286,287]. 

Except for yet controversial results from short-baseline neutrino experiments, that 
point towards the existence of extra sterile neutrino states [with m ~ 0(1) eV] [288,289], 
all the data from the above experiments are consistent with the simplest extension of 
the standard electroweak model needed to accommodate nonzero neutrino masses and 
mixings. Namely, with a scenario where the three known flavor eigenstates are 

mixed with the three mass eigenstates ^ 1 , 2 , 3 , through a unitary matrix 1/, which in 
terms of one-particle neutrino states \v), is defined as (see, e.g., [257]) 

( 16 ) K ) = X^K>- 

i= 1 

A common parametrization for the matrix U is [290]: 

/ 1 0 0 \ / ci 3 0 si 3 e~ u \ / C 12 S 12 0 \ 

(17) U = 0 C 23 523 0 1 0 -S12 C 12 0 

\ 0 —S23 C23 ) \ — S\3e l5 0 Ci 3 / \ 0 0 1 ) 

with Cij = cos Oij and Sij = sin 0^-, Oij being the mixing angles and 5 E [0, 2tt] being the 
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CP-violating phase ( 8 ). 

The current neutrino phenomenology implies that the three-neutrino mass spectrum 
{rrii}i =j 2,3 is made by a “doublet” of relatively close states and by a third “lone” neutrino 
state, which may be either heavier than the doublet (“normal hierarchy,” NH) or lighter 
(“inverted hierarchy,” IH). Typically, the lightest (heaviest) neutrino in the doublet is 
called v\ (y^) and the corresponding mass squared difference is defined as 

(18) Sm 2 = m 2 — m 2 > 0 

by convention. The Sm 2 is traditionally named the solar mass squared difference. The 
lone state is then labeled as z/ 3 , and the physical sign of m 2 — m 2 2 distinguishes NH from 
IH. The second independent squared mass difference Am 2 , called also atmospheric mass 
squared difference, is 


(19) 


Am 2 


so that the two hierarchies (NH and IH) are simply related by the transformation 
+Am 2 —)> — Am 2 . Numerically, it results that Sm 2 <C Am 2 . 

The latest solar, reactor and long-baseline neutrino oscillation analyses indicate the 
following ±2a ranges for each parameter (95% C.L.) taken from [291] (see also [292,293] 
for other updated analyses) 


( 20 ) 


Am 2 = (2.43±g;il) x 10“ 3 eV 2 , 
Sm 2 = (7.54±g;||) x 10“ 5 eV 2 , 
sin 2 0 i 2 = ( 3 .O 8 S 0 . 34 ) x 10 \ 
sin 2 0 23 = (4.37±“ 3 ) x 1CT 1 , 
sin 2 0 13 = (2.34tg;$) x 1(T 2 . 


Remarkably, the 1 — 3 mixing angle 6^13 has been the last one measured in 2012 by 
the Daya Bay [294] and RENO [295] reactor experiments. These recent measurements 
confirmed and greatly strengthened the significance of early hints suggested by the long- 
baseline v^-Ve experiments T 2 K [296] and MINOS [297] as well as by the Double Chooz 
reactor experiment [298], especially when analyzed in combination with other oscillation 
data [299,300]. The measurement of a “large” value of 0 13 has siginifcantly reduced the 
ambiguity in characterizing the oscillated SN neutrino oscillations. Notably, the latest 
global analysis hints towards a non-zero CP-violation around S ~ IAtt at > la level [291], 
see also [292,293] ( 9 ). 

Within the three-neutrino scenario, the most important unsolved problems require 
probing the CP-violating phase S , the mass hierarchy, and the absolute neutrino masses. 


( 8 ) We neglect here extra-phases possible if neutrinos are Majorana particles, since they are 
not relevant in oscillations. 

( 9 ) It has been shown that in the context of SN neutrino oscillations, CP-violation effects are 
negligible [301-303]. Therefore, we will not include them in the following discussion. 
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In this context, a high-statistics detection of Galactic SN neutrinos could give a unique 
help to solve some of these open issues. 

4‘2. Equations of motion for supernova neutrinos. - The treatment of neutrino mixing 
is well understood in terms of the propagation of a beam of particles in vacuum or in a 
medium, and neutrino flavor transitions have been described by means of a Schrodinger- 
like equation [304]. Neutrinos free streaming beyond the neutrinosphere also interact 
among themselves ( neutrino self-interactions). As pointed out in seminal papers [1- 
3,167,305-316], in the deepest regions of a SN (as well as in the Early Universe) the 
neutrino gas is so dense that neutrinos themselves form a background medium leading 
to intriguing non-linear effects in the flavor distribution. Neutrino-neutrino interactions 
could trigger large self-induced flavor conversions in the deepest SN regions [5,6]. Under 
these circumstances, neutrinos emitted with different energies would be locked to oscillate 
in a collective fashion. 

A natural treatment to properly characterize the flavor evolution in this situation 
requires the formalism of the neutrino density matrix [317] 


( 21 ) 


£?p,x — 


Qee Qen Qer 
Qlie Qfir 

Qre Qt/i Qtt 


The diagonal entries of this matrix are the usual occupation numbers, whereas the off- 
diagonal terms encode phase information related to the oscillations. An analogous ex¬ 
pression holds for the antineutrino density matrix, p p x . See also [318-322] for similar 
derivations and [323-326] for recent developments on the formalism. 

In SNe one is concerned with the spatial evolution of the neutrino fluxes in a quasi¬ 
stationary situation. Therefore, the matrices p P)X do not explicitly depend on time, so 
that the EoMs reduce to the Liouville term involving only spatial derivatives [317,324] 


Vp • V x £> PjX = -i[ftp ac , <?p, x ] - £p.*] > 


( 22 ) 


v p ' Vx^ 


P-X 


T i [D. 


vac - i 

P 5 Cp,xJ 


■ , £* P , X 


In the flavor basis the first term in the Hamiltonian at the right-hand-side (r.h.s.) of 
Eq. (22) for ultrarelativistic neutrinos represents the vacuum oscillation term, 


(23) = U^-U^ . 

2p 

The U matrix describes the mixing [Eq. (17)] and M 2 is the squared mass matrix which, 
except for a common term proportional to the identity matrix and irrelevant for os¬ 
cillations, is parametrized, in terms of the solar Sm 2 [Eq. (18)] and atmospheric Am 2 
[Eq. (19)] mass-squared differences, as [290] 

(24) M 2 = diag (ml,ml,ml) = ,+^-,±Am 2 ^ . 

The sign db in front of Am 2 refers to NH (+) and IH (-), respectively. One can associate 
two vacuum oscillation frequencies to these mass-squared differences, i.e., 


(25) 




Sm 2 

~2E 


and ujh 


Am 2 
2 E 
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The vacuum energy can be affected by a shift (similar to the photon index of refrac¬ 
tion) when neutrinos propagate in a medium. This effect is induced by their forward 
scatterings with the medium constituents [304]. The refractive effect for SN neutrinos is 
described by the second term in the Hamiltonian at the r.h.s. of Eq. (22) 

(26) ^P,X = ^xL + a/2Gf J ( 2^)3 (Tu> x _ £q,x) (1 — v q ' v p) • 


The first term at the r.h.s. in Eq. (26) represents the ordinary refractive matter effect. In 
particular, for typical energies of SN neutrinos [E ~ 0(10) MeV] the only relevant process 
is due to charged-current interactions of electron neutrinos v e (or antineutrinos D e ) with 
the background e± (the neutral current is flavor conserving and therefore equal for all 
neutrino flavors). This is the well-known MSW effect, first pointed out by Wolfenstein, 
and by Mikheyev and Smirnov [261]. This term is represented by A x = — n e +) 

with L = diag(l, 0, 0) in the weak interaction basis. 

Neutrino-neutrino interactions, dominant in the deepest SN regions, make an addi¬ 
tional contribution to the refractive energy shift, represented by the second term at the 
r.h.s. of the Eq. (26). This term is proportional to the neutrino density matrix p PjX 
that in the presence of mixing has also off-diagonal elements, giving rise to “off-diagonal 
refractive indices” as first pointed out by Pantaleone [1-3]. The main complication in the 
v-v refractive term for SN neutrinos is the angular factor (1 — v q • v p ) = (1 — cos# pq ) 
coming from the current-current nature of the weak interactions, where v p = p/p is the 
neutrino velocity. This angular term averages to zero for an isotropic v gas. However, for 
the non-isotropic neutrino emission from the SN core this velocity-dependent term would 
not average to zero, producing a net current that leads to a different refractive index for 
neutrinos that propagate on different trajectories. This is at the origin of the so-called 
“multi-angle effects” [6,327-329]. Remarkably, while in an isotropic neutrino gas the 
self-induced effects would lead to a collective behavior in the flavor evolution [7], in an 
anisotropic case this is not guaranteed. Indeed the multi-angle term in some cases chal¬ 
lenges the collective behavior of the flavor evolution, leading to flavor decoherence under 
certain conditions, with a resultant flux equilibration among electron and non-electron 
(anti)neutrino species [327-329]. 

The multi-angle flavor evolution described by the partial differential equations, Eq. (22), 
has not been solved in its full complexity until now. However, a few years ago the first 
large-scale multi-angle simulations were developed within the so-called “bulb model” [5, 
8,328], whose geometry is represented in Fig. 21. In this framework it is assumed that 
neutrinos are emitted uniformly and half-isotropically (i.e., with all the outward-going 
modes occupied, and all the backward going modes empty) from the surface of a spherical 
neutrinosphere, like in a blackbody. Moreover, it is assumed that there is azimuthal sym¬ 
metry of the neutrino emission at the neutrinosphere and that all the physical conditions 
in the star only depend on the distance r from the center of the star. 

Under the assumption of azimuthal symmetry, one can reduce the general partial dif¬ 
ferential equations [Eq. (22)] into ordinary differential equations, projecting the evolution 
along the radial direction. Then 


(27) 


„ d 

'p-V x -MV dr , 

p ->■ (E ~ |p|, u e= sin 2 Or) , 

£p,x ^ Qp,u,r i 
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Fig. 21. - Geometrical picture of the neutrino bulb model. Courtesy of A. Marrone. 


where Or is the emission angle relative to the radial direction at the neutrinosphere, 
r = R. 

The distribution matrices Q p ,o, r ( or Qp,u,r) are n °t especially useful to describe a spher¬ 
ically symmetric system because they vary with radius even in absence of oscillations. A 
quantity that is conserved in absence of oscillations is the total flux matrix [328] 


(28) 


J r = Airr A 


/ 


d 3 p 

j2^)3 9p ’ 9 ’ rVr 


that in terms of g p , u ,r reads 


( 29) Jr = r d j\ui pur = . 

J 0 Jo ’ Jo Jo 2 7T 

This implies that the trace of J p?n?r is r-independent. For the sake of brevity, we address 
the interested reader to [328] for the explicit form of the EoMs for the J Pjn?r matrices. 

4*3. Polarization vectors and equations of motion for a 2u system. - Often SN neutrino 
flavor evolution can be characterized in a 2v flavor scenario, associated with Am 2 and 
# 13 , in which the electron state v e mixes with the non-electron one v x . We will comment 
later about possible three-flavor effects associated with Sm 2 and #12 [330,331]. 

In the 2v framework, one can expand the 2x2 matrices in the EoMs in terms of the 
unit matrix I and the Pauli matrices a. One can introduce polarization vectors P P:U ,r 
from the expansion of the flux matrices J P;U?r . One of the possible definitions suitable in 
the SN context is the following [328] 

1 _ F(pe) P ,u T F{yx) P , u . F{ye) F'ipx) -p 

Jp,u,r — ^ T ~ *^p,u,r ' & 1 
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(30) 


F(v e )p,u + F(y x ) 


p,u 


J p,u,r 


F(v e ) - F(v x )- 

2 *^p,u,r ' & 


where the quantities F(u a ) p ^ u are the number flux of u a emitted from the neutrinosphere 
with energy E = |p| and emission angle u = sin 2 Or. The total number flux of flavor a is 
a simple integral F(v a ) = f dp f duF(i/ a ) p , u . Note that the normalization of the second 
term in Eq. (30) is the same for v and v and it is the difference between the z/ e ’s and v x s 
fluxes. We factorize the flux of each flavor as F(is a ) p , u = N Ua x fv a {p) x <Pv a (u). The 
v number N V(x — L Ua / (E Ucx ) is expressed in terms of the v luminosity L V(x and of the 
v average energy (E Ua ) of the different species. The function f Ua (p) is the normalized 
neutrino energy spectrum (f dpf Ua (p ) = 1). An example of a quasi-thermal energy 
spectrum, widely used in SN literature is given in Eq. (9). Finally, p V(x {u) represents the 
v angular distribution that we assume flat in u and flavor and energy independent, i.e. 
(u) = 1. 

One can expand the other quantities in Eqs. (22)-(23)-(26) as 

^p aC = +cj p B • a) , 

(31) L = | (I + L - cr) , 

where ujq = ra 2 2 +u7j, and the vacuum oscillations are determined by the mass difference 
Am 2 and vacuum mixing angle #13 [7]: 

uj p = Am 2 /2p , 

(32) B = (sin26>i3,0,=F cos 26 >i 3 ) , 


with the upper (lower) sign refers to normal (inverted) mass hierarchy. For small # 13 , 
B =pz. Moreover, in an ordinary medium composed by electrons and positrons, L = e z . 

Using the definition of Eqs. (30)-(32) one can write the EoMs [Eq. ( 22 )] for the 
polarization vectors [328] 


8 P 

u r r p,u,r 


- (cjpB T A r L) 


+ Mr 


J dq J du r {\ 


q,u,r A q 


Pq,u',r) 


1 


^r,u^r,u' 


8 r F P: u, r 


-(— cjpB T A r L) 


(33) + fi r J dq J du'(Pg tU > tr 


Pg,n',r) 


1 


^r,u^r,u' 


- 1 


- 1 


X P r 


X P 


p,u,r 5 


where we introduced the self-interaction potential 


(34) 


= V2 G, 


F(v e ) ~ Fjvx) 

47 r r 2 


and the radial neutrino velocity 


( 35 ) 
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The initial conditions for the polarization vectors at the neutrinosphere are the following 


(36) 


p;,u( R ) = 

Ku( r ) = 


F iye)p,u F (h'x'jpu 
Pipe) ~ F{P X ) 
F{ye)p, U F {yx)p,u 

F(p e ) ~ F(P X ) 


The presence of the trajectory-dependent factor u n?r in the neutrino-neutrino interac¬ 
tion term in Eq. (33) leads to the multi-angle effects in the self-induced oscillations. In 
this regard, it has been understood [327] that the flavor asymmetry is crucial to assess 
the impact of multi-angle effects. In particular, in the presence of a sufficiently large flux 
hierarchy among different species, multi-angle decoherence would be suppressed. 

4*3.1. Matter suppression of self-induced flavor conversions . Remarkably also the mat¬ 
ter term in Eq. (33) is affected by multi-angle effects due to the factor u n?r . It has been 
shown in [332] that matter effects play a sub-dominant role in the development of self- 
induced flavor conversions when 


(37) n e - - n e + < n„ e - n Dx . 

Conversely, if n e - — n e + n„ e — n „ x , the multi-angle matter effects could produce a 
large spread in the oscillation frequencies for neutrinos travelling on different trajectories, 
blocking the self-induced flavor conversions. Finally, when n e - — n e + ~ n„ e — np x , 
matter-induced multi-angle decoherence in the neutrino ensemble may occur. According 
to realistic SN models, the matter density is expected to dominate over the neutrino one 
during the early accretion phase (see Sec. 4*4). In this situation, dedicated studies have 
been performed in [333-338] using inputs from recent hydrodynamic SN simulations. 
It has been found that the large matter term inhibits the development of collective 
oscillations at early times for iron-core SNe. Only in the case of low-mass O-Ne-Mg SNe, 
where the matter density is never larger than the neutrino one, the matter suppression is 
not complete and partial self-induced flavor conversions are possible at early times [333, 
334,338]. 

The matter suppression of self-induced effects for iron-core SNe makes relatively easy 
the prediction of the oscillation effects during the accretion phase, since the neutrino 
fluxes will be processed by the only MSW matter effects. As a consequence, as we will 
discuss in Sec. 4*11, the detection of the SN neutrino signal at early times is particu¬ 
larly relevant for the mass hierarchy discrimination. Nevertheless, for all types of SN 
progenitors self-induced effects still remain crucial during the later cooling phase, when 
the matter density decreases continuously (at post-bounce times t p b > 1 s) and becomes 
sub-dominant with respect to the neutrino density. 

As final caveat, we remark that the matter suppression of collective oscillations studied 
in [333-338] has been characterized referring to spherically symmetric one-dimensional 
SN models. However, as discussed in Sec. 2, large deviations from a spherical neu¬ 
trino emission can be generated by hydrodynamical instabilities such as SASI or LESA 
especially during the accretion phase. Before the explosion sets in, during the standing- 
accretion-shock phase, the recently discovered LESA instability [102] is responsible for 
a neutrino lepton number flux emerging primarily in one hemisphere (see Sec. 2). The 
asymmetry between v e and P e can be very small and even negative in some regions of the 
star, becoming potentially responsible for flavor instabilities [328]. LESA also influences 
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the electron density and therefore the matter background felt by neutrinos. Moreover, 
the neutrino flavor ratio can also be affected by SASI that seems to coexist with LESA 
in more massive SN progenitors (see Sec. 2) and that is responsible for wild shock os¬ 
cillations and a time-dependent directional bias of the neutrino emission. Exploratory 
work conducted within a simplified setup [339] confirmed the suppression of flavor con¬ 
versions within the shock front radius in the presence of LESA. However, such results 
need to be tested within a more realistic scenario and in the presence of SASI. Besides 
the relevance for phenomenological purposes, the confirmation of the possible suppres¬ 
sion of self-induced flavor conversions will also be crucial for the modeling of the neutrino 
propagation in the hydro dynamical SN simulations. 

4*4. Supernova potential profiles. - Before discussing the solution of the neutrino EoMs 
for the dense SN neutrino gas, we find it worthwhile to show snapshots of the interaction 
potentials appearing in the r.h.s. of Eq. (33) in Fig. 22. We use the data from a 27 
M 0 SN simulation for different post-bounce times (see Sec. 2 ). As from Eq. (34), the 
neutrino-neutrino potential is 


(38) 


= V2G, 


F(v e ) ~ F{vx) 

47r r 2 


1 ( Ly e _ L Vx 

4jrr 2 \(E Pe ) (E Px ) 


= 7.0 x 10 5 km -1 (-E±- 
V \^e / 


L Px \ 15 MeV 
(E Px )) 10 52 erg/s 


10 km^j 


where the number fluxes of the different species is a are expressed in terms of the neutrino 
luminosities L Va and of the average energies (E Ua ). As shown in Fig. 22, the v-v potential 
decreases as the SN cools. 

Matter effects on SN neutrinos [Eq. (33)] depend on the potential 


(39) 


/ Y 

X r = V2GFn e (r) = 1.9 x 10 6 km -1 x ( — 

U.O 


10 10 g/cr 


encoding the net electron density n e = n e - — n e +, where Y e = Y e - — Y e + is the net 
electron fraction and p is the matter density. The numerical values of p r and A r from 
the previous two expressions are quoted in km -1 . 

Figure 22 shows that the SN electron density profile is non-monotonic, time-dependent, 
and presents an abrupt discontinuity, corresponding to the position of the forward shock- 
front. In fact, the shock-wave, while propagating outwards at supersonic speed, leaves 
behind a rarefaction zone, and creates a high-density front with a sharp density drop 
(down to the static value) at ( 10 ) 

(40) rf w d = forward shock radius . 

We also note that additional features can appear such as the reverse shock propaga¬ 
tion [340]. For example, after the explosion, a neutrino-driven baryonic wind can develop 


( 10 ) Note that in hydrodynamical simulations, like the ones shown in Fig. 22, the forward shock- 
front is typically softened because of the limited numerical resolution. Therefore it has to be 
steepened by hand when adopted for studying the neutrino flavor evolution. 
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Fig. 22. - Snapshots of SN potentials for different post-bounce times (1.0 —7.0 s) for a 27 M© SN 
progenitor (see Sec. 2). The profile at 0.2 s is an illustrative case for a typical condition before 
shock revival. The matter potential A r is drawn with thin curves, while the neutrino potential 
Hr with thick ones. The horizontal bands represent the vacuum oscillation frequencies relevant 
for the MSW resonant conversions associated with Am 2 (uh) and 8 m 2 (cjl), respectively (see 
the text for details). 


and collide with the (slower) SN ejecta, thus triggering a second (reverse) shock at 

(41) r rev = reverse shock radius , 

which propagates (at lower velocity) behind the forward one [340]. Neutrinos may thus 
encounter two subsequent density discontinuities, leading to significantly different spec¬ 
tral features with respect to the case of a single discontinuity. 

One expects that the matter term would lead to resonant flavor conversions via the 
MSW effect [261] when the matter potential is of the order of the vacuum oscillation 
frequencies [Eq. (25)] c jh,l [341], i.e. 

(42) A r ~ wl,h • 

These oscillation frequencies are represented by the two horizontal strips in Fig. 22 for 
a neutrino energy range E G [5 — 50] MeV. Note that resonant flavor conversions should 
be expected for r G [10 3 ,10 5 ] km. 

Comparing the neutrino and matter density profiles, we realize that in the deepest 
SN regions n v n e (r < 10 3 km, see t = 1.0 s in Fig. 22), except during the accretion 
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phase (t = 0.2 s in Fig. 22). When the neutrino density dominates over the matter 
one (as during the cooling phase), self-induced flavor conversions would develop without 
any hindrance. This schematic investigation suggests that intriguing effects should be 
expected in the SN neutrino flavor evolution, with an interplay between self-induced and 
matter terms. Moreover, self-induced oscillations usually develop at lower r than the 
MSW flavor conversions, so that their effects effects could be studied separately for this 
progenitor. 

4*5. Self-induced flavor conversions: Single-angle approximation. - In order to simplify 
the numerical treatment of the EoMs [Eq. (33)] and to develop analytical interpretations 
of neutrino self-interactions, it has been often adopted in the literature the so-called 
“single-angle” approximation [5,8,328]. The main idea is to assume a single angular 
mode representative of all the neutrino ensemble. For a blackbody emission, in which 
all of the angular modes are equally occupied, it is natural to take neutrinos emitted 
at Or = 7r/4 (i.e., uq = 0.5) as the representative ones ( n ). In this case, there is no 
integration over u and all the radial velocities appearing in the EoMs are set to 


(43) 


^ r,uo 



The EoMs in the single-angle approximation read 


(44) 


d r P p — + A r L + /i*D] x Pp , 

<9 r P p = [—cjpB T A r L T yU r D] x Pp , 


where the self-interaction term only depends on the radial neutrino self-interaction strength 

F(v e ) - F(v x ) R 4 


(45) 


Mi 


= V2Gi 


4 t tR 2 


1 


R 2 


1 


— Mr 


2r 2 l-i£ 


that at large distances from the core declines as p,* ~ i7 4 /2r 4 , and on the difference of 
the total polarization vectors D = P — P, where P = f dpP p (and analogously for P). 
In order to simplify the notation we remove the subscript r in the polarization vectors. 

Assuming that self-induced effects are not suppressed by the multi-angle effects asso¬ 
ciated with a dominant matter term, it has been shown that in the single-angle approx¬ 
imation, when A r c j p , the EoMs can be studied in a co-rotating frame with angular 
velocity A r L [7-9]. After this rotation, the only notable matter effect would be the 
suppression of the in-medium mixing angle with respect to the vacuum one. 

The EoMs [Eq. (44)] for A r = 0 imply conservation of the scalar product 


(46) D B = const = 0*6- • z = Ffr - = Te , 

Fp e ~ Fv x 

where e indicates the flavor asymmetry. Equation (46) corresponds to the conservation of 
the lepton number [7], and implies pair conversions v e V e —x v x V x [7]. In the co-rotating 
frame D • L is the corresponding conserved quantity. 


62 



Supernova Neutrinos: Production, Oscillations and Detection 


Neutrinos 



Antineutrinos 



E(MeV) 


E(MeV) 


Fig. 23. - Initial SN neutrino and antineutrino fluxes with F® e : F® e : F® x — 2.40 : 1.60 : 1.0. 
Neutrinos are plotted in the panel on the left, while antineutrinos are on the right panel. Electron 
species (e) are plotted with continuous curves, while non-electron species (x = /i, r) with dashed 
curves. 


4*6. Synchronized vs. bipolar oscillations and spectral splits. - We will start to char¬ 
acterize the SN neutrino flavor conversion phenomenology, first discussing self-induced 
effects associated with the neutrino-neutrino interactions in the single-angle approxima¬ 
tion [Eq. (44)]. Assuming to be in a co-rotating frame such that A r = 0 (see Sec. 4*5), 
the EoMs for the flavor evolution of the dense SN neutrino gas, even in the single-angle 
approximation, present a rich phenomenology since neutrino-neutrino interactions are 
strongly dependent on the ordering and hierarchy of the SN neutrino fluxes of the differ¬ 
ent species. Since the ordering of the SN fluxes changes significantly during the different 
post-bounce phases, one should expect a rich phenomenology as a function of the time 
after the bounce [5,8,331,343-347]. 

In this Section, we take as illustrative example a “classical” case, where the flavor 
dynamics of the polarization vectors under the self-induced effects is analytically well- 
understood through an analogy with a gyroscopic pendulum in flavor space [7-9,348]. 
This behavior is realized for cases presenting a large excess of electron neutrinos over 
non-electron species, i.e. F® e F® e F® x (as expected, e.g., during the accretion 
phase). In the pendulum analogy, the generic motion of the polarization vectors in flavor 
space is a combination of nutation and precession. The neutrino-antineutrino asym¬ 
metry e [Eq. (46)] acts as a “spin” for the system with a large asymmetry inducing a 
precession-like motion around Co> p B, leaving only the nutation motion when the asymme¬ 
try is vanishing [7]. The role of the inertia moment in the pendulum analogy is played 
by the inverse strength of neutrino-neutrino interaction fi* [Eq. (45)]. The neutrino mass 
hierarchy sets the behavior of the system. In NH, the pendulum starts in downward 
(stable) position and, assuming a small mixing angle, stays nearby, i.e. only manifests 
small flavor changes. Conversely, the pendulum starts in upward (unstable) position 
in IH—the misalignment being of 0(0 13 )—and it is subject to the maximal flavor re¬ 
versal pair-conversions v e v e -F v x v x [7], that conserve the lepton number [Eq. (46)]. 


) A formal derivation of the single-angle approximation has been presented in [342]. 
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Fig. 24. - Single-angle simulation in inverted mass hierarchy for an initial flux ordering F„ e : F® e : 
F® x = 2.40 : 1.60 : 1.0. Modulus (dashed curves) and z-component of the global polarization 
vector (continuous curves) of neutrinos (P) and antineutrinos (P), as a function of radius. 


It has been shown that large flavor asymmetry could block possible multi-angle deco¬ 
herence effects [328]. Therefore, the single-angle approximation is well justified in this 
case [8,349]. Figure 23 shows an example of this flux configuration, obtained considering 
representative neutrino fluxes with the following spectral ordering at the neutrinosphere: 
F° e : F? e : = 2.40 : 1.60 : 1.0. 

In order to illustrate the self-induced flavor dynamics, Fig. 24 shows the radial evo¬ 
lution of the modulus P = |P| and ^-component P z of the global neutrino polarization 
vector P (and analogously for the antineutrino vector P) in the nontrivial case of IH for 
the fluxes shown in Fig. 23 (no significant effect occurs in NH). The radial evolution of 
the polarization vectors can be interpreted as follows. Up to ^ 90 km, it is P = P z and 
P = P z : all polarization vectors are “glued” along the vertical axis. In the gyroscopic 
pendulum analogy, this corresponds to a precessing top with a huge spin and negligible 
displacement from the vertical 2 -axis. It just spins in the upward position without falling. 
This behavior is named synchronized oscillation regime [7,312]. At r ^ 90 km, the pen¬ 
dulum falls for the first time and nutations appear, marking the transition from synchro¬ 
nized to the so-called bipolar oscillation regime ( 12 ). The nutation amplitude gradually 
decreases as ~ /i* 1 ^ 2 being driven down by the decreasing v-v potential strength. In the 
spinning top analogy, the relaxation of the pendulum to its downward rest position as 
kinetic energy is extracted by the reduction of the neutrino-neutrino interaction poten- 


( 12 ) Notably, matter effects delay this transition by a few nutation periods with respect to what 
expected in the presence of v-v interactions only [7]. 
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tial and thus the increase of the pendulum moment of inertia. This would lead to large 
self-induced flavor transitions, occurring in the form of pair-conversions v e V e —> u x V x [7] 
that conserve the lepton number [Eq. (46)]. 

As a consequence of the self-induced flavor dynamics, antineutrinos tend to completely 
reverse the polarization vector (P —» — P), so that P z ~ — P asymptotically. Neutrinos 
also try to invert their global polarization vector. Then, |P| decreases for r > 150 km. 
However, the inversion cannot be complete, because of the lepton number conservation 
D z = P z — P z at any r [see Eq. (46)]. We remark that the inversion of the polarization 
vectors corresponds to large flavor conversions occurring at low-radii (r < 200 km) where 
the system would not have exhibited any flavor change in the presence of matter effects 
only. 

A generic outcome of the self-induced flavor conversions is the development of “spec¬ 
tral swaps” among the fluxes of different flavors, marked by “spectral splits” at the 
boundary features of each swap interval [344,346,350-356]. In order to elucidate these 
effects with the specific example we are discussing, we refer to the flavor evolution of the 
total polarization vectors P and P shown in Fig. 24, for the initial fluxes of Fig. 23. One 
would not expect complete flavor conversions in both neutrino and antineutrino sectors, 
due to the partial inversion of P. The initial spectra shown in Fig. 23 give a larger num¬ 
ber of v e over v e and v x . Therefore, if all antineutrinos are swapped, the conservation of 
the net lepton number [8] 

pOO poo 

(47) / dE(F Ve -F Vx )= dE(F Pe -F p J, 

♦'-E'split ^0 

prevents the same thing happening to the neutrino channel. This explains the appearance 
of the spectral neutrino splits in the oscillated fluxes. The previous equation allows us to 
determine the “spectral split” energy E sp n t separating the swapped part of the spectrum 
from the unswapped one. This is shown in Fig. 25, where we plot the initial energy spectra 
for v e (black dashed curves) and v x (light dashed curves) and after collective oscillations 
for v e (black continuous curves) and v x (light continuous curves) at r = 350 km. Neutrino 
oscillated spectra (left panel) clearly show the split effect and also the corresponding 
sudden swap of v e and v x fluxes above E sp u t ~ 10 MeV. The antineutrino spectra (right 
panel) are almost completely swapped with respect to the initial ones. However, they 
also present a spectral split at low-energy E sp u t ~ 3 MeV. This feature is not explained 
by the conservation law expressed by Eq. (47) and its nature is associated with non- 
adiabatic features of the flavor evolution, as discussed in [349]. We address the interested 
reader to [351-353] for a detailed theoretical explanation of the spectral splits through 
the adiabatic solution of the EoMs. 

4*7. Self-induced spectral splits and multi-angle effects. - For a long time the results 
shown in Fig. 25—an almost complete flavor exchange in the antineutrino sector and 
a spectral split in the neutrino channel in IH, no conversion in NH, and sub-leading 
multi-angle effects—have been considered the paradigm of the self-induced effects in SN 
neutrino flavor conversions. However, later it has been realized that this is not the most 
general situation [343]. Indeed, when z/s have a moderate flavor hierarchy of fluxes and 
spectral energies are not large, as expected during the cooling phase, more complicated 
conversion patterns could be associated with the self-induced effects [347]. Multiple 
spectral splits are possible in both NH and IH [343,347] for neutrinos and antineutrinos. 
Under such conditions, also multi-angle and three-flavor effects may play a crucial role. 


65 



Supernova Neutrinos: Production, Oscillations and Detection 



0 .. . 1 L/T. r-~rr^ l 1 

0 10 20 30 40 0 10 20 30 40 50 


E(MeV) E(MeV) 


Fig. 25. - Two-flavor single-angle simulation in inverted hierarchy for an initial flux ordering 
F® e : F® e : F® x = 2.40 : 1.60 : 1.0. Initial energy spectra for v e (black dashed curves) and v x 
(light dashed curves) and after collective oscillations for v e (black continuous curves) and v x 
(light continuous curves) at r — 350 km. 


The first large-scale multi-angle simulations have been performed in 2006 in Ref. [5], 
adopting the neutrino bulb model. After that, different groups developed independent 
multi-angle simulations (see, e.g. [8,328,349,357-359]) that extended the seminal findings 
of [5] and explored the dependence of the flavor evolution on the initial SN neutrino 
fluxes. Surprises and unexpected results were found with a strong dependence on many 
details (e.g. neutrino flavor asymmetries [328,347], angular distributions [360,361], three- 
flavor effects [330,331,345,346,358]). Despite the numerous studies on the subject, at 
the moment a complete picture of the self-induced flavor conversions under multi-angle 
effects in SNe is still missing. Here we show a few examples of the possible behavior due 
to self-induced flavor conversions in a 3v scenario in the presence of multi-angle effects 
during different post-bounce phases: 

- F® & F® x F® e . This is the typical flux ordering expected during the neutron- 
ization phase (see Fig. 9), where the v e flux is strongly enhanced with respect to 
u Xl while the v e flux is strongly suppressed. In this case, bipolar flavor conver¬ 
sions, proceeding through pair productions of v e v e v x v x , are not possible [7]. 
Therefore, only synchronized oscillations occur, with no relevant effect of flavor 
conversion, since the in-medium mixing angle is small. Multi-angle effects are also 
negligible ( 13 ). 

- F® e F® e F ® x . This flux hierarchy is typically expected during the accre¬ 
tion phase (see Fig. 9), where the first part of the hierarchy is caused by the 
deleptonization of the collapsed core, and the second is caused by the absence of 
charged-current interactions for neutrino species other than v e and v e . Figures 26 


( 13 ) A different scenario could be encountered in the case of low-mass SNe with an oxygen-neon- 
magnesium core. For these SN progenitors, the matter density profile can be very steep. The 
usual MSW matter effect occurs within the region of high neutrino densities close to the neutrino 
sphere. Therefore, self-induced flavor conversions will be possible at low-radii [354,362]. 
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E (MeV) E (MeV) 

Fig. 26. - Case with F^ e : F° e : F® x = 2.40 : 1.60 : 1.0. Three-flavor evolution in the single-angle 
case and in inverted mass hierarchy for neutrinos (left panels) and antineutrinos (right panels). 
Upper panels: Initial energy spectra for v e (long-dashed curve) and v x , y (short-dashed curve) 
and for v e after collective oscillations (solid curve). Lower panels: Probabilities P ee (solid 
red curve), P ey (dashed blue curve), P ex (dotted black curve). (Reprinted figure from [358]; 
copyright (2011) by the American Physical Society.) 
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Fig. 27. - The same as Fig. 26 but in the multi-angle case. (Reprinted figure from [358]; copyright 
(2011) by the American Physical Society.) 
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and 27, as examples of this configuration, show the three-flavor multi-angle flavor 
evolution in IH for a flux ordering F® e : F® e : F^ x = 2.40 : 1.60 : 1.0 (already used 
in Figs. 23 and 25). In particular, Fig. 26 refers to the single-angle evolution, while 
Fig. 27 is for the multi-angle case. The initial v e fluxes (dashed curves) and the 
final ones for v e and v x , y are represented in the upper panels. Note that v x<y are 
linear combinations of the v^ r fluxes, defined as v x ^ y = cos # 23 ^ Tsin0 2 3*V [341], 
with 023 — tt/ 4 [see Eq. (20)]. The conversion probabilities P ee ,P ey ,P ex are shown 
in the lower panels. As already discussed, one can see that the final v e flux is almost 
completely swapped with respect to the initial one, while the final v e flux presents a 
peculiar spectral split at E sp u t — 10 MeV, being swapped to v y at higher energies. 
Flavor conversions occur in the 2v (e — y) sub-system associated with Am 2 and #13 
(as shown in Fig. 25). From the comparison with the multi-angle case, we see that 
these effects play a sub-leading role, being suppressed by the large flavor hierarchy 
of the accretion phase [328]. No self-induced flavor conversion occurs in NH, as 
expected. 

The results presented in Figs. 26 and 27 have been obtained adopting an effective 
small mixing angle to simulate the matter effects; other multi-angle studies [333- 
339], conducted under a simplified setup, but including the matter background 
have shown as during the accretion phase the dense matter would dominate over 
the neutrino density, generally suppressing the self-induced flavor conversions (see 
Sec. 4*3). Therefore, it is not clear wether this type of flavor evolution is realized 
in SNe, since probably the flux hierarchy is not large during the cooling phase. 

- F® > F® > F® . This spectral ordering with a less pronounced flavor hierarchy 
among the different species is possible during the cooling phase (see Fig. 9). For 
definitiveness we take F^ e : i 7 ^ : F® x = 0.85 : 0.75 : 1.0. Figures 28-29 are in 
the same format adopted in Fig. 26 and 27 for the accretion phase. In particular, 
Fig. 28 refers to the single-angle case, while Fig. 29 is for the multi-angle evolu¬ 
tion. Differently from the accretion phase, multiple spectral splits are present in 
both neutrino and antineutrino channels. Three-flavor effects are observable in the 
single-angle scheme, while get suppressed in the multi-angle case. Moreover, the 
spectral swaps and splits are less pronounced, due to some amount of multi-angle 
decoherence in the flavor conversions. In this regard, complete decoherence could 
occur further reducing the flavor asymmetry [358]. Finally, we mention that for 
the flux ordering of the cooling phase spectral splits and swaps would occur also in 
NH. 

In Table I we summarize the results on the role of multi-angle effects, 3v effects (asso¬ 
ciated with 5m 2 ) and spectral splits for different SN neutrino fluxes. Matter suppression 
effects (see Sec. 4‘3.1) should be mainly relevant during the accretion phase, see two 
dedicated cases. 

4*8. Multi-azimuthal-angle instability. - As discussed until now, numerical studies of 
self-induced flavor conversions have been performed within the bulb model. However, it 
has been recently questioned if removing some of the symmetries assumed in this model, 
this would trigger new instabilities in the flavor evolution. A crucial assumption of the 
bulb model is the spherical neutrino emission from the neutrinosphere, that reduces to 
a cylindrical symmetry along a given neutrino trajectory (see Fig. 21 ). Recently, by 
means of a stability analysis of the linearized neutrino EoMs [363], it has been pointed 
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Fig. 28. - Case with F° e : F§ e : Ff} x — 0.85 : 0.75 : 1.0. Three-flavor evolution in the single-angle 
case and in inverted mass hierarchy for neutrinos (left panels) and antineutrinos (right panels). 
Upper panels: Initial energy spectra for v e (long-dashed curve) and v x , y (short-dashed curve) 
and for v e after collective oscillations (solid curve). Lower panels: Probabilities P ee (solid 
red curve), P ey (dashed blue curve), P ex (dotted black curve). (Reprinted figure from [358]; 
copyright (2011) by the American Physical Society.) 
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Fig. 29. - The same as Fig. 28 but in the multi-angle case. (Reprinted figure from [358]; copyright 
(2011) by the American Physical Society.) 
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Table II. - Summary of multi-angle effects, effects and spectral spits for different SN neutrino 
fluxes, assuming sub-leading matter effects. The cases “accretion, \ r <C Hr” (“accretion, \ r 
Hr”) stand for absence (presence) of multi-angle matter suppression effects. 


Initial spectral pattern 

Multi-angle effects 

£ra 2 -effects 

Spectral splits 

F° e > F® x > F° e (neutronization) 

no 

no 

no 

F° e > F° e > Ff, x (accretion, A r <C /i r ) 

marginal 

absent 

robust 

F° e > F° e > F% x (accretion, A r > Hr) 

relevant 

absent 

no 

F° x > F° e > Fl (cooling) 

relevant 

present / absent 

smeared 

F° e ~ F° x ~ Fl (cooling) 

strong 

present 

washed-out 


out that removing the assumption of axial symmetry in the v propagation, a new multi- 
azimuthal-angle (MAA) instability could emerge in the flavor evolution of the dense SN 
neutrino gas. The occurrence of this instability has been then clarified with simple toy 
models [364,365]. The presence of MAA effects unavoidably implies the breaking of 
the spherical symmetry in the flavor evolution after the onset of the flavor conversions. 
This would lead to a challenging multi-dimensional problem involving partial differential 
equations. However, assuming that the variations of the global solution in the direction 
transverse to the radial one are small, one can still study the local solution along a 
given line of sight, without worrying about its global behavior. This approach, even if 
not completely self-consistent, allowed to obtain the first numerical solutions of the non¬ 
linear neutrino propagation equations in SNe, introducing the azimuthal angle as angular 
variable in addition to the usual zenith angle in the multi-angle kernel [see Eq. (26)]. 

Adopting this generalized bulb model, it has been shown [366,367] that the pattern of 
the spectral crossings (energies where F Ve = F Vx , and F„ e = F „ x ) is crucial in determining 
the impact of MAA effects on the flavor evolution. For neutrino spectra with a strong 
excess of v e over z/ e , as expected during the accretion phase, new flavor conversions occur 
in NH. This is represented in Fig. 30, where we have considered initial fluxes as the ones 
of Fig. 25. In particular the initial (dashed curves) and final fluxes (continuous curves) 
are shown for e (black curves) and x (light curves) flavors, for neutrinos (left panels) 
and antineutrinos (right panels). Upper panels show the NH case while lower panels 
refer to the IH one. One realizes that MAA effects produce flavor conversions in NH, 
otherwise absent in the azimuthal symmetric model. However, the dominance of matter 
terms during the accretion phase would typically suppress also the MAA effects [338,339]. 
For spectra with a moderate flavor hierarchy, as the ones expected during the cooling 
phase, the growth of the MAA instability should be inhibited. As a result, the oscillated 
spectra are rather close to the ones seen in the azimuthal symmetric case; see Fig. 29 for 
comparison (see also [367]). 
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Fig. 30. - Multi-azimuthal-angle flavor evolution for z/s (left panel) and z^’s (right panel) in NH 
(upper panels) and IH (lower panels) for fluxes with an initial ordering F® e : F® e : F® x = 2.40 : 
1.60 : 1.0. Energy spectra for v e (black dashed curves) and v x (light dashed curves) in absence 
of flavor oscillations and for v e (black continuous curves) and v x (light continuous curves) after 
collective oscillations. (Reprinted figure from [367]; copyright (2014) by the American Physical 
Society.) 


4’9. Spontanueous breaking of space-time symmetries. - The MAA instability has 
shown as self-interacting neutrinos can spontaneously break the symmetries of the ini¬ 
tial conditions. Such insight stimulated further investigations about the validity of the 
solution of the SN neutrino EoMs worked out within the bulb model (see Sec. 4‘2). 

The assumptions of azimuthal symmetry in neutrino propagation and quasi-stationary 
neutrino emission reduce the general partial differential equations [see Eq. (22)] into 
ordinary differential equations describing the stationary spatial evolution of the dense 
neutrino gas along the radial direction. These assumptions, although adopted since the 
earliest papers on the subject, are unjustified. Indeed, it was tacitly assumed that small 
deviations from them do not significantly perturb the flavor evolution. However, it has 
been recently realized that instabilities may grow once initial symmetries are relaxed, 
since self-interacting neutrinos can spontaneously break the translation symmetries in 
time [368-370] and space [371-375]. This implies that the characterization of the self- 
induced effects obtained within the spherically symmetric bulb model should be taken 
cum grano salis. 

A self-consistent solution of the flavor evolution would eventually lead to solve the 


71 














Supernova Neutrinos: Production, Oscillations and Detection 


complete space-time-dependent problem described by the following partial differential 
equations [317] 


(48) 




Qt, p,x 


^[^P,X) Qt, p,x] 5 


where D p x = Dp ac + and an analogous equation exists for antineutrinos [see 

Eq. (22)]. This represents a formidable seven-dimensional problem. First attempts of 
solution have been recently presented in [368,372,373] by Fourier transforming Eq. (48). 
For example, performing a spatial Fourier transform one finds 


(49) 


d 

m +Vp 


•k 


Qt, p,k 


i J d 3 x e * k x [ri PjX 


Qt, p,> 


which represents a tower of ordinary differential equations in t for the different Fourier 
modes Pt,p,k with wavenumber k, which are coupled through the interaction term. This 
technique sheds the basis to study this challenging problem in the SN case. 

In this context, it has been recently pointed out that even tiny space inhomogeneities 
could lead to new flavor instabilities [371,372], developing at small scales, even at large 
neutrino densities, where oscillations are otherwise expected to be suppressed due to syn¬ 
chronization [371]. However, as shown in [375], large neutrino densities in a supernova are 
typically accompanied by a large matter density, which produces the multi-angle matter 
effects (Sec. 4*3.1) that suppresses the low-radii small-scale instabilities. Furthermore 
in [370], it has been shown that removing the assumption of the stationarity of the flavor 
evolution, these results can dramatically change. Indeed, the steady solution is not stable 
and self-interacting neutrinos can generate a pulsating solution with a frequency that ef¬ 
fectively compensates the phase dispersion associated with the large matter term, lifting 
the suppression of the space instabilities at small-scales. If these space-time instabilities 
develop fully and cascade, that paves the way for flavor conversions at large neutrino and 
matter densities. The flavor-content of SN neutrino fluxes will be strongly varying with 
time, perhaps leading to flavor-averaging and have profound consequences for supernova 
explosions and nucleosynthesis. 

The above described very recent results add new challenges in the characterization 
of the self-induced SN neutrino flavor conversions. Remarkably, the final goal would be 
a realistic treatment of the self-induced flavor conversions within the SN environment, 
where large deviations from a spherical neutrino emission can be generated by hydro- 
dynamical instabilities such as SASI or LESA, especially during the accretion phase, as 
discussed in Sec. 2. 

Another assumption that has recently been questioned is that neutrinos are free- 
streaming after the neutrinosphere [376], as even a fraction of neutrinos that occasionally 
scatter outside of the neutrinosphere should produce a small “neutrino halo” [376,377]. 
On the basis of a stability analysis, it has been concluded that the neutrino halo should 
not drastically affect flavor oscillations during the accretion phase [378]; however a self- 
consistent numerical simulation is not available at the moment. The inclusion of the 
halo effect in the numerical simulations would change the nature of the flavor evolution, 
turning it into a boundary value problem instead of a non-stationary initial value one. 
All these open issues require further dedicated work to fully clarify the role of flavor 
instabilities in the interacting neutrino field. 


72 



Supernova Neutrinos: Production, Oscillations and Detection 


4*10. Mikheyev-Smirnov-Wolfenstein matter effect in the wake of the shock-wave. - 
Self-induced effects in SN neutrino oscillations, occurring in the deepest stellar regions 
would eventually die out at r > 0(1O 3 ) km. However, as neutrinos stream through the 
outer layers of the stellar envelope, they would feel ordinary matter effects. There is a 
wide literature on (SN) neutrino oscillations in matter to which we refer the interested 
reader (see, e.g., [304,341,379-382]). As the SN matter potential A r = ) 

declines, neutrinos would eventually encounter Mikheyev-Smirnov-Wolfenstein (MSW) 
resonances when [304,341] 

(50) A r = lvh,l , 

corresponding to the atmospheric Am 2 (P-resonance) and the solar 5m 2 (L-resonance) 
mass-squared differences, respectively. 

Typical regions where we expect H and L resonances have been shown in Fig. 22. We 
realize that the two resonant regions are rather separated due to the hierarchy Sm 2 <C 
Am 2 , so that one can typically factorize the dynamics studying the resonant effects in 
sub-sectors. When one of the conditions in Eq. (50) is fulfilled a resonant amplification of 
the flavor conversions is expected. From the Hamiltonian (in absence of self-interactions), 
one can build the level crossing diagrams for the two mass hierarchies (see Fig. 31) 
showing the neutrino propagation eigenstates from the regions at high-density to the 
vacuum, where A r = 0 [341]. These diagrams allow us to determine in which mass 
eigenstate will emerge a neutrino produced in a given interaction eigenstate. In the case 
of antineutrinos, the effective potential for v e is A r = — v // 2 G j pn e (r). Antineutrinos can 
be represented on the same level crossing diagram, as neutrinos traveling through matter 
with “effectively” negative n e (r). Given the resonance conditions in Eq. (50), the sign 
of the matter potential A r and of the mass-squared splittings, we realize that the H- 
resonance can be satisfied by neutrinos in NH (A r , Am 2 > 0) or by antineutrinos in IH 
(A r ,Am 2 < 0). Therefore, in principle, the neutrino burst is sensitive to the neutrino 
mass hierarchy thanks to the matter effects, associated with the P-resonance. Conversely, 
the L-resonance can be satisfied only by neutrinos in both the mass hierarchies (since 
Sm 2 > 0). 

For typical SN simulations, the matter density profile declines so slowly that the 
neutrino propagation is adiabatic , i.e., each mass eigenstate in Fig. 31 remains the same. 
However, this condition is violated at the (forward and reverse) shock-fronts where, 
due to the abrupt density change (see Fig. 22), strong non-adiabatic conversions occur. 
Neglecting self-induced flavor conversions, the neutrino flux arriving at Earth can be 
expressed in terms of energy-dependent v e survival probabilities P ee (E) at the shock- 
front [341]: 

(51) F Ue = P ee (E)F? e (E) + [1 - Pee(E)]Fu x (E) ■ 

An analogous expression exists for v e with survival probabilities P ee (E). In particular, 
considering for simplicity a static SN matter profile and a complete adiabatic propagation 
we have (P ee , Pee) = (0, cos 2 9 12 ) in normal mass hierarchy, and (P ee , P ee ) = (sin 2 # 12 , 0) 
in inverted mass hierarchy [341]. 

Corrections to a pure adiabatic neutrino propagation are expressed in terms of level 
crossing probabilities among the instantaneous eigenstates in matter at the resonance 
point [304]. As shown in Fig. 22 , the real SN density profile is non-monotonic and 
time-dependent, so that multiple resonances would occur along it. For #13 as large as 
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Fig. 31. - Three-flavor level diagram for neutrino propagation eigenmodes, relevant for neutrinos 
streaming from a SN core [341] for normal hierarchy (left) and inverted hierarchy (right). (Figure 
taken from [20] with permission.) 


recently measured, the JT-resonance is adiabatic, except at the shock-fronts, where the 
crossing probability Ph = Ph( Am 2 , # 13 ) between the eigenmodes z/ 3 ?m and 1 ^ 2 ,m (see 
Fig. 31) would be extremely non-adiabatic, giving Ph — 1. Then, as z/s propagate at 
larger radii, they would eventually encounter the L-resonance between the 1 / 2 , m and zq ;m 
states, associated with the ((bn 2 , # 12 ) sub-sector. The L-resonance intercepts the shock- 
front only at relatively late times (t p b > 10 s) and it is never strongly non-adiabatic. 
Therefore, sub-leading effects related a level crossing probability Pl 0 are typically 
neglected. 

The survival probability P ee of SN neutrinos at Earth (neglecting Earth matter cross¬ 
ing) is related to the crossing probability Ph at the shock-front [383]: 


(52) 


f sin 2 6 12 P h (u, NH), 

I cos 2 <? 12 (V, NH), 

| sin 2 612 (iA IH), 
[cos 2 9 12 P H (17, IH). 


It is clear as Ph can modulate the (otherwise constant) survival probability of v e in 
NH and of v e in IH, thus providing an important handle to solve the current hierarchy 
ambiguity. 

In order to compute Ph , one has to numerically integrate the flavor evolution equa¬ 
tions [Eq. (44)] in the basis of the instantaneous matter eigenstates. Analytical prescrip¬ 
tions have been presented in [383,384]. In Fig. 32, we show the crossing probability Ph 
as a function of energy E for the SN matter density profile of Fig. 22 at post-bounce 
time t p b = 5 s. We closely follow the approach of Ref. [383] to compute Ph (the in¬ 
terested reader is referred to this reference for further details). The crossing probability 
has a typical top-hat structure, jumping from Ph ~ 0 (adiabatic regime) to Ph ~ 1 
(extremely non-adiabatic regime) when the resonance condition is satisfied across the 
forward shock-front. A more complicated pattern would emerge when resonances occur 
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Fig. 32. - Neutrino crossing probability Ph as function of the neutrino energy E for the SN 
matter potential of Fig. 22 at post-bounce time t p b = 5 s. 


on both the forward and the reverse shock fronts [246,340]. It is expected that the 
transient violation of the adiabaticity condition, when neutrinos cross the shock-fronts, 
would emerge as an observable modulation of the neutrino signal. This signature could 
be particularly useful to follow in “real-time” the shock-wave propagation, as well as 
to probe the neutrino mass hierarchy. This opportunity has been widely discussed in 
literature [246,340,383,385-388]. 

A realistic characterization of matter effects, during neutrino propagation across the 
SN shock-wave, must also take into account stochastic density fluctuations, inhomo¬ 
geneities of various magnitudes as well as correlation lengths in the ejecta layer in the 
wake of the shock front. These fluctuations are a result of hydrodynamic instabilities 
between the proto-neutron star and the SN shock during the very early stages of the 
SN explosion. They lead to large-scale explosion asymmetries and turbulence in a dense 
shell of shock-accelerated ejecta, which subsequently also seed secondary instabilities in 
the outer shells of the exploding star (e.g., Refs. [135,389-394]). 

Neutrino flavor conversions in a stochastic matter background have been subject of 
intense investigations both in a general context [395-401] and specifically in relation to SN 
neutrinos [402-409]. It is expected that stochastic matter fluctuations of sufficiently large 
amplitude may suppress flavor conversions and lead to Ph — 1/2 when the suppression is 
strong [402]. Therefore, the spectral properties of the fluctuations are very important for 
understanding the neutrino signal emerging from a core-collapse SN. At the moment there 
is no unanimous consensus about the impact of matter fluctuations on the SN neutrino 
flavor conversions. A recent study based on two-dimensional SN simulations has shown 
a modest damping of the neutrino crossing probabilities [410]. However, further analysis 
would be mandatory when high-resolution three-dimensional simulations will become 
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Fig. 33. - Observable SN v e spectra in a 40-kton LAr TPC for benchmark fluxes as in Fig. 28 
in both the mass hierarchies. (Figure adapted from [356]; courtesy of B. Dasgupta.) 


available. 

4*11. Observable signatures of supernova neutrino flavor conversions. - After dis¬ 
cussing in detail how self-induced and matter effects would process the SN neutrino 
fluxes, we present some of the possible observable signatures of flavor conversions in SNe 
and their sensitivity to the neutrino mass hierarchy. For definitiveness, we will refer to a 
Galactic SN at a distance d = 10 kpc from the Earth. We address the interested reader 
to Sec. 3 for a detailed discussion on SN neutrino detection techniques. 

(a) Self-induced spectral splits. As discussed in this Section, in recent years the picture 
of SN neutrino oscillations, based only on the MSW matter effects, has undergone a 
change of paradigm by the insight that the refractive effects of neutrinos on themselves 
are crucial. Observationally, the most important consequence of the self-induced flavor 
conversions is a swap of the v e and v e spectrum with the non-electron species v x and 
v x in certain energy intervals, and the resultant spectral splits at the edges of these 
swap intervals. Some of the spectral splits could occur sufficiently close to the peak 
energies to produce significant distortions in the SN neutrino signal, observable in the 
large underground dectectors. As an example, we show in Fig. 33 the observable v e signal 
in a 40-kton LAr TPC, while Fig. 34 refers to the v e signal in a 20-kton liquid scintillator 
detector. Both the mass hierarchy cases are shown. The initial fluxes in these figures are 
the same as in Fig. 28 (during the cooling phase) where self-induced conversions are not 
matter-suppressed. The self-induced flavor evolution has been characterized in a single¬ 
angle and three-flavors scenario (as in Fig. 28). The MSW effect is calculated according 
to Eq. (51) and (52), where it is assumed complete adiabatic propagation ( Ph = 1) 
neglecting possible shock-wave effects (see [356] for details). 

Concerning the electron spectrum produced by v e in a LAr TPC (Fig. 33), the observ- 
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Fig. 34. - Observable SN v e spectra in a 20-kton liquid scintillator detector for benchmark fluxes 
of Fig. 28 in both the mass hierarchies. (Figure adapted from [356]; courtesy of B. Dasgupta.) 


able flux is mostly due to the initial F® and then it does not present any special spectral 
feature in both the mass hierarchies. Conversely, for the v e signal in a scintillator detec¬ 
tor (Fig. 34), the observable positron spectrum would be mostly due to F® e for E < E sp n t 
and F® x at higher energies (see also Fig. 28). This would produce a bimodal positron 
spectrum, with two peaks produced by the two initial antineutrino distributions. The 
result is similar in both the mass hierarchies. This example shows that spectral splits are 
potentially identifiable during the cooling phase, if the average energies and luminosities 
of non-electron fluxes are sufficiently large. 

(b) Neutrinos from the SN neutronization burst. The v e neutronization burst is a 
particularly interesting probe of flavor conversions, since it can be considered almost as 
a “standard candle” being independent of the mass progenitor and nuclear EoS [140]. As 
discussed in Sec. 4*7, self-induced effects are not operative on the SN v e neutronization 
burst because of the large excess of v e due to the core deleptonization. Therefore, the v e 
flux would be only affected by MSW effects. At the very early post-bounce times relevant 
for the prompt burst, MSW flavor conversions would occur along the static progenitor 
matter density profile. We expect that the observable u e flux at Earth would be [Eq. (51)] 

(53) F Ve = F° x (NH) , 

(54) F Ve = sin 2 0 12 F° + cos 2 6 12 F° x (IH) . 

The observation of the neutronization peak would indicate the inverted mass hierarchy, 
while we do not expect its detection in NH. 

Concerning the possibility to detect the v e neutronization burst, while Mton class 
WC detectors measure predominantly the v e flux using inverse beta decay, they are also 
sensitive to the subdominant v e channel via elastic scatterings on e~ (see Sec. 3 ). On 
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Fig. 35. - Neutronization events rate per time bin in a 560-kton WC detector for 27 M© and 
11 M© SN progenitors (see Sec. 2'4) at d = 10 kpc for both NH (continuous curve) and IH 
(dashed curve). 


the other hand, a large LAr TPC will make the cleanest identification of the prompt v e 
burst through its unique feature of measuring v e charged-current interactions, enabling 
to probe oscillation physics during the early stage of the SN explosion. We show the 
observable v e neutronization burst from 11 M 0 and a 27 M 0 SN progenitors with LS EoS 
(see Sec. 2‘4) in a 560-kton WC detector (Fig. 35) and in a 40-kton LAr TPC (Fig. 36); 
the neutronization burst peak is clearly visible in IH in both detectors. Moreover, from 
these two figures one realizes that the variation in the progenitor mass has a sub-leading 
impact in the features of the signal. 

The detection of the SN neutronization burst has also been proposed to constrain 
nonstandard scenarios, like Lorentz invariance violation [411], neutrino decay [412], os¬ 
cillations into light sterile neutrinos [413], neutrino-antineutrino oscillations mediated by 
a neutrino magnetic moment [414]. All these scenarios would lead to a suppression of 
the neutronization burst. 

(c) Rise time of the neutrino signal. The rise time of a Galactic iron-core SN V e 
light curve, observable in large underground detectors, can provide a diagnostic tool for 
the neutrino mass hierarchy. Due to the combination of matter suppression of collective 
effects at early post-bounce times (see Sec. 4*3) and the presence of the ordinary MSW 
effect in the outer layers of the SN, one expects that the observable v e flux at Earth 
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Fig. 36. - Neutronization events rate per time bin in a 40-kton LAr TPC for 27 M© and 11 M© 
SN progenitors (see Sec. 2"4) at d = 10 kpc for both NH (continuous curve) and IH (dashed 
curve). 


would be given by 

(55) F Pe = cos 2 6 1 2 F° + sin 2 0 12 j£ (NH) , 

(56) F Pe = Fl (IH) . 

It is clear that the F„ e flux at the Earth would basically reflect the original F® flux if IH 
occurs, or closely match the F® e flux in NH. As from state-of-the-art simulations [415,416], 
temporal profiles of the original V x and V e fluxes appear quite different in the accretion 
phase, and relatively model-independent. In particular, the V x signal rises faster than 
the V e one. Figure 37 shows the early post-bounce evolution of luminosities (left panels), 
and mean energies (right panels) for a set of eleven ID simulations with progenitor of 
different masses obtained by the Garching group. Due to this feature, it would be possible 
to distinguish the two mass hierarchies from the rise time on (9(100) ms scale. 

Figure 38 shows the expected overall signal rate R(t) in the IceCube detector for a 
Galactic SN (see Sec. 3 for details about the SN neutrino detection technique in IceCube). 
We refer to the case of a 15 M© SN progenitor also adopted for illustration in Fig. 37. 
The NH cases are shown with continuous curves, while the IH cases with the dashed ones. 
The right panel overlies the error size using 2 ms bins with typical error estimates from 
the photomultiplier background noise. The difference between the observed neutrino 
light curve in the NH and IH is evident. For the NH case, a relatively longer hump in 
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Fig. 37. - Early postbounce evolution of luminosities (left panels), and mean energies (right 
panels) for a set of eleven ID simulations with progenitors of different masses as obtained by 
the Garching group. Quantities for u e , z7 e , and v x are shown in the top, middle and bottom 
panel, respectively. The vertical line indicates the early timescale (100 ms). (Figure taken 
from [68,415].) 


the signal, associated with the accretion, is clearly visible. While in IH, the light curve 
has a sudden rise. 

An important issue is the dependence of this signature is from theoretical uncertain¬ 
ties, most notably the progenitor structure, EoS, and numerical schemes. The models 
tested in [239] provide quite a satisfactory test that the uncertainties associated with 
progenitor structure do not spoil the viability of the method. Concerning numerical 
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Fig. 38. - Supernova signal in IceCube assuming a distance of 10 kpc, based on a simulation for 
a 15 M© SN progenitor from the Garching group. In the right panel it is illustrated the error 
size using 2 ms bins with typical error estimates from the photomultiplier background noise. 
(Figure adapted from [239].) 


schemes, the same signature has been found also within the hydro dynamical simulations 
of the Basel/Darmstadt group [34] and independently with the simulations of the Tapir 
group [416]. On the other hand, a firm conclusion on the other uncertainties requires 
further studies. Therefore, given the potential importance of the rise time signature in 
shedding light on the unknown neutrino mass hierarchy, it would be mandatory to further 
explore the robustness of this feature with more accurate simulations. 

(d) Earth matter effect. It is widely known that if SN neutrinos should reach the 
detector from “below” [234], the Earth crossing would induce an energy-dependent mod¬ 
ulation in the neutrino survival probability [341,417]. The appearance of the Earth effect 
depends on the neutrino fluxes and on the mixing scenario. The accretion phase is par¬ 
ticularly promising to detect Earth crossing signatures because the absolute SN v flux is 
large and the flavor-dependent flux differences are also large. Due to the matter suppres¬ 
sion of the collective oscillations during the accretion phase, SN neutrino initial fluxes 
will be processed by the only MSW effect. In this situation the oscillated SN v fluxes 
at Earth (before Earth crossing) are given by Eq. (51) in terms of the energy-dependent 
v e survival probabilities P ee (E) (or P ee (E) for P e ). The Earth effect can be taken into 
account by just mapping the non-vanishing P ee —>> 1 — P 2 e and P ee —> 1 — P^e, where P^ e 
is the i >2 v e transition probability for neutrinos propagating through the Earth, and 
analogously for P^ e [417]. For large # 13 , the neutrino fluxes at Earth for NH are [341]: 

(57) F® = (1 - P 2e )Fl + P 2e Fl and F® = F^ , 
while for IH 

(58) F® = Fl and F® = (1 - P 2 e )F° + P 2e F^ , 
where F® indicates the neutrino fluxes after Earth crossing. 
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Fig. 39. - Observable signal E 2 F„ e with (continuous curve) and without (dotted curve) Earth 
crossing. 


The analytical expressions for P 2e and P 2e can be calculated for the approximate 
two-density model of the Earth [418]. When neutrinos traverse a distance L through the 
mantle of the Earth, these quantities assume a very simple form [341,417] 

(59) P 2e = sin 2 0i2 + sin 20^ sin(20^ - 20 12 ) sin 2 f 5 ™_ Sm ^f l) , 

\ 4P sin 2(J r {2 J 

(60) P 2e = sin 2 0i2 + sin 20^ sin(20^ - 20 i2 ) sin 2 L ^j > 

where 6 ^ and O 7 ^ are the effective values of #i 2 in the Earth matter for neutrinos and 
antineutrinos respectively [382]. The Earth crossing induces a peculiar oscillatory sig¬ 
nature in the neutrino energy spectrum. From Eqs. (57)—(58), it results that the Earth 
matter effect should be present for antineutrinos in NH and for neutrinos in IH, providing 
a potential tool to distinguish between these two cases. 

Since the typical event rate associated with an inverse beta decay process is oc 
P 2 P®(P)—the detection cross section being cr oc E 2 —we show in Fig. 39 the P 2 P® 
as representative of the observable Earth-modulated signal, where we have taken as ini¬ 
tial v fluxes the ones in Fig. 28. Earth matter effects could be measured in a single 
detector, if it has enough energy resolution and statistics to track the wiggles in the 
observed energy spectrum, induced by the neutrino oscillations in the Earth. In this 
context, a Fourier analysis of the SN neutrino signal has been proposed as a powerful 
tool to diagnose this modulation, identifying the peak associated with the Earth crossing 
in the power spectrum [419]. 
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The observability of the Earth matter effect largely depends on the neutrino aver¬ 
age energies and on the flavor-dependent differences between the primary spectra. In 
this regard, recent SN simulations indicate lower average energies than previously ex¬ 
pected [34,415] and a tendency towards the equalization of the neutrino fluxes of differ¬ 
ent flavors during the cooling phase (at t p b > 1 s). Motivated by these new inputs, an 
updated study on the observability of this effect at large next-generation underground 
detectors (i.e., 0.4-Mton WC, 50-kton scintillation and 100-kton LAr detectors) was per¬ 
formed in Ref. [420]. It has been found that the detection of the SN neutrino Earth 
matter effect could be more challenging than expected from previous studies. Remark¬ 
ably, it was argued that none of the proposed detectors shall be able to detect the Earth 
modulation for the neutrino signal of a typical Galactic SN at 10 kpc. This should be 
observable in a 100-kton LAr detector for a SN at few kpc, while all three detectors 
would clearly see the Earth signature for very close-by stars only (d rsj 0.2 kpc). 

The Earth effect could also produce a modification in the SN v e light curve measured 
by the IceCube neutrino telescope. Therefore, IceCube could detect the Earth effect by 
the relative difference in the temporal signals with a high-statistics Mton WC detector, 
if only one of the two detectors is shadowed [421]. However, using the neutrino fluxes as 
from recent SN simulations also IceCube, used as co-detector together with a Mton-class 
WC detector, seems not to be able to detect any sizable variation in the SN neutrino 
event rate for any Galactic SN. 

Synopsis of signatures of flavor conversions. To summarize, we present a synopsis 
of the discussed signatures of flavor conversions in SN neutrinos and their sensitivity 
to the neutrino mass hierarchy in Table III. We focus on early time neutrino signal 
(neutronization or accretion phase) where we assume that the self-induced effects are 
matter suppressed. Moreover, we take the MSW matter effects along a progenitor static 
matter density profile. 

Table III. - Different detectable signatures of SN neutrino oscillations at early times and sen¬ 
sitivity to the mass hierarchy. 


Mass Hierarchy 

Pee 

P 

1 ee 

v e burst 

v e rise time 

Earth effects 

NH 

0 

COS 2 #12 

absent 

long 

h e 

IH 

sin 2 #i 2 

0 

present 

short 

v e 


4T2. Outlook. - The dense SN core represents a unique laboratory to probe neu¬ 
trino flavor mixing in high-density conditions. Indeed, within a radius of a few hundred 
kilometers, the neutrino gas is so dense to become a “background to itself.” In these 
conditions, we have shown how the neutrino flavor evolution equations become highly 
non-linear, leading to surprising and counterintuitive phenomena, like neutrinos and an¬ 
tineutrinos of different energies collectively oscillating to some other flavor. Another 
surprising result is the partial or complete energy-dependent flavor exchange, namely 
spectral swapping or spectral splits, as neutrinos travel out from a dense central region 
to the low neutrino density of the outer layers. During the past few years, our under¬ 
standing of non-linear neutrino oscillations has seen substantial progress. Analytical and 
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numerical works by several groups have shed light on the basic picture of non-linear neu¬ 
trinos oscillations, although still many open issues remain to be clarified. In particular, 
it has been appreciated that self-induced flavor conversions are related to instabilities in 
the flavor space. While the self-induced flavor evolution is a non-linear phenomenon, the 
onset of these conversions can be examined through a standard stability analysis of the 
linearized EoMs [422] (see also [329]). Moreover, interesting conceptual studies have been 
also carried on to clarify some of the aspects related to the collective flavor dynamics in 
SNe [423-427]. 

The recent insight that interacting neutrino fields can spontaneously break symmetries 
inherent to the initial conditions in the flavor evolution together with hints of asymmetric 
neutrino emission from the the most recent SN simulations in 3D will require a critical re¬ 
investigation of the previous results obtained within the bulb model. Seminal studies in 
this directions have just started. Furthermore, it has been speculated that in the transi¬ 
tion region between the neutrino scattering-dominated regime (at high-densities) and the 
oscillation-dominated one (at lower densities) spin-flavor oscillations could take place, al¬ 
lowing for coherent transformations between neutrinos and antineutrinos [326,428,429]. 
A possible role of anomalous neutrino-antineutrino correlations in the development of 
the self-induced oscillations has also been discussed [430-432]. Remarkably, new physics 
scenarios could have a strong impact on the self-induced flavor conversions. In this 
regard, strong effects on collective oscillations have been found due to spin-flip transi¬ 
tions [433,434], triggered by a v magnetic moment 100 times larger than the Standard 
Model one, assuming typical magnetic fields in the SN envelope (B ^ 10 10 — 10 12 G). 
Moreover, the presence of non-standard flavor changing interactions among SN neutrinos 
and the background fermions [435] or among neutrinos themselves [436] would dramati¬ 
cally affect the self-induced flavor changes. The putative existence of extra light sterile 
neutrino families is also expected to modify the active neutrino fluxes and the above 
described picture of neutrino self-interactions [168,170,171]. 

An appropriate characterization of these effects would motivate additional work to 
investigate their impact on the observable neutrino signal [92,169,437], on the r-process 
nucleosynthesis in SNe [171,438,439], and on the neutrino energy transfer to the stalled 
shock wave [91,92,102]. 

Matter effects on SN neutrino oscillations can imprint peculiar signatures on the 
observable neutrino signal with a strong sensitivity to the mass hierarchy. Neutrino 
flavor conversions on the wake of the shock-wave propagation represent an intriguing 
possibility to follow the SN dynamics in real time through the neutrino signal. An 
accurate prediction of the observable oscillation signatures imprinted on the SN neutrino 
signal is extremely timely now, in relation to the intense experimental activity for the 
realization of large-volume detectors for low energy neutrino astronomy. 
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5. — Diffuse Supernova Neutrino Background 

Author: I. Tamborra 


While the occurrence rate of a Galactic SN is only < 3 per century [440,441] and the 
probability of detecting neutrinos from the next Galactic explosion will be challenged 
by the location of the SN with respect to the detectors on Earth, on average one SN 
is exploding every second somewhere in the Universe. The Diffuse Supernova Neutrino 
Background (DSNB) is the total flux of neutrinos and antineutrinos with MeV energy 
emitted by all SNe in our Universe. It is expected to be isotropic and stationary and it 
will provide us with a unique portrait of the SN population. Although not yet detected, 
the DSNB discovery prospects are excellent. The DSNB detection will be crucial to 
test our current understanding of the SN dynamics and SN redshift distribution. In 
this section, we report on the present status of the DSNB searches and its theoretical 
uncertainties. 

5T. Main ingredients and DSNB upper limits. - The DSNB flux depends on the SN 
distribution according to their progenitor mass M and redshift £ and on the neutrino 
energy spectra. For each neutrino flavor u a (a = e, e, u or r) with observed energy E, 
the DSNB is defined as [442]: 


(61) 




.(/-') „ 
rid 


, /“^rriax pZ 

- dM 

0 JMq J 0 


dz 


Psn(z, M)F U(x (E(l + 2 ?),M) 
V^m(1 + z) 3 + Da 


where c is the speed of light, Hq = 70 km s -1 Mpc -1 the Hubble constant, Dm = 0.3 
and Da = 0.7 are the matter and dark energy fractions of the cosmic energy density. 
F V(x (E, M ) is the oscillated u a energy flux for a SN progenitor with mass M introduced 
in Sec. 4*3. The comoving SN rate (SNR) is labelled as psN and it is expected to be 
larger between z = 0 and z max ~ 5 (see Sec. 5’1.1). The redshift correction of the 
energy in Eq. (61) is due to the fact that each neutrino emitted from a SN at redshift 
z will have energy E' = E(1 + z) being E the neutrino energy observed on Earth. 
As a consequence, neutrinos coming from high- 2 : are pushed towards lower energies; 
therefore the DSNB flux is dominated by the z < 1 contribution [442]. As it will be 
discussed in the next section, the SN mass distribution is such that the least massive 
stars (with M ~ Mq — 8 M©, being Mq the minimum progenitor mass necessary to have 
a core-collapse SN) give the larger contribution to the DSNB. The DSNB is then weakly 
dependent on M max ~ 125 M© (tentative upper limit for the occurrence of core-collapse 
SNe [442]) and 2 : max . 

Since first studies on the DSNB and on its possible detection [443-447], our knowl¬ 
edge on the topic has increased considerably (see, e.g., Refs. [448, 449] for dedicated 
review papers). From the experimental standpoint, a milestone was reached when 
Super-Kamiokande placed an upper limit on the v e component of the flux [450] ~ 

0.1-1. cm -2 s -1 above 19.3 MeV threshold), excluding some theoretical models [25]; 
while the upper limits on the other neutrino flavors above 19.3 MeV are less strong: 
(,K e < 73.30-154 cm _2 s _1 , <$> Vvk +v T < (1-1.4) x 10 3 cm -2 s -1 and 0p M +p T < (1.3- 
1.8) x 10 3 cm -2 s -1 [451]. The current most stringent upper limits have been placed by the 
Super-Kamiokande experiment on the v e component of the flux above 17.3 MeV thresh¬ 
old, (f)v e < 2.8-3.0 cm -2 s -1 at 90% C.L. [22,452,453]. More recently, Super-Kamiokande 
released a new upper limit by performing an analysis with a neutron-tagging technique: 
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4V e < 5-30. cm - 2 s - 1 MeV -1 for neutrino energies between 13.3 MeV and 17.3 MeV [454]. 
As we will discuss later, the most promising energy window relevant for the DSNB de¬ 
tection against background is 11 MeV < E < 40 MeV [28,452,455,456]. 

5*1.1. Cosmic supernova rate. One of the key ingredients to compute the DSNB flux 
is the distribution of core-collapse SNe (i.e., Type II SNe and the subdominant Type Ib/c 
SNe) with the redshift 2 and progenitor mass M. The SNR, psn(z,M), is proportional 
to the Star Formation Rate (SFR), Psf(z ), through the initial mass function (IMF), 
'ip(M): 


(62) 


Psn(z,M) 


'ip(M) 


Psf{z) ; 


therefore normalization uncertainties on the IMF do not affect Psn{z , M). The distribu¬ 
tion of stellar masses is assumed to follow a universal IMF, as the conventional Salpeter 
scaling law: !p(M) oc M -2,35 for stellar masses M between 0.1 M 0 and 100 M 0 [457], 
although even other IMF scaling laws have been suggested, such as an intermediate one 
[458] or a shallow one [459]. 

Recent studies [460,461] suggest that the SFR can be fitted by a piecewise function 
of the redshift, for example [462]: 


(63) 


Psf{z) oc 


(1 + z) 5 z< 1 
(1 + z) a 1 < z < 4.5 , 
(1 + z)'* 4.5 < 2 


with S = 3.28, a = —0.26, 7 = —7.8, normalized such that /^ max dM psn{ 0,M) = 
1.5 x 10 -4 Mpc - 3 yr -1 . Note that the DSNB dependence on psf{z > 1) is weak as the 
DSNB is mainly determined by the contribution coming from z < 1. The SNR defined 
as in Eq. (63) is a growing function of the redshift implying that SNe were most frequent 
in the past (see also Fig. 40). Given the definition of the SNR (via the IMF), the DSNB 
is dominated by the low-mass threshold Mq, although this quantity is difficult to predict 
accurately from theory since stellar properties change rapidly between 6-10 M 0 . On the 
other hand, the upper limit mass M max is less important because of the fast decline of 
the IMF with M. 

In order to estimate the DSNB, one could rely on direct measurements of the SNR, 
instead of using Eq. (62). However, while the shape of the SNR is basically well known, 
its measurements are governed by normalization errors [56,460]. Reference [56] pointed 
out that the normalization from direct SN observations was lower than that from SFR 
data by a factor ~ 2, this discrepancy is known as the “supernova-rate problem.” The 
reasons for such a mismatch may be manifold. It might be that a large number of SNe 
are actually dim (either due to dust obscuration or being intrinsically weak), or that 
systematic changes are needed in our understanding of stars and SN formation. 

Measurements of the SNR have greatly improved in the last few years [463-467]. 
Figure 40 shows the SNR from the most recent cosmic measurements in comparison with 
the one predicted from the SFR assuming that all SNe with M > 8 M 0 yield optical 
SNe [56,460,468]. It is clear as the SNR problem has been cured at low z where we 
now have a better understanding of the dust extinction in the host galaxies and of the 
SN luminosity function, but the SNR is still affected by large uncertainties at high z. It 
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Fig. 40. - Cosmic supernova rate as a function of the redshift. The SFRs from [56, 468] are 
shown as continue lines after extrapolating the SNR from the SFR as from Eq. (62), assuming 
that progenitors with mass 8 M© < M < 50 M© produce optical SNe. Measured SNRs from 
recent work are also shown. (Reprinted figure with permission from [466]; copyright (2012) by 
the American Astronomical Society.) 


might also be that the fraction of “invisible” SNe (i.e., not optically visible SNe) increases 
with the redshift [465,469]. 

The SNR estimation will hopefully further improve within the next few years. In 
fact SNR measurements will come from synoptic surveys which will scan the full sky 
and probe the SN population with high sensitivity (see, e.g., discussion in Ref. [470] and 
references therein). Such surveys are expected to pin down the error on the normalization 
of the SNR and, as a consequence, on the DSNB normalization. On the other hand, note 
that the forthcoming detection of the DSNB could provide independent constraints on 
the SFR [471]. 

5'1.2. Time and progenitor dependence of the neutrino energy fluxes. The dependence 
of the DSNB on neutrino oscillation physics and SN progenitors has been discussed 
in Ref. [473]. As shown in Sec. 2 and Fig. 37, the neutrino luminosities and spectra 
depend on the compactness of the progenitor star [472]. Moreover, the luminosities of 
the different flavors are almost equal during the cooling phase, while L Ve ,L Ve > L v T 
during the accretion phase. The mean energies are {E 9ej1/ ) > (. E Ve ) during the cooling 

phase (see Fig. 9 and Sec. 2 '4). 

Concerning the neutrino oscillation scenario (see Sec. 4), while the MSW effect is well 
understood analytically and it occurs far away from the neutrinosphere, neutrino self¬ 
interactions crucially depend on the initial hierarchy among neutrinos of different flavors, 
number of energy crossings and on the neutrino mass hierarchy. During the accretion 
phase (t p b < 1 s), complete or partial multi-angle matter suppression of neutrino self¬ 
interactions occurs because of the high matter potential [333,334,336,337]. While during 
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the cooling phase, multiple spectral splits might appear according to the neutrino mass 
hierarchy and the number of crossings in the non-oscillated spectra. However, the shape 
of the spectral splits is smeared and their size reduced due to the similarity of the non- 
oscillated flavor energy spectra and to v-v multi-angle effects, partial neutrino flavor 
conversions occur and the flavor transition become gentler as time increases [473]. 

As discussed in Sec. 4, it is possible to factorize the effects of v-v interactions and the 
MSW resonances for sake of simplicity, as they occur in well separated spatial regions. 
Denoting with Ff^ the fluxes after the collective oscillations, the unchanged will 
undergo the traditional MSW conversions after v-v interactions. The fluxes ( F v& and 
Ffy e ) reaching the Earth after both the collective and MSW oscillations for NH and IH 
are" [330,474-476]: 


(64) 

(65) 

( 66 ) 
(67) 


F™ = sin 2 0 12 [1 - P c (F£ a ,F$ a ,E)](FS' - + F° Uy , 

= cos 2 e 12 P c {F^,F^,E)(Fl - F° Vy ) + F° Uy , 

F™ = sin 2 9 12 P c (FZ e ,F? e ,E)(F2 e - F° Uy ) + F° y , 

F™ = cos 2 0 12 [1 - Pc{Ft a ,F$ a ,E)]{F* a - ) + F° . 


In Eqs. (64)-(67) we have considered multi-angle v-v interactions driven by Am 2 and 
$13 between v e and v y , while the other flavor ( v x ) does not evolve, the effects on self- 
induced flavor conversions induced by the third flavor being a negligible correction for 
our purpose [10,344-346]. 


5*2. Dependence of the DSNB from flavor oscillations and SN progenitor properties. 
- Adopting the SNR as from the SFR [Eq. (62)] and time-(and progenitor-) dependent 
oscillation physics described as in Sec.5'1.2, the dependence of the DSNB on the neutrino 
oscillation physics and on the SN progenitors has been studied in Ref. [473], adopting 
long-term SN simulations for three progenitors with masses M = 18,10.8,8.8 M 0 [34] 
(with Shen EoS [150]). The features of the neutrino signal of these three progenitors have 
been extended to the whole estimated SN mass range ([Mq, M max ] = [8 M 0 ,125 M 0 ]) and 
time-dependent oscillated neutrino spectra have been included in the DSNB computation 
implementing MSW and multi-angle v-v interactions as a function of the post-bounce 
time. A quantitative estimation of the dependence of the time-integrated F Ua (E) on 
oscillation effects and on the stellar population is provided in Table IV [473]. Notice 
that, as we will discuss later, the upper limit of the SN mass progenitor as well as the 
neutrino fluxes might change including failed SNe [477]. 

Figure 41 shows the DSNB, E 2 x [with Q V(x as in Eq. (61)], as a function of the 
energy for v e ( v e ) on the top (bottom). In the plots on the right, the region of interest 
for Super-Kamiokande detection (E > 17.3 MeV) is zoomed. A variation of 10-20% due 
to the neutrino mass hierarchy is visible at E ~ 20 MeV. In order to investigate the 
DSNB dependence from the SN cooling and progenitor mass, a case using t p b = 0.5 s as 
representative of all post-bounce times and the 10.8 M 0 SN progenitor as representative 
of the whole SN population is shown (thin lines) in comparison to the one computed 
including all SN mass models (thick lines). 

No signature due to v-v interactions is visible in Fig. 41, since neutrino self-interactions 
are suppressed during the accretion phase and they occur at different energies for each 
post-bounce time during the cooling [473]. Note that Ref. [34] does not include nucleon 
recoil effects responsible for reducing the differences among the mean energies of differ¬ 
ent flavors during the cooling phase [33,64,478]. Therefore, the adopted SN inputs are 
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Table IV. - Mean energies for the time-integrated neutrino flux (E) in normal and inverted 
hierarchy. The listed cases are: (a) F„ a (E) obtained considering t p b = 0.5 s as representative 
post-bounce time and one SN progenitor (10.8 Mq) as representative of the whole stellar popu¬ 
lation, no oscillations; (b) F Ua (E) obtained including time-dependent fluxes for the 10.8 M 0 SN 
progenitor, no oscillations; (c) F Ua (E) with time-dependent fluxes for the 10.8 M 0 SN progeni¬ 
tor, MSW effects included; (d) F U(X {E) with time-dependent fluxes for the 10.8 M 0 progenitor, 
MSW and neutrino self-interactions included; (e) F U(X {E) obtained including the whole stellar 
population, time-dependent fluxes, MSW and neutrino self-interactions. (Adapted from [473]J 
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Fig. 41. - Left panels: DSNB as a function of the neutrino energy. Right panels: Zoom of 
the DSNB flux in the region of Super-Kamiokande detection (E > 17.3 MeV). The v e (p e ) 
flux is plotted in black (red). The solid (dashed) line represents the DSNB in IH (NH). The 
thick lines represent the DSNB including the three SN progenitors and the thin line shows 
the DSNB obtained considering the 10.8 M 0 SN progenitor as representative of the whole 
stellar population. (Reprinted figure from [473]; copyright (2012) by the Institute Of Physics 
Publishing.) 


useful to estimate the largest oscillation effects that could realistically affect the DSNB. 
Note that any flavor conversion occurring during the accretion phase will have a larger 
impact on the total DSNB, the luminosities during the accretion phase being larger than 
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the cooling ones. While, during the cooling phase, the fluxes are similar and any flavor 
conversion scarcely affect the DSNB. 

Summarizing the impact of the different ingredients on the DSNB: MSW effects are the 
largest source of variation of the DSNB with respect to the case without oscillations (50- 
60%), the mass hierarchy induces variations of ~ 20%; neutrino-neutrino interactions, 
time-dependence of the neutrino energy spectra over ~ 10 s and stellar population are 
responsible for a smaller variations 5-10%) [473]. Astrophysical uncertainties in the 
determination of the DSNB are instead larger: Besides the SN rate problem discussed 
above, the SNR is affected by a normalization error of ~ 25%, although such error 
is expected to be dramatically reduced within the next decade. Another potentially 
small source of error on the DSNB is the nuclear EoS [479]. Moreover, Refs. [480,481] 
pointed out as SN models with longer shock revival time produce a resultant higher 
DSNB event rate, since the larger the gravitational energy of the accreted matter, the 
larger is the resultant emitted neutrino energy. In principle, the shock revival time 
might be progenitor dependent and therefore introduce a non-negligible variation of the 
high-energy tail of the DSNB. 

Although theoretical studies have pointed out that the DSNB is sensitive to the SN 
progenitor mass, EoS, variations of the shock revival time, neutrino mass hierarchy, and 
oscillation physics, we stress that the DSNB carries integrated population information. 
In view of our current knowledge of the flavor oscillation physics and stellar dynamics, 
it will therefore be extremely difficult to extract information on each of the mentioned 
aspects from a forthcoming DSNB detection. In this sense, the high statistics of neutrino 
events potentially provided by the next galactic explosion will be more useful to constrain 
the SN engine and to learn about its dynamics. As we will discuss in the next section, the 
DSNB might also be affected by the presence of invisible SNe (probably also responsible 
for the mismatch between the measured SNR and the one extracted from the SFR) and 
the DSNB detection will be an extremely powerful instrument to learn about the whole 
SN population. 

5*3. Contribution from invisible supernovae. - Besides stellar core-collapse events that 
trigger SN explosions, a significant fraction (some 10-30%, possibly even more; e.g., [56- 
59,180,181,482-484]) could lead to black holes (BH) “directly,” i.e., after a relatively short 
(fractions of a second to seconds) accretion phase of a transiently existing hot neutron 
star. Such “failed SNe” are expected to produce only faint electromagnetic emission by 
stripping their loosely bound hydrogen envelopes [59,485,486]. The dim emission and 
red color make them hard to detect at extragalactic distances. Their neutrino emission, 
however, could account for an important contribution to the DSNB. In Sec. 2*5 neutrino 
signals of such events were presented as predicted by recent hydrodynamic simulations. 
The Super-Kamiokande Collaboration [26, 27] has provided constraints on the neutrino 
flux from failed SNe, setting a limit about a factor of four higher than predicted in 
Ref. [477]. 

Black hole formation may also occur with considerably longer delay, if the SN mecha¬ 
nism does not release enough energy to unbind the whole star. In this case a fair fraction 
of the stellar matter can fall back to the neutron star (on time scales of minutes to hours) 
and can push the neutron star beyond the mass limit for the black-hole formation. Such 
“fallback SNe” would radiate SN-like, though possibly sub-luminous, electromagnetic 
displays [487]. Their neutrino signals, however, are initially indistinguishable from those 
of forming neutron stars in ordinary core-collapse SNe. Whether the potentially massive 
fallback of matter to the neutron star produces observable neutrino emission is presently 
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not clear because of the lack of detailed studies, which are difficult due to the long time 
scales involved and the likely importance of asymmetries and angular momentum (but 
see Ref. [488]). Since the rate of fallback SNe in the present-day Universe is predicted to 
be low [181], we will not consider them in our discussion of the DSNB. 

Taking into account the current uncertainties, we will focus on the contributions from 
core-collapse and failed SNe to the DSNB in the following. We will see that the neutrino 
emission of failed SNe may indeed not be negligible but instead could make the more 
luminous component of the observable DSNB. 

Neglecting the progenitor mass dependence of the neutrino fluxes as well as of the SN 
rate, Eq. (61) can be generalized as 

(68) $„ a (E) = $™(E) + $™(E) = 

= ^ f Z - dz /W^/bh^T*^ 1 + *)) + /ns F* s (E(1 + z))] 

Hq Jo (1 + W “I - 

with /ns (/bh) the fraction of core-collapse (failed) SN, such that /ns + /bh = 1, where 
we have subdivided the whole SN family in two major categories, neglecting the fallback 
SN class (that could easily fall in one or the other case according to the progenitor mass) 
and assumed that both these families follow the SFR redshift dependence. 

The expected diffuse flux has been discussed in Ref. [477], assuming typical neutrino 
spectral parameters for the NS case as in Ref. [31] (i.e., Lp e = Lp x = 5 x 10 52 erg/s, 
(Ep e ) = 15 MeV, (Ep x ) = 18 MeV, otp e = 3.5, a Pe = 2.5) and modeling the BH case as in 
Ref. [176] for a typical failed SN progenitor of 40 M 0 with stiffer Shen (S) EoS [150] (i.e., 
Lp e = 12.8 x 10 52 ergs, Lp x = 4.9 x 10 52 erg/s, (Ep e ) = 23.6 MeV, (Ep x ) = 24.1 MeV) and 
softer Lattimer-Swesty (LS) EoS [90] (i.e., Lp e = 4.5 x 10 52 ergs, Lp x = 2.2 x 10 52 erg/s, 
(Ep e ) = 20.4 MeV, (Ep x ) = 22.2 MeV). 

The uncertainty of the fraction of failed SNe has been parametrized with /ns £ 
[0.78 — 0.91], corresponding to an upper limit for neutron-star-forming collapses varying 
in the interval 24-50 M 0 . Moreover, since the role of neutrino oscillations (especially of 
neutrino self-induced flavor conversions) has not been explored in this context, only MSW 
effects have been included: Fp e (E) = P ee F® e (E) + (1 — P ee )F®- x (E), with P ee = 0,0.68 
(see Tab. Ill, with P ee = 0.68 coming from the dependence of the survival probability 
from #12). The corresponding DSNB is shown in Fig. 42. 

Assuming flavor conversions play a non-negligible role (i.e., P ee = 0.68), the largest 
contribution from failed SNe comes from the stiffer EoS (S EoS case). In the latter case, 
the 4 >bh enhances the DSNB for E > 20 MeV, improving the DSNB detection chances 
(compare dashed with continue lines in Fig. 42) [477]. More recently, similar conclusions 
have been reached by adopting the progenitor dependence as from hydrodynamical sim¬ 
ulations of a range of SN progenitors and redshift dependence of the BH-formation rate 
based on the metallicity evolution of galaxies [481]. 

Although the largest contribution to the DSNB comes from z < 1, the flux from 
failed SNe from higher redshifts might still be substantial [465,469,489]. In this sense, 
the DSNB detection could directly allow to investigate the black-hole forming collapse, 
constrain its energetics and cosmological rate even beyond z ~ 1. 

If the neutrino emission is larger for failed SNe than for ordinary SNe as expected, 
then the signal increase in the detectors will be significantly larger in the high-energy tail 
of the DSNB spectrum [28,174-176,477,481,490]. This is promising in terms of neutrino 
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Fig. 42. - Diffuse flux of v e from core-collapse SNe (failed SNe) as a function of the energy plotted 
with dashed (solid) curves for S and LS EoS and survival probability P = 0.68 (including MSW 
effects, left panels) and P ee = 0 (no oscillations, right panels). The fraction of failed SN is 
/bh = 22% (thick curves) and /bh = 9% (thin curves). (Reprinted figure with permission 
from [477]; copyright (2009) by the American Physical Society.) 


detection chances and it would be a unique opportunity to study such SN progenitors, 
being otherwise optically invisible [491,492]. 

As discussed in Ref. [470], we can already constrain the visible and invisible SN rate 
at 2 = 0 on the basis of current data under the assumption that the shape of the SNR 
is known. It emerges that the allowed region for invisible SNe is indeed non-zero, but it 
cannot be arbitrarily large. Future observations will be extremely helpful on restricting 
the allowed region for visible SNe. 

Assuming that in the next few years, upcoming SN surveys will be able to pin down 
to 5% the uncertainty on the SNR, the variation of the one-year detection rate in Super- 
Kamiokande for different fractions of /bh is discussed in Ref. [470] by approximating the 
neutrino energy spectrum with a Fermi-Dirac distribution with zero chemical potential 
and total energy emitted for each flavor to L Ua = 0.5 x 10 53 erg, (E a ) = 12.6 MeV 
for the failed SNe, while L Va = 0.9 x 10 53 erg, (E Va ) = 23.6 MeV for core-collapse 
SNe. Figure 43 shows the expected one-year detection rate in Super-Kamiokande for 
/bh = /invis = 0,10,40%. It is evident as, increasing the fraction of invisible SNe, a 
higher detection rate is expected especially in the high energy tail of the energy spectrum. 
Note as, any indirect experimental constraint of this kind on the the NS/BH formation 
ratio will be also useful to constrain the still disputed mechanism(s) by which stars 
explode as SNe. 

We assumed that “failed” SNe are a fraction of the total SNR through this Section. 
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Fig. 43. - One-year neutrino detection rate in Super-Kamiokande as a function of the observed 
neutrino energy for three different fractions of invisible SNe (/bh = /invis = 0,10,40%). The 
band thickness of the curves is given by 5% of the uncertainty on the SNR as expected from 
upcoming SN surveys. (Reprinted figure with permission from [470]; copyright (2010) by the 
American Physical Society.) 


However, a first attempt to model the rate of failed SNe as a function of the redshift has 
been done in Ref. [469], based on gamma-ray burst observations which could originate 
from core collapses yielding rapidly rotating BHs. Reference [469] seems indeed to find 
a much higher fraction of failed SNe at z > 1 than locally. See also Refs. [181,481] for 
dedicated discussions. 

5‘4. Detection perspectives. - The DSNB detection is one of the main goals of neutrino 
astrophysics. However, the DSNB detection is strongly affected by backgrounds, which 
reduce the detectable energy window to 11 MeV < E < 40 MeV [28,452]. Figure 44 
shows the expected DSNB, in the presence of standard core-collapse (blue line) and failed 
(brown line) SNe and their sum (black line), as discussed in Sec. 5‘3. The expected 
atmospheric, reactor and solar neutrino backgrounds are also plotted. 

Water Cherenkov detectors, such as Super-Kamiokande, being mostly sensitive to F e , 
have atmospheric and reactor neutrinos as major backgrounds. The atmospheric neu¬ 
trinos fall within the same DSNB energy range and have an isotropic distribution in 
space, therefore it is very difficult to discriminate the DSNB signal from the background. 
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Fig. 44. - Expected diffuse flux of v e from core-collapse SNe in blue and failed SNe (S EoS, 
Pee — 0.68 and /ns = 0.78) in brown. Their sum is in black (see Sec. 5'3). The atmospheric 
and reactor fluxes are shown for the Kamioka (solid, gray) and Homestake (dashed, red) sites. 
Since the v e and v e atmospheric fluxes are very similar, only one is plotted. (Reprinted figure 
with permission from [489]; copyright (2012) by the American Physical Society.) 


However, the DSNB is larger than the atmospheric background around 30-40 MeV (see 
Fig. 44), thus restricting the experimental detection window to this range [246]. Reac¬ 
tor neutrino events could instead be distinguishable from their direction [493], but no 
dedicated study concerning the DSNB exists yet. Liquid argon detectors are mostly 
sensitive to z/ e , and therefore solar and atmospheric neutrino fluxes are the main back¬ 
grounds [494]. Solar neutrinos prevent studying the DSNB below 18 MeV for all the 
above-mentioned detectors (see Fig. 44); however they can be subtracted very effectively 
since their direction and flux can be well determined( 14 ). Clearly the detection-energy 
windows would be larger for an enhanced DSNB (i.e., including the contribution from 
failed SNe), and therefore a larger /bh would be advantageous for signal to background 
discrimination. 

An enhanced signal discrimination over the background for the DSNB detection in 
Super-Kamiokande will be reached by dissolving gadolinium in water [207] (see Sec. 3); 
this option is currently being tested through the Evaluating Gadolinium’s Action on 
Detector Systems (EGADS) facility [495] and the Super-Kamiokande Collaboration has 
just approved the future enrichment of Super-Kamiokande with Gd, the so-called SuperK- 
Gd project [496]. This would result in a reduction of the background by a factor of ^ 5 
for invisible muons and by at least an order of magnitude for spallation [184,185,207,246]. 
The planned liquid scintillator detector, JUNO, should also be able to separate DSNB 
neutrino events from the background very efficiently, offering a great opportunity to 


( 14 ) Note that solar neutrinos might constitute a non-minor impediment to the DSNB detection 
in LAr detectors possibly due to uncertainties in the cross section and poor pointing. 
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Fig. 45. - Smeared event rates as a function of detected energy based on an optimistic estimation 
of the DSNB obtained assuming a 27 M© SN progenitor (see Sec. 2*4) as representative of the 
SN population and MSW effects only. The rates are shown for a 560-kton WC detector (top, 
blue curves), a 20-kton liquid scintillator (center, green curves) and a 40-kton LAr detector 
(bottom, red curves). MSW oscillations under NH and IH assumptions are included. 


detect the DSNB [223,224,497]. 

In order to give an idea of the expected number of events for the detectors described 
above (see also Sec. 3), Fig. 45 shows the smeared event rates as a function of the detected 
energy for a 560-kton WC detector (blue curve), a 20-kton liquid scintillator (green curve) 
and a 40-kton LAr detector (red curve). The rates have been computed including the 
dominant reaction channels (i.e., IBD in water and scintillator detectors, and v e CC 
interactions on 40 Ar for the LAr detector) and assuming the 27 M 0 SN progenitor 
introduced in Sec. 2‘4 as representative of the whole SN population; for simplicity, v-v 
interactions have been neglected. Note as the presented rates may be optimistic because 
of the high NS mass of the adopted 27 M 0 SN progenitor. 

5‘5. Outlook. - The study of the DSNB is still affected by several theoretical uncer¬ 
tainties, preventing us from a precise forecasting of the expected signal. Assuming that 
synoptic surveys will be able to pin down the uncertainty on the SNR up to 5% and that 
the neutrino mass hierarchy will be known within the next decade, what could we learn 
from the DSNB detection? 

• Constraints on the stellar population. The DSNB receives contributions from all 
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families of SN progenitors and it will be an independent test of the global SN rate. 

• Constraints on the fraction of core-collapse and failed SNe. The DSNB will be 
a precious instrument to test our current knowledge of the core-collapse physics, 
especially for the failed SN class only detectable through its neutrinos. Constraints 
on the NS/BH formation ratio may also indirectly help to decipher the still disputed 
mechanism (s) by which stars explode as super novae. 

• Constraints on the neutrino emission properties. As discussed in Ref. [498], with the 
existing current upper limits on the v e diffuse background from Super-Kamiokande 
and the few events from the SN 1987A, it is possible to exclude some regions of the 
Lp e -(E„ e ) parameter space. Such bounds will surely be improved by real data. 

• A precise estimation of the DSNB is also relevant for the background modeling in 
direct dark matter searches [499,500]. 

• The DSNB detection could be useful to constrain neutrino decay models [501,502], 
neutrino electric or magnetic transition moments [503] and the existence of light 
scalar bosons [504]. 
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6. — Conclusions and Perspectives 

Neutrinos play a crucial role during all stages of stellar collapse and explosion. The 
emission of electron neutrinos produced by electron captures accelerates the implosion 
of the unstable, degenerate stellar core and mediates the neutronization of its matter. 
The absorption of electron neutrinos and antineutrinos above the neutrinospheric lay¬ 
ers initiates and powers the blast wave of the explosion and determines the trans-iron 
nucleosynthesis by setting the neutron-to-proton ratio in the innermost, neutrino-heated 
ejecta. The emission of neutrinos and antineutrinos of all flavors drives the cooling evo¬ 
lution of the newly formed neutron star that is left behind at the center of the supernova 
explosion. In order to summarize the expected neutrino emission characteristics, we 
show in Fig. 46 the source luminosities L v (upper panels), average energies (E u ) (middle 
panels), and spectral shape parameters a (lower panels) of z/ e , v e and v x for the shock- 
breakout, accretion and cooling phases of a 27 M 0 supernova progenitor simulated by 
the Garching group. 

The detection of two dozen electron antineutrinos from SN 1987A in three under¬ 
ground experiments splendidly confirmed our basic theoretical picture of stellar collapse 
and the birth of neutron stars. However, the signal statistics were too poor to yield 
detailed information of the explosion mechanism, which is still one of the most nagging 
problems in stellar astrophysics. The solution of this riddle is of fundamental importance 
for a better understanding of the origin of neutron stars and black holes and for the defi¬ 
nition of the role of supernovae in the cosmic cycle of element formation. In this context 
only neutrinos (and gravitational waves) can serve as direct probes of the physics taking 
place at the center of the explosion. A well resolved time and energy dependent neutrino 
signal from a future Galactic event will therefore provide a benchmark of unprecedented 
value for supernova physics. 

Hydrodynamical supernova simulations have now reached the multi-dimensional front 
including a full and sophisticated treatment of the neutrino reactions that play a role in 
the stellar core. The first three-dimensional supernova models of that sort have shown 
that neutrinos carry imprints of hydrodynamic processes occurring around the newly 
formed neutron star during the accretion phase. In particular, quasi-periodic modulations 
of the neutrino emission reflect dynamical mass motions that precede and enable the onset 
of the explosion. Rapid declines of the neutrino luminosity and non-monotonicities in 
the evolution of the spectral parameters may reveal the composition-shell structure of 
the progenitor and can indicate the end of the accretion phase (see Fig. 46). The spectral 
evolution of the neutrino signal is determined by the neutrinospheric conditions and thus 
carries radius and mass information of the contracting proto-neutron star. The late-time 
evolution of the neutrinos radiated during the cooling phase has a great potential to 
yield important insights into the still incompletely understood properties of matter at 
supranuclear densities in the interior of the nascent neutron star. 

While multi-dimensional supernova modeling seems to be on a good track to confirm 
the viability of the neutrino-driven explosion mechanism within the present uncertain¬ 
ties of input physics and initial conditions, an empirical validation of the theoretical 
concept is still missing. A Super-Kamiokande-size detector for electron neutrinos like 
DUNE, in combination with existing (Super-Kamiokande, IceCube, etc.) and proposed 
(Hyper-Kamiokande, JUNO) water and liquid scintillator detectors for electron antineu¬ 
trinos, would provide excellent time and/or spectral resolution for a Galactic supernova. 
Therefore, they would offer a fundamentally new quality for supernova research through 
neutrino measurements. Resolving the prompt shock-breakout burst of electron neutri- 
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Fig. 46. - Post-bounce evolution of the neutrino signal from a neutron star with baryonic (final 
gravitational) mass of ~1.776 (~1.592) M® formed in the explosion of a 27 Mq progenitor. The 
panels show the luminosities (top, in units of Bs -1 = 10 51 ergs _1 ), mean energies (ratios of 
energy fluxes to number fluxes; middle ), and spectral shape parameters ( bottom ) for v e (solid 
lines), v e (dashed), (dotted) and (dash-dotted) during the shock-breakout burst of v e 
(left column), accretion phase (middle column), and proto-neutron star cooling phase (right col¬ 
umn). The ID simulation was performed with the nuclear EoS of Lattimer & Swesty [90] (using 
the incompressibility modulus of K — 220 MeV) and took into account self-energy shifts in the 
[3 -interactions of free nucleons and a mixing-length treatment of proto-neutron star convection. 
The steep decline of the v e and v e luminosities at ^0.2 s is associated with the decrease of the 
mass-accretion rate when the edge of the stellar silicon layer falls through the stalled shock. The 
explosion was artificially initiated at 0.5 s and accretion finally ends at ~0.8s after bounce. 


nos would allow to determine the distance to the supernova even without electromagnetic 
detection and would define a precise reference point for the instant of core bounce. The si¬ 
multaneous information in both electron neutrino and antineutrino sectors could set con¬ 
straints on the equation of state (e.g., the baryonic symmetry energy) at neutrinospheric 
conditions and, in particular, it would help to answer the long-standing question whether 
supernovae are sources of r-process elements or of certain proton-rich isotopes produced 
by the neutrino-proton process. The direct comparison of electron neutrino and antineu¬ 
trino measurements could also serve to clarify the existence of a dipolar lepton-emission 
asymmetry (LESA) that is a recent, unexpected (though yet unconfirmed) discovery by 
the first three-dimensional supernova models with important implications for supernova 
nucleosynthesis and pulsar kicks. Moreover, a multitude of complex and not yet fully 
understood flavor transformation effects can modify the neutrino signal on the way from 
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the source at the supernova center to the detectors on Earth. Supernova neutrino data 
would provide extremely valuable information to decipher the evolution of the neutrino 
flavor field, which is indispensable for disentangling the source properties from effects 
imprinted on the signal by flavor transformations. 

The role of astrophysical messengers played by supernova neutrinos is largely associ¬ 
ated with the signatures imprinted on the observable neutrino burst by the supernova 
core dynamics and by the flavor conversions occurring deep inside the star. Within a ra¬ 
dius of a few hundred kilometers from the neutrinosphere, the neutrino field is so dense to 
become a “background to itself,” making the neutrino flavor evolution highly non-linear 
and leading to surprising and counterintuitive collective phenomena, when the entire neu¬ 
trino system oscillates coherently as a single collective mode. The rich phenomenology 
associated with these non-linear flavor dynamics is still in its infancy and many unex¬ 
pected results have been found in the past years. It is therefore mandatory to further 
investigate the role of flavor instabilities in the interacting neutrino field. Directions of 
further studies include: The characterization of the self-induced oscillations in models 
with reduced symmetries compared to what is usually assumed in the neutrino light-bulb 
model; investigations of the role of flavor conversions in the light of multi-dimensional 
hydrodynamic supernova simulations and their impact on the explosion dynamics and 
nucleosynthesis; or the role of residual scatterings and the possibility of self-induced 
spin-flavor transitions. 

Matter effects on supernova neutrino oscillations could imprint peculiar signatures 
on the observable neutrino signal with a strong sensitivity to the still-unknown neutrino 
mass hierarchy. Neutrino flavor conversions in the wake of the shock-wave propagation 
could also represent an intriguing possibility to follow the supernova dynamics in real time 
through measurements of the neutrino signal. The impact of the matter turbulence on the 
shock-wave signature would require further investigations. Studies in these directions will 
be mandatory once future, high-resolution supernova simulations will become available. 

Current and planned large underground neutrino detectors offer unprecedented op¬ 
portunities to study supernova and neutrino properties through high-statistics signals 
from different interaction channels. However, the exciting possibilities associated with 
the next Galactic explosion are in tension with the low Galactic supernova rate. The su¬ 
pernova rate in our Galaxy and within the Magellanic Cloud is estimated to be between 
one and several per century, but even with a high assumed rate we might have to wait a 
long time for the next Galactic supernova. In contrast, the DSNB, i.e., the flux of neu¬ 
trinos and antineutrinos coming from all past core-collapse supernovae in our Universe, 
is a guaranteed signal to study supernovae by their neutrinos in the forseeable future. 
The perspectives for a detection of the DSNB within the next decade are promising, 
especially with the just approved SuperK-Gd project (i.e., gadolinium-enhanced Super- 
Kamiokande detector) and with JUNO as a planned liquid scintillator experiment. Such 
a measurement would push the frontier of neutrino astronomy literally towards the edges 
of our Universe. 

In conclusion, the detection of supernova neutrinos represents the next frontier of low- 
energy neutrino astrophysics. Supernovae are celestial laboratories where neutrinos, after 
escaping from the highly opaque core, play a crucial role in the mechanism of the stellar 
explosion and where they might allow us to gain new insights into fundamental neutrino 
properties. In spite of remarkable progress, a satisfactory understanding of the formation 
of the neutrino emission as well as of the flavor conversions in this exotic environment is 
still lacking. Further theoretical and experimental work is therefore needed in order to 
get prepared for exploiting the wealth of information that the next Galactic supernova 
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explosion is going to provide. 
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Acronyms 


Table V. - List of the most used acronyms. 


ID 

One-dimensional 

2D 

T wo-dimensional 

3D 

Three-dimensional 

BH 

Black Hole 

CC 

Charged Current 

CEz/NS 

Coherent Elastic Neutrino-Nucleus Scattering 

DSNB 

Diffuse Supernova Neutrino Background 

EoM 

Equations of Motion 

EoS 

Equation of State 

ES 

Elastic scattering 

IBD 

Inverse Beta Decay 

IMB 

Irvine-Michigan-Brookhaven 

IMF 

Initial Mass Function 

IH 

Inverted Hierarchy 

LAr TPC 

Liquid Argon Time Projection Chamber 

LESA 

Lepton Emission Self-sustained Asymmetry 

l.h.s. 

left-hand-side 

LSD 

Large Scintillator Detector 

LVD 

Large Volume Detector 

LS EoS 

Lattimer-Swesty Equation of State 

MAA 

Multi-azimuthal-angle 

MSW 

Mikheyev-Smirnov-Wolfenstein 

NC 

Neutral Current 

NH 

Normal Hierarchy 

NS 

Neutron Star 

r.h.s. 

right-hand-side 

S EoS 

Shen et al. Equation of State 

SASI 

Standing Accretion Shock Instability 

SFR 

Star Formation Rate 

SN 

Supernova 

SNe 

Super novae 

SNEWS 

SuperNova Early Warning System 

SNR 

Supernova Rate 

WC 

Water Cherenkov 
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